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II  ABSTRACT 


The  cyclotron  resonance  ina-  uon  between  ULF-VLF  waves  and  trapped  particles 
may  be  strongly  enhanced  in  the  magnetosphere  by  the  injection  of  either  cold 
plasma  or  energetic  particle  beams.  A  variety  of  natural  conditions  and  injection 
are  investigated  numerically  using  the  linear  theory.  Both  satellite  and  rocket 
injection  of  cold  plasma  provide  significant  enhancements  of  amplification. 
Conventional  hot  electron  beams  may  also  amplify  narrow  frequency  bands,  if  the 
beam  geometry  can  be  effectively  altered  at  injection.  Heavy  ion  beams  appear 
to  offer  attractive  propagation  conditions  as  well  as  stimulate  amplification. 
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by 

Harold  B.  Liemohn 

SUMMARY 

The  geomagnetic  cavity  surrounding  the  earth,  called  the  magneto¬ 
sphere,  contains  an  enormous  amount  of  energy  in  the  form  of  electro¬ 
magnetic  waves,  trapped  energetic  particles,  and  the  static  magnetic 
field.  The  waves  and  particles  are  in  a  state  of  dynamic  equilibrium 
as  they  interact  with  each  other  through  various  plasma  instabilities. 

Our  understanding  of  this  equilibrium  state  has  progressed  to  the  point 
where  it  may  be  possible  to  stimulate  wave  amplification  at  ULF  (^0.1-10  Hz) 
and  VLF  M-100  kHz)  by  disturbing  the  equilibrium  through  catalytic 
injections  of  plasma  clouds  or  particle  beams.  Such  a  prospect  has  great 
significance.  First,  it  offers  a  direct  test  of  the  fundamental  physical 
processes  that  are  believed  to  control  the  state  of  the  magnetosphere. 
Second,  controlled  injections  might  be  used  to  stimulate  artificial 
amplification  of  VLF  nr  ULF  emissions  on  command  and  allow  them  to  be  used 
as  a  wide-area  communication  system. 

One  of  the  dominant  plasma  instabilities  in  the  magnetosphere  is  the 
cyclotron-resonance  interaction  between  energetic  particles  and  electro¬ 
magnetic  waves.  The  condition  for  the  interaction  to  occur  is  that  the 
motion  of  the  particle  Doppler  shifts  the  wave  frequency  to  its  local 
cyclotron  frequency.  The  net  exchange  of  energy  between  a  band  of  waves 
and  a  group  of  particles  depends  upon  the  shape  of  the  phase-space  dis¬ 
tribution  of  the  particles  and  the  local  frequency  parameters  of  the  plasma 
medium.  The  interaction  is  undoubtedly  responsible  for  significant  wave 
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amplification  and  particle  precipitation  in  certain  circumstances  and 
appears  to  be  the  process  responsible  for  maintenance  of  an  equilibrium 
configuration  of  waves  and  particles  in  the  magnetosphere. 

These  qualities  make  this  interaction  a  logical  first  choice  for 
investigation  as  a  means  of  stimulating  artifical  ampl it ication  of  VLF 
and  ULF  waves  in  the  magnetosphere.  The  interaction  is  sensitive  to 
both  the  background  magnetoplasma  parameters  and  the  hot  energetic  particle 
distribution.  Modification  of  the  interaction  can  be  achieved  by  the 
introduction  of  localized  clouds  of  plasma  that  reduce  the  phase  velocity 
of  the  waves  and  cause  enhanced  amplification  by  lowering  the  resonance 
velocity.  Alternatively,  direct  injection  of  hot  energetic  particle  beams 
can  also  appreciably  alter  the  local  amplification  characteristics  of  the 
beam.  These  two  methods  of  modifying  the  cyclotron-resonance  interaction 
are  explore d  in  this  document  for  a  variety  of  cases. 

Many  different  methods  for  injecting  plasma  clouds  or  particle  beams 
nave  been  modelled  quantitatively  to  assess  their  relative  merits  as 
a_means  of  inducing  strong  wave  amplification.  The  injection  of  a  dense 
jet  of  barium  that  travels  upward  along  the  geomagnetic  field  toward  the 
equatorial  region  produces  significant  amounts  of  enhanced  VLF  amplifica¬ 
tion  when  the  jet  is  within  30°  of  the  geomagnetic  equator  some  thirty 
i  nutes  after  injection.  Injection  of  a  geosynchronous  lithium  cloud 
stimulates  amplification  of  both  VLF  and  ULF  waves,  but  the  latter  are 
limited  to  frequencies  below  the  equatorial  lithium  gyrofrequency 
M.|  Hz).  The  extent  of  the  lithium  amplification  effect  depends  criti¬ 
cally  on  the  the  natural  conditions  in  the  magnetosphere,  particularly  the 
intensity  of  the  proton  ring  current.  Amplifications  of  many  tens  of 
decibels  appear  to  be  feasible  with  injections  of  only  1-2  kgm  of  material. 
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Beam  injections  were  included  in  the  study  because  they  offer  the 
opportunity  for  prompt  localized  interactions  with  the  waves  that  can 
produce  readily  identifiable  signatures  of  the  amplified  signal.  Electron 
beams  amplify  VLF  whereas  proton  beams  amplify  ULF  signals.  Unfortunately, 
the  beam  diameter  from  a  conventional  particle  gun  on  board  a  satellite 
or  rocket  is  too  small  for  effective  amplification.  An  adequate  interaction 
of  the  beam  with  a  wave  front  requires  a  cross  section  comparable  to  several 
wave  lengths.  Nevertheless,  if  such  a  broad  beam  is  achievable  in  the 
near  future,  appreciable  amounts  of  amplification  are  possible  and  their 
natural  frequency  signature  is  quite  distinct.  Furthermore,  modulation 
of  such  beams  offers  the  opportunity  to  transmit  information  instead  of 
merely  enhancing  background  noise. 

Injection  of  low-energy  beams  of  heavy  ions  appears  to  be  practical  for 
stimulation  of  VLF  amplification.  For  example,  a  beam  of  43  eV  cesium 
ions  has  a  gyroradius  that  is  comparable  to  a  VLF  wave  front.  When  such 
a  beam  travels  up  a  field  line  it  drags  along  a  neutralizing  low-energy 
electron  beam  which  changes  the  local  propagation  conditions.  The  ampli¬ 
fication  of  VLF  by  such  a  beam  appears  to  be  on  the  order  of  tens  of 
decibels  at  frequencies  near  trie  local  electron  gyro  frequency.  Again 
the  characteristic  frequency  signature  is  present  and  beam  modulation 
can  further  enhance  t'  9  detectability  of  induced  amplification. 

Several  conclusions  can  be  drawn  from  the  study.  Amplifications 
of  40-60  decibels  are  indicated  for  several  cases  of  interest.  Such 


enhancements  of  VLF  or  ULF  noise  levels  are  easily  detectable  with 
ground  or  satellite  receivers  located  along  the  injection  field  lines. 

One  of  the  most  promising  cases  is  lithium  injection  at  geosynchronous 
altitude  where  broad  regions  at  the  surface  of  the  earth  would  be 
illuminated  by  ULF  noise  in  the  band  0.4  -  0.5  Hz.  At  VLF  the  low- 
energy  cesium-electron  beam  works  well  as  an  amplifier  and  also  pro¬ 
vides  a  fluxtube  waveguide  for  the  electromagnetic  energy.  Although 
the  electron  and  proton  beams  appear  to  generate  amplification,  the 
inte  action  region  is  severely  limited  with  present  particle  gun  tech¬ 
nology  and  the  total  energy  transferred  to  the  waves  is  quite  small. 
Plasma  injections  into  the  more  dense  inner  magnetosphere  are  relatively 
ineffective  due  to  the  local  propagation  characteristics. 

The  study  is  based  on  the  linear  theory  for  the  cyclotron-resonance 
interaction.  Since  some  of  the  linearity  assumptions  are  not  valid  in  the 
injection  region,  the  quantitative  results  are  open  to  question.  Never¬ 
theless,  the  linear  results  are  believed  to  provide  a  reasonable  basis 
for  selecting  injection  methods  that  appear  to  be  effective  for  stimu¬ 
lating  amplification.  As  might  be  anticipated,  the  application  of  the 
nonlinear  theory  is  much  more  involved  and  requires  extensive  computer 

time.  Nevertheless,  a  thorough  nonlinear  study  of  promising  cases  is 
recommended. 
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Pitch-angle  distribution  models  H0  at  the  geomagnetic  equator.  22 
These  anisotropic  model s  .including  a  loss  cone  cutoff,  are  used 
for  both  electrons  and  protons  in  the  magnetosphere. 

Typical  electron-energy  distributions  in  the  magnetosphere.  23 

Both  the  cold  background  plasma  and  the  hot  energetic 

particles  are  illustrated  to  show  their  relative  magnitudes 

and  locations  on  the  energy  scale.  Some  of  the  calculations 

are  performed  for  energy  distributions  of  the  form  E"n  that 

are  typical  of  many  observations.  The  substorm  distribution 

'A'  was  measured  on  the  ATS-5  satellite  at  geosynchronous 

orbit  (DeForest  and  Mcllwain,  1971). 

Increased  amplification  exponent  for  VLF  waves  due  to  a  26 

natural  enhancement  in  the  cold  plasma  density.  The  propaga¬ 
tion  at  L  =  4  is  assumed  to  be  Inside  the  plasmasphere  where 
the  cold  plasma  Is  in  hydrostatic  equilibrium.  For  this 
case,  the  hot  plasma  is  assumed  to  have  an  energy  and  pitch 
angle  distribution  of  the  form  *5  sin  a  (see  Figs.  1  and  2). 

Increased  net  path  amplification  of  VLF  waves  due  to  a  natural  27 
enhancement  in  the  cold  plasma  density.  The  propagation  path 
at  L  *  4  is  assumed  to  be  Inside  the  plasmasphere  where  the 
cold  plasma  is  in  hydrostatic  equilibrium.  For  this  case,  the 
hot  plasma  is  assumed  to  have  an  energy  and  pitch  angle  distri¬ 
bution  of  the  form  E~1 * 5  sin  a  (see  Figs.  1  and  2). 

Rocket  injection  of  a  barium  shaped-charge  along  the  L  =  4  30 

field  line.  The  barium  jet  is  assumed  to  consist  of  1  kgm 
of  fully  ionized  gas.  The  trajectory  is  assumed  to  be  in 
the  pjasmatrough,  where  N  =  1  electron/cm_3  and  N  varies 
as  R-3. 

Amplification  exponent  for  VLF  waves  propagating  through  the  32 
barium  shaped-charge  injection  along  l  =  4.  The  hot  electrons 
are  assumed  to  have  the  distribution  E_2  sin^a  (see  Figs.  1 
and  2). 

Net  path  amplification  for  VLF  waves  propagating  through  the  33 
barium  shaped-charge  Injection  along  L  =  4.  The  hot  electrons 
are  assumed  to  have  the  distribution  E”2  sin2a  (see  Figs.  1 
and  2). 

Injection  of  lithium  from  a  point  source  (cannister)  in  36 

geosynchronous  orbit  (L  =  6.6).  One  kgm  of  total  ionization 
is  assumed  to  be  diffusing  spherically  from  the  source  with 
a  diffusion  coefficient  D  -v  1  Km'  /s.  The  ambient  background 
plasma  is  assumed  to  be  in  the  plasmatrough  where  N0  •  1  cm- 3. 
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Amplification  exponent  for  VLF  waves  propagating  through  the 
geosynchronous  lithium  injection.  The  exponent  is  sharply 
peaked  In  the  equatorial  enhancement  which  has  a  scan  diameter 
of  about  6400  km  at  50  minutes  after  Injection.  The  hot  electrons 
are  assumed  to  have  the  distribution  E’T,5s1n  a. 

Net  path  amplification  for  VLF  waves  propagating  through  a 
geosynchronous  lithium  injection.  The  enhancement  Is  entirely 
attributable  to  the  lithium  bubble.  The  hot  electrons  are 
assumed  to  have  the  distribution  F-l-5  sin 

Amplification  exponent  for  VLF  waves  propagating  through  a 
geosynchronous  lithium  bubble.  The  hot  electron  distribution 
Is  based  on  Model  'A'  for  a  typical  geomagnetic  substorm 
(see  Fig.  2). 


Net  path  amplification  of  VLF  waves  propagating  through  a 
geosynchronous  lithium  bubble.  The  hot  electron  distribution 
Is  based  on  Model  'A'  for  a  typical  geomagnetic  substorm. 

Electron  beam  model  distributions  for  pitch  angle,  Hq,  and 
differential  energy  spectrum,  dN/dE.  The  beam  distributions 
are  sharply  peaked  In  pitch  angle  and  energy  assuming  a 
particle  gun  deploys  the  beam  from  a  satellite  or  very  hlqh 
altitude  rocket. 

Amplification  exponent  for  VLF  waves  propagating  through  an 
electron  beam  along  the  field  line  at  L*  4  inside  the 
plasmapause.  The  electron  beam  parameters  are  an  equatorial 
pitch  angle  of  55°  and  a  particle  beam  energy  of  10  keV.  See 
the  text  for  an  explanation  of  the  shape  of  these  curves. 

Net  path  amplification  for  VLF  waves  propagating  through  an 
electron  beam  at  L  ■  4  Inside  the  plasmapause.  Equatorial 
pitch  angles  for  the  beam  are  25°,  55°,  and  85°  and  the  beam 
energy  Is  10  keV.  See  the  text  for  an  explanation  of  the 
shape  of  these  curves. 


Schematic  diagram  of  magnetic-gravitational  trapping  of  low 
energy  Ions  on  an  Intra hemisphere  trajectory.  Injection  of 
these  low  energy  Ions  above  the  Ionosphere  can  produce  a 
sharp  density  enhancement  along  a  portion  of  the  flux  tube. 

Trapping  zones  for  warm  Ion  Injection  above  the  Ionosphere. 
The  end  points  of  particle  orbits  for  the  various  Ions  are 
displayed  as  a  function  of  energy.  The  Injection  level  Is 
at  an  altitude  of  1600  km  along  the  field  line  at  L  ■  4. 
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at  L  *  3  5.  The  neutral  lithium  jet  of  approximately  1  kgm  is 
expected  to  follow  a  ballistic  trajectory  across  the  magnetic 
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lithium  shaped-charge  Injection  along  L  ■  3.5.  The  hot  protons 

are  assumed  to  have  the  distribution  E"2  sin^a  (see  Figs.  1 
and  21). 
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diameters  of  1280  and  6400-km.  The  hot  protons  are  assumed  to 
have  the  distribution  E'1**5  sin  a. 

26  Net  path  amplification  for  ULF  waves  propagating  through  a  72 

geosynchronous  lithium  Injection  (see  Fig.  8).  The  enhancement 

Is  entirely  attributable  to  the  lithium  bubble.  The  hot  protons 
are  assumed  to  have  the  distribution  E~ 1  * 5  sin  a. 
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Amplification  exponent  for  ULF  waves  propagating  through  a  73 

geosynchronous  lithium  Injection  (see  Figure  8).  The  ex¬ 
ponent  Is  sharply  peaked  In  the  equatorial  enhancements  which  have 
mean  diameters  of  1280  and  6400  km.  The  hot  proton  distri¬ 
bution  Is  based  on  model  'B'  for  the  ring  current  of  a  typical 
geomagnetic  substorm  (see  Figure  21). 

Net  path  amplification  for  ULF  waves  propagating  through  a  74 

geosynchronous  lithium  Injection.  The  enhancement  is  entirely 
attributable  to  the  lithium  bubble.  The  hot  proton  distribution 
is  based  on  model  B‘  for  the  ring  current  of  a  typical  geo¬ 
magnetic  substorm  (see  Figure  21). 

Amplification  exponent  for  ULF  waves  propagating  through  a  proton  79 
beam  along  the  field  line  at  L=4  Inside  the  plasmapause.  The 
beam  model  has  the  same  distributions  as  for  electrons  (see 
Figure  13).  The  beam  has  an  equatorial  pitch  angle  of  55° 
and  a  particle  beam  energy  of  10  keV. 

Net  path  amplification  for  ULF  waves  propagating  through  a  proton  80 
beam  along  the  field  line  at  L*4  Inside  the  plasmapause. 

Equatorial  pitch  angles  for  the  beam  are  25°,  55°,  and  85°  and 
the  beam  energy  Is  10  keV. 
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>  STIMULATED  AMPLIFICATION  OF  VLF  AND  ULF  WAVES 

IN  THE  MAGNETOSPHERE  BY  LOCALIZED  INJECTIONS 
OF  PLASMA  CLOUDS  AND  PARTICLE  BEAMS 

,  INTRODUCTION 

There  is  an  ongoing  military  and  civilian  requirement  to  improve 
the  reliability  of  global  radio  communications.  Long  distance  transmis¬ 
sion  of  information  is  interrupted  frequently  by  the  diurnal  excursions 
of  ionospheric  propagation  conditions  and  occasionally  by  strong  geo¬ 
magnetic  disturbances  that  cause  exceptionally  intense  ionization.  To 
some  extent  these  vagaries  are  overcome  at  VLF  (3-30  kHz)  where  the  long 
wavelengths  are  less  affected  by  changing  conditions  in  the  ionosphere. 

At  lower  frequencies,  notably  ELF  (3-3000  Hz)  and  ULF  (below  3  Hz),  the 
signals  readily  propagate  on  a  global  scale  with  relatively  negligible 
attenuation.  At  these  latter  frequencies,  of  course,  the  information 
band  width  is  severely  limited.  Furthermore,  efficient  generation  of 
such  signals  is  exceedingly  difficult  due  to  the  extremely  long  free- 
space  wavelengths. 

Since  VLF,  ELF,  and  ULF  noise  bands  occur  naturally  in  the  mag¬ 
netosphere  surrounding  earth,  it  is  appropriate  to  inquire  about  v/ays  of 
exciting  this  energy  on  command.  The  electromagnetic  waves  are  apparently 
quided  by  the  geomagnetic  field  to  ionospheric  alitudes  where  they  subse¬ 
quently  illuminate  wide  surface  areas  around  the  base  of  the  field  line. 
In  the  maqnetosphere,  VLF  and  ULF  noise  is  believed  to  be  caused  by  the 
cyclotron-resonance  interation  with.. energetic  (1-100,  keV)  electrons  and 
protons.  In  fact,  the  dynamic  equilibrium  state  between  waves  and 


energetic  particles  in  the  magnetosphere  is  believed  to  be  controlled 
by  this  dominant  Interaction  (Kennel  and  Petschek,  1966).  Furthermore, 
observed  amplification  of  VLF  and  ULF  whistlers  is  explained  by  this  inter¬ 
action  (e.g.f  Cornwall,  1966;  and  Liemohn,  1967).  The  amount  of  energy 
exchange  that  may  occur  depends  critically  on  such  parameters  as  the 
local  ambient  plasma  density,  the  local  geomagnetic  field  strength,  and 
the  shape  of  the  energetic  particle  distribution.  Thus,  it  is  appropriate 
to  explore  ways  of  altering  some  of  these  parameters  which  might  catalytically 
enhance  the  amplification  process.  Such  experiments  would  be  valuable  as 
a  quantitative  test  of  the  cyclotron-resonance  interaction  theory  and  subse¬ 
quently  provide  a  basis  for  an  operational  communications  system. 

The  physics  of  the  cyclotron-resonance  interaction  is  really 
simple  to  understand,  although  its  mathematical  description  in  magneto- 
spheric  applications  is  relatively  complicated.  In  essence,  it  is  a  simple 
resonance  trapping  of  selected  charged  particles  in  the  electromagnetic 
potential  of  the  wave.  Although  all  particles  are  perturbed  by  the  circu¬ 
larly  polarized  wave  fields,  there  is  a  select  group  of  particles  with 
velocities  that  Doppler  shift  the  propagation  frequency  to  the  vicinity  of 
their  local  gyrofrequency.  The  electromagnetic  forces  slow  down  some  of  these 
particles  and  speed  up  others  to  bring  them  closer  to  the  resonance  fre¬ 
quency  and  a  net  exchange  of  energy  occurs  between  each  monochromatic  wave 
and  the  relevant  particles.  If  more  particles  give  up  energy  than  gain,  the 
wave  is  amplified  and  otherwise,  attenuated.  The  amount  of  energy  ex¬ 
change  depends  on  the  shape  of  the  particle  distribution  in  the  vicinity 
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of  the  resonance  velocity  band.  The  local  parameters  of  the  medium  which 
prescribe  the  phase  velocity  of  the  wave  also  influence  the  extent  of  the 
interaction  by  establishing  minimum  resonance  energies  for  the  particles. 

In  order  to  estimate  its  effect  in  the  magnetosphere,  the  inter¬ 
action  must  be  evaluated  for  a  n/ide  range  of  propagation  frequencies, 
field-line  paths,  and  particle  distributions.  Several  cases  that  are 
typical  of  magnetospheric  conditions  have  been  investigated  (e.g.,  Cocke  and 
Cornwall,  1967;  Liemohn,  1967  and  Ho  and  Liemohn,  1972)  for  a  wide  range 
of  conditions.  It  was  concluded  that  the  magnetosphere  behaves  like  a 
c  jnvective  traveling-wave  amplifier  at  both  VLF  and  ULF,  and  that  observable 
amounts  of  VLF  and  ULF  amplification  occurred  over  a  wide  range  of  fre¬ 
quencies  but  varied  appreciably  with  geomagnetic  activity.  A  considerable 
body  of  indirect  evidence  has  been  gathered  in  support  of  the  cyclotron- 
resonance  amplification  theory  (Cocke  and  Cornwall,  1967;  Liemohn,  i969; 
Cornwall  et  al . ,  1971;  Thorne  and  Kennel,  1971;  Rosenberg  et  al . ,  1971). 
Unfortunately,  there  have  been  no  direct  experimental  observations  to  verify 
the  quantitative  aspects  of  the  cyclotro  resonance  theory. 

Thus  the  possibility  of  a  controlled  experiment  to  enhance  VLF 
or  ULF  amplification  by  a  known  amount  has  a  two-fold  purpose.  First, 
it  can  be  used  to  test  the  quantitative  conclusions  that  have  been  drawn 
from  the  theoretical  analyses.  Second,  successful  artificial  simulation 
of  wave  amplification  in  the  magnetosphere  is  an  essential  step  toward 
its  ulilization  in  communication  systems. 

There  are  essentially  two  known  methods  for  artificial  stimulation  of 
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the  cyclotron-resonance  interaction  in  the  magnetosphere.  The  first 
method  involves  injection  of  dense  plasma  clouds  of  low-energy  electrons 
or  ions  (Brice,  1970;  and  Cornwall  &  Schulz,  1971)  into  regions  of  the 
magnetosphere  where  the  ambient  plasma  density  is  low  but  energetic 
particles  abound.  Selective  placement  of  these  clouds  reduces  the  phase 
velocity  of  the  VLF  or  ULF  waves  and  thereby  reduces  the  minimum  particle 
energy  necessary  for  resonance.  Generally  there  are  many  more  hot  particles 
at  the  reduced  resonance  energy  so  that  appreciable  amplification  enhance¬ 
ments  are  anticipated.  Since  the  plasma  clouds  do  not  enter  directly  into 
the  amplification,  their  injection  is  viewed  as  a  catalytic  process  whereby 
the  increased  wave  energy  is  derived  from  the  natural  energetic  particle 
population. 

The  second  method  of  artificially  stimulating  the  amplification 
involves  direct  injection  of  hot  plasma  beams  by  rocket  or  satellite 
particle  guns.  Beam  densities  comparable  to  the  ambient  natural  particle 
densities  along  the  flux  tube  are  apparently  achievable  with  current 
technology  (Hendrickson  et  al . ,  1970;  Hess  et  al.,  1971;  and  Cartwright 
&  Kellogg,  1971).  In  this  case  the  particle  beam  merely  adds  a  known 
distortion  in  the  energetic  particle  distribution  along  the  field  line, 
and  the  enhanced  wave  amplification  is  derived  entirely  from  the  energy 
of  the  injected  beam.  This  method  of  stimulation  is  severely  limited  by 
the  low  energy-conversion  efficiency  and  the  relatively  small  beam  diameter. 
Beams  are  attractive,  however,  because  they  do  not  rely  significantly  on 
the  ambient  conditions  in  the  magnetosphere,  and  their  injections  can  be 
modulated  to  transmit  information. 


In  this  document  selected  cases  of  plasma  cloud  and  beam  in¬ 
jections  are  analyzed  quantitatively  using  the  linear  theory  for  the 
cyclotron-resonance  interaction.  The  linear  theory  has  several  advan¬ 
tages  for  an  initial  investigation  but  it  suffers  certain  inadequacies 
that  necessitate  non  linear  investigations  to  corroborate  the  results. 

The  principle  advantage  of  the  linear  expressions  is  that  their  computer 
evaluation  along  magnetospheric  propagation  paths  is  relatively  inex¬ 
pensive.  Furthermore,  the  computer  program  was  generated  some  years 
ago  (Liemohn,  1967).  Computer  programs  for  the  complete  non  linear 
options  have  been  prepared  (Ossakow  et  al.,  1972;  Cuperman  and  Salu,  1972 
Denavit,  private  communication),  but  they  currently  require  about  100 
times  more  computer  time  for  each  uniform  element  of  the  magnetospheric 
path.  Thus,  the  linear  theory  provides  a  useful  means  of  exploring  a 
variety  of  parameter  regimes  to  optimize  promising  injection  schemes 
before  they  are  subjected  to  a  detailed  non  linear  investigation. 

The  next  section  of  the  paper  summarizes  the  linear  theory  for 
the  cy<  lotron-resonance  interaction  and  introduces  relevant  propagation 
properties  of  the  magnetosphere.  Specific  injection  schemes  for  stimu¬ 
lating  amplification  of  VLF  and  ULF  are  treated  separately  in  the  follow¬ 
ing  two  sections.  The  final  section  summarizes  the  important  results  of 
the  study  and  draws  specific  conclusions  about  future  research. 
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CYCLOTRON-RESONANCE  THEORY 

The  fundamental  linear  theory  for  the  cyclotron-resonance  inter¬ 
action  and  its  application  to  VLF  and  ULF  propagation  in  the  magneto¬ 
sphere  is  summarized  in  this  section.  Only  the  key  equations  are  presented 
here;  for  a  complete  derivation  of  the  theoretical  exDressions,  the  reader 
is  referred  to  the  early  literature  (e.g.,  Scarf,  1962;  Kennel  and  Petschek, 
1966).  The  solutions  of  the  dispersion  equation  for  circularly  polarized 
VLF  and  ULF  waves  lead  naturally  to  a  discussion  of  the  propagation 
characteristics  in  the  magnetosphere.  Some  subtle  aspects  of  the  matha- 
matical  solutions  are  described  in  order  to  alert  the  reader  to  the 
limitations  of  the  theoretical  expressions.  Finally,  models  of  the 
plasma  medium  encountered  in  the  magnetosphere  during  disturbed 
and  quiescent  conditions  are  described.  These  models  provide  a  quanti¬ 
tative  basis  for  estimating  the  effectiveness  of  various  methods  of 
stimulating  the  interaction. 

Dispersion  Equation.  In  order  to  allow  tractible  mathematical 
expressions,  several  assumptions  are  built  into  the  theoretical  descrip¬ 
tion  presented  here.  First,  it  is  assumed  that  the  wavelengths  are  suf¬ 
ficiently  small  that  the  propagation  can  be  described  in  a  localized 
region  where  the  magnetoplasma  has  a  uniform  static  field  and  a 
homogeneous  particle  density  N.  Second,  the  analysis  is  restricted  to 
the  special  case  where  the  propagation  vector  k  is  parallel  to  B  so  that 
the  waves  are  necessarily  circularly  polarized.  When  k  is  not  parallel 
to  B,  Landau  damping  is  introduced  by  the  presence  of  a  longitudinal 
electric  field  component  (Kennel  and  Thorne,  1967;  and  Kennel  and  Wong, 
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1967).  Third,  the  transfer  of  energy  between  the  waves  and  particles 
is  assumed  to  be  sufficiently  small  and  localized  that  non-linear  effects 
can  be  ignored.  Implications  of  these  specific  assumptions  will  be  dis¬ 
cussed  below. 


Circularly  polarized  monochromatic  waves  are  described  by  the  Fourier 
component  form  exp  (ikz-iwt)  where  ui  is  the  propagation  frequency  of 

interest.  For  the  foregoing  conditions,  the  dispersion  equation  for  the 
wave  characteristics  in  the  cyclotron-resonance  interaction  may  be  derived 
from  the  coupled  Vlasov  and  Maxwell  equations  (e.g.  Montgomery  and  Tidman, 

1 964) 


where 


and 


vc  =  (<d  +  wc)/k 


(1) 

(2) 

(3) 


1  0 

F  ’  ( Vj_  ,  v„ )  are  the  phase-space  distributions  of  the  charged  particles 

with  velocity  components  perpendicular  and  parallel  to  B.  The  quantities 
‘"pi’6  =  ( 4ttN ^  ’ e  e2/m  1,e)  f  and  u^1  ,e  =  (+  e)  B/m1  ,e  c  are  the  plasma 
and  cyclotron  frequencies  of  the  medium.  The  particle  labels  for  the  ions 
(i)  and  electrons (e)  are  generally  suppressed  throughout  the  remainder 
of  the  paper  to  simplify  notation. 

The  resonance  speed  vc  is  the  particle  velocity  component  necessary  U 
Doppler  shift  the  wave  frequency  to  the  local  cyclotron  frequency. 

Only  those  particles  with  Vy  near  vc  participate  in  the  energy  exchange. 

The  plus  and  minus  signs  in  vQ  refer  to  right-hand  and  left-hand 
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circularly  polarized  modes,  respectively.  Clearly,  there  are  four 
possible  combinations  of  interactions  between  waves  and  particles.  By 
convention  the  right-hand  electron  VLF  and  the  left-hand  ion  ULF  reson¬ 
ance  which  require  the  wave  and  resonance  particles  to  rotate  in  the  same 
sense  and  travel  in  opposite  directions  are  termed  normal  interactions. 

The  right-hand  ion  and  left-hand  electron  resonances,  which  require  the 
particles  to  overtake  the  wave,  thus  reversing  their  apparent  sense  of 
rotations  as  required  for  resonance,  are  called  anomalous.  In  the  analyses 
presented  subsequently,  only  the  normal  interactions  are  considered  be¬ 
cause  they  occur  at  much  lower  energies  (0.1-100  keV)  than  the  anomalous 
interactions  (above  1  meV).  At  anomalous  interaction  energies,  there  simply 
are  an  insufficient  number  of  particles  to  warrant  investigation.  Thus, 
there  is  only  one  plasma  component  that  contributes  to  the  summations  in 
equation  (1),  depending  on  the  mode  of  interest. 

The  dispersion  equation  (1)  determines  the  propagation  characteristics 
k(w)  or  oi(k)  for  the  electromagnetic  wave.  In  general  it  has  several 
solution  trajectories  in  the  complex  hyperspace  (ay  +  ica. ,  k^  +  i ki ) , 
depending  on  the  form  of  F.  The  space-time  conditions  in  a  given 
problem  determine  the  appropriate  locus  of  acceptable  solutions.  For  an 
initial -value  problem,  k  must  be  real  and  u>  is  allowed  to  be  complex, 
whereas  a  boundary-value  problem  requires  real  o>  and  complex  k.  Instab¬ 
ilities  inherent  in  the  system  are  prescribed  by  characteristic  roots 
of  equation  (1)  that  have  greater  than  0  or  k^  less  than  0  corresponding 
to  wave  amplitude  growth  in  the  representation  used  here  (with  ay,  kr 
greater  than  0).  When  such  roots  are  present,  they  normally  dominate  the 
solution.  In  some  applications,  however,  these  roots  may  be  suppressed  by 
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imposed  space-time  conditions  on  the  system. 

In  many  cases  of  physical  interest  such  as  those  presented  here, 
the  local  boundary  conditions  remain  unspecified,  and  the  dispersion 
equation  alone  is  the  only  information  available  about  the  nature  of  the 
wave-particle  system.  When  instabilities  occur  in  this  circumstance,  it 
is  mportant  to  ascertain  whether  an  unstable  system  is  Inherently  a  non- 
convective  runaway  oscillator  described  by  and  k  real,  or  a  convective 
traveling-wave  amplifier  described  by  a  real  and  k,  <  0.  The  answer  is 
found  in  a  sophisticated  mathematical  criterion  for  hot  plasma  instabilities 
(Derfler,  1967  and  1970).  The  criterion  has  been  applied  to  several  cases 
ot  interest  in  the  magnetosphere  (Ho  and  Liemohn,  1972)  and  it  has  been 
concluded  that  instabilities  in  the  magnetosphere  are  generally  convective. 
The  specific  cases  presented  below  have  not  been  investigated  on  the  basis 
of  the  instability  criterion,  but  they  generally  fall  into  the  class  of 
interactions  which  are  deemed  to  be  convective. 

Before  deriving  solutions  of  equation  (1)  for  complex  k  and  real 
u  ,s  necessar)'  to  define  the  form  of  the  distribution  function  F. 
Experimental  observations  of  particles  in  the  magnetosphere  reveal  that  F 
can  be  divided  into  a  very  low  energy  cold  plasma  part  FCQLD  for  particles 
below  10  eV,  and  a  hot  plasma  part  FH0J  which  describes  all  the  energetic 
particles  above  10  eV.  The  sum  of  these  two  distributions  is  normalized 
to  unity  so  that  the  total  plasma  density  N  is  explicitly  Included  in  „ 

In  order  to  solve  the  integral  dispersion  equation,  it  is  necessary 
to  make  an  additional  assumption  that  the  amplitude  changes  slowly  with 
wavelength, which  can  be  expressed  in  the  form  I  kf|  «kr.  Such  a  condition 
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usually  applies  in  the  quiescent  magnetosphere,  but  may  be  violated  during 
artificial  stimulation  of  the  interaction  by  plasma  injections.  When  the 
condition  applies,  the  imaginary  part  of  vc  is  very  small  and  equation  (1) 
may  be  expanded  in  a  Taylor  series  in  the  complex  v  plane  (Jackson, 

1960).  By  equating  real  and  Imaginary  parts,  the  propagation  characteristics 
to  lowest  order  in  /kr  are  given  by 


cVr 


U) 


ck, 


2c2k* 


V  v. 

uj  (w  u>c ) 

2  up  /  dvx 

1  ,e  y  o 


^HOT  ^Vj-’ 


(4) 


(5) 


In  the  latter  expression  only  F^  is  substituted  in  equation  (2)  to 
obtain  IH0T,  and  vc  is  the  real  part  of  equation  (3).  Expressions  (4) 
and  (5)  describe  the  local  phase  velocity  (refractive  Index  n)  and  amplitude 
exponent  for  VLF  or  ULF  waves  in  the  magnetosphere. 


In  some  applications  where  wave  growth  is  extremely  rapid,  I  k^ I  >  kr 
and  a  proper  description  requires,  the  nonlinear  theory  (Cuptrman,  1972; 
Bud'ko,  et  al.,  1972).  If  |k^|>>kr,  then  the  propagation  characteristics 
may  be  described  by  asymptotic  expansions  (e.g.,  Jackson,  1960)  of  the 
integral.  It  Is  Important  to  emphasize  that  these  analyses  are  not  the 
same  as  quasi-1 Inear  theory  (Kennel  and  Petschek,  1966;  Gendrin,  1968) 
where  repeated  interaction  with  the  same  particles  distorts  F^q-j..  Such 
distortion  is  only  Important  for  extended  interaction  regions,  where 
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1  6 

the  longitudinal  drift  in  one  bounce  period,  4  v/wc’  v ,  small  compared 

i  c. 

to  the  characteristic  wavelength,  2n/k  (evaluated  at  0.5  »  *  ).  For 

r  C- 

magnetospheric  conditions  of  Interest  here,  these  quantities  are  comparable 
for  10  keV  electrons  and  protons  (3  km  and  150  km,  respectively).  However, 
in  the  applications  only  narrow  meridian  wedges  several  wavelengths  across 
are  stimulated  by  injection.  Thus  quasl-linear  theory  is  not  applicable 
near  the  edge  of  the  injection  region  where  the  particles  enter,  but  may 
be  significant  where  they  exit. 

Physical  Interpretation.  Certain  physical  attributes  of  the  char- 
a  teristics  deserve  further  elaboration.  First,  they  can  be  expressed  In 
much  simpler  form  where  their  physical  interpretation  Is  more  transparent. 
Second,  the  effects  of  introducing  a  plasma  cloud  or  a  stream  of  hot 
particles  are  readily  evident  from  the  analytical  form  of  trig  expressions. 
Finally,  the  effect  of  off-axis  propagation  (k^  not  pa-allel  to  B)  is  ex¬ 
amined  qualitatively  to  assess  its  effect  on  stimulated  amplification  and 
the  energy  budget  of  the  magnetosphere. 

The  expression  for  the  phase  velocity  in  equation  (4)  can  be  simpli¬ 
fied  considerably  for  frequencies  just  below  the  electron  gyrofrequency 
(VLF)  and  the  ion  gyrofrequency  (ULF).  By  neglecting  terms  of  order  me/m^ 
the  index  of  refraction  has  the  two  forms 
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These  are  well-known  expressions  for  VLF  and  ULF  whistler-mode  propa¬ 
gation  which  display  the  characteristic  frequency  dispersion  that  Is 
observed  experimentally  In  the  magnetosphere. 

Regimes  of  wave  growth  (k^O)  and  decay  (ki  >0)  for  anisotropic 
distributions  can  be  formally  derived  from  equation  (5).  For  this  purpose 
It  Is  necessary  to  Introduce  the  concept  of  a  particle  ^tch  angle  a 
between  the  velocity  and  the  magnetic  field. 


a  *  tan  (-Vjl  /v|( ) 

Then  equation  5  can  be  expressed  In  the  form 

^’*]b 


(8) 


(9) 


For  an  Isotropic  distribution,  a  -  0  so  that  Its  summation  term  In 
equation  (9)  Is  always  positive  and  the  wave  Is  attenuated  by  that  term. 
The  special  class  of  anisotropies  *  s1nm  a  gives  a  particularly 
simple  result.  In  this  case,  a  *  m/2  and  the  normal  Interaction  terms 
yield  amplification  when  u/ 1  wc|<m/(m  +  2).  This  latter  anisotropic 
behavior  Is  characteristic  of  natural  distributions  In  the  magnetosphere 
at  VLF  and  ULF. 


The  physics  of  the  interaction  may  be  explained  in  terms  of  particle 
trapping  in  the  apparent  potential  well  of  the  wave  fields.  Although  all 
the  particles  in  the  plasma  are  subjected  to  the  electromagnetic  wave 
fields,  only  those  particles  with  vM  sufficiently  close  to  vc  actually 
participate  in  the  exchange  of  energy.  The  forces  tend  to  drive  v)( 
toward  vc  but  because  the  circularly  polarized  wave  has  a  helical  electro¬ 
magnetic  potential  well,  the  vxB.  forces  also  change  the  perpendicular 
particle  velocity  vj,  .  The  velocity  vector  of  each  particle  is  shifted 
slightly  to  bring  the  particle  Into  phase  with  the  potential  well  and 
this  process  leads  to  a  net  exchange  of  particle  kinetic  energy  and  electro 
magnetic  energy.  As  noted  above,  the  net  exchange  depends  rather  critical 
on  the  shape  of  the  phase-space  distribution,  particularly  the  pitch  angle 
distribution.  The  question  of  whether  more  particles  lose  or  gain  energy 
is  prescribed  by  the  derivatives  in  the  Integrands  of  equations  (5)  0r  (10) 
which  are  evaluated  at  Vn  *  vc.  In  other  words,  the  slope  of  the  distri¬ 
bution  function  is  a  measure  of  the  number  of  particles  that  are  giving  up 
energy  compared  to  those  that  are  gaining  it.as  they  shift  their  velocity 
to  follow  the  helical  potential  well.  The  distortion  of  F^q-j.  by  this 
process  is  ignored  in  the  linear  theory  for  reasons  which  were  discussed 
at  the  end  of  the  preceding  sub-section. 

Possible  ways  to  enhance  VLF  or  ULF  wave  amplification  In  the 
magnetosphere  are  now  readily  Identifiable  in  the  foregoing  theoretical 
equations.  The  expression  for  the  amplitude  exponent  k^  depends 
on  a  variety  of  parameters  and  distributions,  which  might  be  modified 
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locally  to  enhance  the  cyclotron-resonance  interaction.’  The  most  obvious 
means  of  stimulating  the  interaction  is  to  change  the  hot  plasma  dis¬ 
tribution  f^gy  Amplification  is  achieved  by  dumping  a  large  amount  of 
energetic  particles  near  the  equator  so  that  sinma  pitch  angle  dis¬ 
tributions  are  formed.  This  type  of  injection  is  feasible  by  nuclear 
bursts.  The  early  tests  which  were  detonated  in  the  upper  atmosphere 
created  pitch  angle  distributions  of  the  form  cosma  that  strongly  absorb 
VLF  and  IJLF  energy.  In  fact,  the  only  signals  detected  from  these  tests 
were  apparently  the  initial  electromagnetic  pulse  generated  by  the  local 
current  system  at  the  source.  A  second  method  for  causing  substantial 
amplification  is  particle  beam  injection.  Such  a  beam  would  be  additive 
to  the  hot  plasma  distribution  F^.  Its  contribution  would  only  occur 

when  the  parallel  component  of  the  beam  velocity  is  near  v  So  that  the 

c 

beam  can  participate  in  the  resonance  interaction.  In  the  following 

sections  electron  beams  are  investigated  as  a  source  of  VLF  amplification 

and  proton  beams  are  considered  at  ULF. 

Another  relatively  simple  procedure  for  stimulating  amplification 

is  to  enhance  the  background  plasma  density  in  a  local  region  where 

natural  amplification  is  occurring  (Brice,  1970;  Cornwall  and  Schulz, 

1971).  Such  an  enhancement  alters  the  index  of  refraction  (see  equations 

6  and  7)  for  the  modes,  reducing  the  phase  velocities  appreciably.  Since 

1/2 

the  resonance  speed  vc  in  equation  (3)  is  proportional  to  N"  ,  this 
increase  in  density  allows  lower  energy  particles  to  participate  in  the 
resonance  interaction.  If  the  hot  plasma  distribution  has  more  particles 
at  these  lower  energies,  the  resonance  interaction  is  sharply  stimulated, 
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and  considerable  amplification  is  antirin  *  A  * 

„  , .  ,  anticipated  from  catalytic  Injections  of 

Plasma.  Several  specific  cases  of  plasma  Injections  are  modeled  in 
the  following  sections  to  Illustrate  these  possibilities. 

These  foregoing  theoretic.,  expressions  are  strictly  valid  only 
for  propagation  vectors  jc  parallel  to  the  local  magnetic  field.  But  their 
applicability  to  propagation  conditions  in  the  magnetosphere  is  expected 
to  extend  far  beyond  this  narrow  subclass  of  propagation  conditions. 

Since  waves  that  a,,  not  parallel  to  the  local  field  are  subject  to  con¬ 
ventional  Landau  damping  (Kennel.  1966),  only  those  signals  with  wave 
t annals  near  the  local  magnetic  field  actually  survive  propagation 

tough  the  magnetosphere.  Consequently,  those  VLF  and  ULF  whistler 
s.Snals  that  are  detected  by  ground  stations  and  satellites  are  likely 

1°  ^  “nti,ned  1n  the  C'aSS  °f  Potion  vectors  i„cl„ded  1n  ^ 
foregoing  theoretical  description  Tn 

n„,  .  ,  description.  I„  this  sense  the  theory  of  circularly 

polarized  waves  is  quite  useful  nn  u 

SefUl-  0n  the  ottier  "and.  it  must  be  recalled 

that  the  wave  normals  change  their  direction 
+  .  airection  as  they  traverse  the  mag- 

e  osphere  (Kitamura  and  Jacobs.  1967)  so  that  much  of  the  wave  energy 
actually  is  ultimately  lost  through  the  Landau  process.  An  estimate  of  the 

fraction  of  energy  that  is  gained  from  particle  amplification  in  one  region 

or  the  magnetosphere  and  lost  bv  thi' 

.  y  attenuation  elsewhere  has  not  been 

Process  plays  an  Important 
a  redistribution  of  energy  throughout  the  magnetosphere  during 

—  disturbance.  It  must  also  play  a  significant  role  in 

e  quasi- steady-state  conditions  that  describe  the  dynamic  equilibrium 
e  ween  waves  and  particles  during  geomagnetical ly  quiet  periods.  This 
difficult  problem  needs  further  study. 
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Magnc tospherio  Models.  The  magnetosphere  consists  of  a  neutral 
(hydrogen)  plasma  imbedded  in  the  geomagnetic  field.  The  field  is  approxi¬ 
mately  a  dipole  subject  to  minor  internal  distortions  aita  a  major 
boundary  distortion  due  to  the  solar  wind.  The  cold  hydrogen  plasma  is 
distributed  throughout  this  geomagnetic  cavity  but  it  has  a  distinctive 
spatial  distribution  and  is  subject  to  considerable  fluctuations  with  geo¬ 
magnetic  activity.  In  addition  there  is  a  hot  plasma  component  of  con¬ 
siderably  lower  density  that  consists  of  energetic  (1-100  keV)  electrons  and 
protons  that  follow  trapped  orbits  along  the  geomagnetic  field  lines. 

Above  the  ionosphere  the  medium  is  assumed  to  be  slowly  varying  electro- 
magnetically  in  the  sense  that  the  wave  length;  A  satisfies  the  condition 

I  v  a  |=  2ttc {  vnl  /n2  to  «  1  (12) 

for  all  frequencies  of  interest.  On  this  basis,  the  propagation  character¬ 
istics  of  a  uniform  medium  that  are  given  in  the  preceding  section  are 
generally  applicable  locally  in  the  magnetosphere.  The  theory  of  VLF 
and  ULF  whistler  propagation  through  the  magnetosphere  is  based  on  the 
condition  (12).  At  VLF  the  whistler  wave  packet  propagates  along  the 
flux  tube  field  line  from  one  hemisphere  to  the  other  in  the  right-hand 
mode  of  circular  polarization  (e.g.,  Helliwell,  1965).  Similarly,  the  ULF 
waves  propagate  between  hemispheres  in  the  left-hand  mode  (Jacobs  and 
Watanabe,  1964).  Much  of  the  whistler  guidance  along  the  flux  tubes  is 
attributed  to  field-aligned  columns  of  enhanced  ionization  (Smith  et  al.. 
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1960;  Smith,  1961).  The  tine  delay  between  successive  hops  of  a  whistler 
signal  between  hemispheres  and  its  characteristic  frequency  dispersion  are 
described  very  accurately  by  path  integrals  of  the  electromagnetic  group 
velocity  for  a  locally  uniform  medium. 

Similarly,  the  amplitude  exponent  for  a  uniform  homogeneous 
medium  is  assumed  to  describe  thewave  growth  or  decay  along  the  propa¬ 
gation  path.  The  net  path  amplification  or  absorption  of  power  in  a 
wave  that  propagates  through  the  slowly  varying  magnetosphere  is 
given  by 

A  (decibels)  =  -10  log  exp  f  2  k.  (a>,s)  ds  (13) 

10L  •'path  1  J 


where  the  path  is  along  a  dipole  field  line  segment  from  the  interaction 
region  to  the  receivers  on  the  ground, and  k^  is  defined  by  the  integral  in 
equation  5.  Evaluation  of  the  double  integral  contained  iri  expression  (12) 
for  a  variety  of  injection  models  is  the  primary  objective  of  this  research. 

Without  significant  loss  of  quantitative  accuracy  the  geomagnetic 
field  may  be  approximated  by  a  pure  dipole  in  this  application.  Thus  the 

signal  path  is  defined  by  the  equation  for  a  field  line 

_  2 

R  =  L  cos  a  (14) 

where  R  is  the  radial  distance  in  earth  radii,  L  is  the  geomagnetic  equatorial 
intercept  in  earth  radii  and  A  is  geomagnetic  latitude.  The  electron-cyclotron 
frequency  along  this  path  has  the  form 


0)  =0) 
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(1+3  sin2 a)1/2  /cos6; 
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where  is  its  equatorial  value  (-5.5  x  10  L  radians/sec). 

At  L=4  this  corresponds  to  an  electron  gyrofrequency  of  13.7  kHz  and  a 
proto i  gyrofrequency  of  7.5  Hz. 

The  thermal  background  plasma  is  relatively  slowly  varying  along  a 
given  field  line  but  is  subject  to  a  large  discontinuity  across  the  plasma- 
pause  in  the  vicinity  of  L  =  3-5  (Carpenter  1966,  Carpenter et  al., 

1969).  Inside  the  plasmapause  the  thermal  plasma  is  in  hydrostatic  equil¬ 
ibrium  whereas  outside  the  density  varies  approximately  as  the  geomagnetic 
field.  For  the  analyses  presented  here,  the  plasma  frequency  is  modeled 
by  the  following  expressions, 

to®  =  cdpQ  exp  (3/R  -  3/L)  in  plasmasphere  (14) 

=  (dpQ  (L/R)3  in  plasmatrough  (15) 

where  uy  is  its  equatorial  value  (5.65  x  104  N.  radians/sec) 
pc  0 

and  the  hydrostatic  model  has  a  temperature  of  1250°  K.  For  modeling 

3 

purposes,  the  equatorial  density  is  assumed  to  be  NQ  =  232  particles/cm 
in  the  plasmasphere  which  gives  an  electron  plasma  frequency  of  137  kHz. 

3 

In  the  plasmatrough  NQ  =  1  particle/cm  which  corresponds  to  a  frequency 
of  9  kHz. 

In  order  to  separate  the  distribution  function  into  a  cold  com¬ 
ponent  and  a  hot  component,  the  frequency  band  of  interest  must  be 
somewhat  below  the  local  gyrofrequency  for  the  resonant  particles.  In 
practice  this  requires  w^0.9  |  |  for  resonance  energies  above  10  eV. 

For  the  applications  presented  here,  this  arbitrary  energy  limit  establishes 


the  division  between  FCqLq  and  F^gy.  At  these  energies  the  mean 
free  path  of  the  hot  plasma  particles  are  so  long  that  they  follow 
conventional  first-order  trapped  orbits  in  the  magnetosphere.  The 
cold  plasma,  however,  behaves  l'ke  a  Maxwellian  gas  similar  to  its 
diffusive  source,  the  ionosphere. 

According  to  Liouville's  theorem,  the  hot  plasma  distribution, 

F-'HgT,  is  prescribed  everywhere  in  terms  of  its  equatorial  distribution 
Since  trapped-particle  data  is  generally  in  the  form  of  separate  pitch- 
angle  and  energy  distributions,  the  models  for  are  assumed  to  have 
’he  form 

fhot  *  g(v)  H(a)  (16) 

Along  the  flux  tube  G  is  obviously  invariant  because  particle  energy 
is  conserved  in  the  geomagnetic  fie^d,  but  H  is  not  constant.  Using 
the  Liouville  theorem  and  the  conservation  of  magnetic  moment,  it  is 
easily  proved  that 

H(a)  =  U0  (sin’1  [(Bo/B)1/2  sin  a])  (17) 

where  UQ  is  the  equatorial  distribution. 

F^g-p  is  normalized  to  the  experimental  data  in  the  following  manner. 

The  pitch-angle  distribution  is  arbitrarily  normalized  to  unity, 

W  2 
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and  the  speed  distribution  is  based  on  experimental  differential  density 
spectra,  d'l/dE,  which  are  related  as  follows: 

d'l/'lQ  dE  =  G  v2  dv/dE  (19) 

The  experimental  data  for  energy  and  pitch  angle  distributions  of  interest 
in  this  application  are  introduced  in  the  following  sections  of  the 
document. 

A  computer  program  for  evaluation  of  the  cyclotron-resonance  inter¬ 
action  has  been  prepared  based  on  the  foregoing  theoretical  expressions. 

A  complete  listing  of  the  program  and  Sts  sub-routines  is  included  in 
an  Appendix.  The  program  is  designed  to  handle  both  electron  interactions 
with  VLF  and  proton  interactions  with  ULF  separately.  The  output  of  the 
program  includes  listings  of  the  amplification  exponent  k.  as  a  function  of 
geomagnetic  latitude  along  the  propagation  path  and  the  net  path  amplifi¬ 
cation  A  for  a  range  of  frequencies.  In  general,  these  two  quantities 
completely  specif}  the  quantitative  characteristics  of  the  interaction. 

A  variety  of  quantitative  models  for  stimulation  of  wave  amplification 
are  treated  in  the  following  sections. 


VLF  AMPLIFICATION 


The  artificial  stimulation  of  VLF  amplification  in  the  magnetosphere 
is  relatively  easy  to  achieve  because  of  the  variety  of  injection  materials 
and  conditions  that  are  available.  However,  the  energetic  electron 
distribution  that  is  responsible  for  amplification  changes  its  shape 
(Frank,  1967,  1968)  depending  upon  the  state  of  disturbance  in  the 
magnetosphere.  Thus,  it  is  important  to  choose  the  proper  time  for  injecti 
of  piasma  clouds  that  can  induce  amplification  via  the  cyclotron-resonance 
interaction.  Several  cases  of  interest  are  modelled  quantitatively  in 
the  following  subsections. 

Before  discussing  specific  cases  it  is  necessary  to  introduce  the 
models  for  the  energetic  electron  pitch  angle  and  energy  distributions. 

The  pitch  angle  distributions,  H  ,  are  assumed  to  be  of  tie  form  *  sinma 
which  are  typical  of  experimental  observations  in  the  magnetosphere 
(O'Brien,  1963).  In  practice  the  model  distributions  are  distorted  at. 
small  pitch  angles  in  order  to  introduce  an  atmospheric  loss-cone  cutoff 
at  an  equatorial  pitch  angle  of  10°.  Some  cases  of  interest  are  illustrated 
in  Figure  1  where  the  isotropic  distribution  corresponds  to  m  =  0  and 
the  two  anisotropic  cases,  m  =  1  and  m  =  2,  are  typical  of  the  energetic 
electrons  in  the  magnetosphere. 

As  indicated  in  the  preceding  theoretical  discussion,  the  energy 
distribution  may  he  subdivided  into  a  cold  and  a  hot  component.  The 
cold  background  plasma  is  relatively  more  dense  and  is  assumed  to  have 
a  Maxwellian  energy  distribution  with  a  mean  energy  on  the  order  of  0.1  eV. 


DISTRIBUTION  H 
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PITCH  ANGLE  a  (deg) 


Fig.  1  Pitch  Angle  distribution  models  H  at  the  geomagnetic  equator. 
These  anisotropic  models  including  a  loss  cone  cutoff  are  used 
for  both  electrons  and  protons  in  the  magnetosphere. 
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a?e  DprfnrmpH  ?"  he  energy  scale-  Some  of  the  calculation 
a  e  performed  for  energy  distributions  of  the  form  E-n  that 

are  typical  of  many  observations.  The  substorm  distribute 
A'  was  measured  on  the  ATS-5  satellite  at  geosynchronous 
orbit  (DeForest  and  Mcllwain,  1971).  yeesyncnronous 
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It  controls  the  propagation  characteristics  of  VLF  waves  but  transfers 
negligible  energy  to  the  waves.  The  hot  electrons  In  the  energy  range 
from  0.1  to  100  keV  are  responsible  for  the  resonance  energy  exchange 
at  VLF  In  the  magnetosphere.  Typical  electron  distributions  which  ere 
used  here  are  illustrated  In  Figure  2.  The  substorm  distribution  "A" 
is  based  on  satellite  observations  (DeForest  and  Mcllwaln,  1971)  located 
at  L  =  G.G  in  a  geosynchronous  orbit.  The  distributions  E“n  where  n  «  1  or  2 
are  fairly  good  approximations  to  conditions  throughout  the  region  inside 
the  geosynchronous  orbit.  Including  the  transition  across  the  plasmapause. 

Natural  Density  Enhancements.  Ourlng  disturbed  conditions  in  the 
magnetosphere  the  location  of  the  plasmapause  can  shift  appreciably  In  Its 
L-shell  location  causing  local  density  variation*  of  1-2  orders  of  magnitude. 
Such  changes  can  cause  major  modifications  In  the  energy  exchange  between 
VLF  waves  and  energetic  electrons  which  results  In  a  major  revision  of  the 
dynamic  equilibrium  state  of  the  magnetosphere  (Brice  and  Lucas,  1971). 

The  specific  case  chosen  for  analysis  Is  assumed  to  occur  Inside  the 
plasmasphere  where  the  plasma  density  is  assumed  to  Increase  by  a  relatively 
small  amount.  However,  the  increase  In  amplification  may  account  for  the 
VLF  emission  and  whistler  activity  that  Is  observed  In  the  plasmasphere. 

The  cold  plasma  is  modeled  by  the  hydrostatic  equilibrium  described  by 
equation  (14).  A  propagation  path  along  L  =  4  Is  assumed  for  Illustration 

_3 

with  an  ambient  plasma  density  of  «  232  cm  and  an  Inhancement  to 

N  =  733  cm'3. 

0 


The  hot  electron  distribution  is  assumed  to  have  the  form  E"1,5  sin  a 
corresponding  to  the  models  given  in  Figures  1  and  2.  The  change  in  the 
amplification  exponent  due  the  density  enhancement  Is  illustrated  in 
Figuro  3  for  two  frequencies  of  interest.  Although  the  change  in  the 
amplification  exponent  appears  relatively  small,  the  net  path  amplification 
is  Increased  by  about  5  db  (decibels)  as  shown  In  Figure  4. 

These  figures  illustrate  some  important  features  of  the  amplification 
exponent  and  the  net  path  arpl Ification  at  VLF.  First,  the  exponent, 
varies  rapidly  with  geomagnetic  latitude  and  depends  strongly  on  the 
°cal  plasma  and  cyclotron  frequency  parameters  as  well  as  the  propagation 
frequency.  Its  value  at  the  geomagnetic  equator  is  not  a  good  indicator 
of  its  overall  Average  throughout  the  region  of  interaction.  The 
amplification  actually  peaks  sharply  some  10c  to  20°  from  the  geomagnetic 
equator.  At  the  higher  frequencies  the  signal  is  sharply  attenuated  near 
the  geomagnetic  equator  due  to  the  interaction  with  lower  energy  particles 
which  can  enter  the  resonance  band.  The  Interaction  decreases  sharply 
at  some  30°  to  40°  from  the  geomagnetic  equator  because  the  Doppler  shift 
becomes  too  large  for  even  very  energetic  particles  to  contribute 
appreciably  to  the  amplification. 

The  shape  of  the  net  path  amplification  depicted  in  Figure  4  is 
typical  of  the  natural  growth  that  can  be  expected  at  VLF.  The  location 
of  the  frequency  at  which  maximum  amplification  occurs  and  the  location 
of  the  upper  cutoff  where  amplification  ceases  depends  primarily  on  the 
s^ape  of  the  pitch  angle  distribution,  i.e.,  exponent  n.  The  relationship 


AMPLIFICATION  EXPONENT, 


26, 


Fig.  3  Increased  amplification  exponent  for  VIF  waves  due  to  a  natural 
enhancement  In  the  cold  plasma  density.  The  propagation  at 
L  *  4  Is  assumed  to  be  Inside  the  plasmasphere  where  the  cold 
plasma  Is  In  hydrostatic  equilibrium.  For  this  case,  the  hot 
plasma  Is  assumed  to  have  an  energy  and  pitch  angle  distribution 
of  the  form  sin  «  (see  Figures  1  and  2). 
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Fig.  4  Increased  net  path  amplification  of  VLF  waves  due  to  a  natural 
enhancement  in  the  cold  plasma  density.  The  propagation  path 
at  L  =  4  is  assumed  to  be  inside  the  plasmasphere  where  the 
cold  plasma  is  in  hydrostatic  equilibrium.  For  this  case, 
the  hot  plasma  is  assumed  to  have  an  energy  and  pitch  angle 
distribution  of  the  form  E-1 . 5  sin  a  (see  Figures  1  and  2). 
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between  the  upper  cutoff  frequency  and  this  anisotropy  index  m  was 
discussed  following  equation  (9)  above.  The  cutoff  location  is  shifted 
somewhat  in  the  inhomogeneous  magnetosphere  where  the  effect  is  smeared 
out  by  the  variable  parameters  of  the  background  medium.  The  magnitude 
of  the  amplification  in  general  depends  on  the  density  of  the  energetic 
particles  available  to  exchange  energy  with  the  wave  at  the  resonance 
velocity.  The  shape  of  the  curve  for  amplification  depends  on  the 
shape  of  the  energy  distribution,  that  is,  the  value  of  n.  As  n  increases 
the  bandwidth  of  amplification  tends  to  increase  with  it.  Of  course 
all  of  these  characteristic  variations  depend  on  both  m  and  n  in  a 
complicated  way  which  is  difficult  to  separate.  More  complete  discussions 
of  these  dependencies  on  the  shape  of  the  distribution  function  can  be 
found  in  the  earlier  literature  (e.g. ,  Liemohn,  1967). 

Barium  Jet  in  the  Plasmatrough.  The  technology  for  rocket  injection 
of  barium  has  progressed  to  the  point  where  jets  can  je  propelled  along 
the  geomagnetic  field  lines  by  shaped  charges  located  in  the  nose  cone  of 
the  rocket  payload.  Weight  limitations  on  the  launch  vehicle  limits  the 
amount  of  barium  to  approximately  one  kgm  which  corresponds  to  5  x  1024 
barium  atoms.  Energy  from  the  shaped  charge  as  well  as  the  sunlight 
rapidly  ionizep  this  barium  so  that  the  ions  and  electrons  are  trapped 
along  the  geomagnetic  field.  Apparently  plasma  instabilities  occur  in  the 

Shaped  charge  injections  of  barium  were  carried  out  in  a  series  of 
experiments  from  Alaska  and  Hawaii  in  March  and  October,  1972,  by 
Dr.  Milton  Peek  of  the  Los  Alamos  Scientific  Laboratories.  The  jet 
of  neutral  barium  was  quickly  ionized  and  followed  the  geomagnetic 
field  line  to  the  opposite  hemisphere  from  its  injection  region.  The 
purpose  of  these  experiments  was  to  map  the  geomagnetic  field  by 
following  the  optical  emission  from  the  cloud  of  barium  ions. 
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Charnel  cloud  and  cause  stations  that  reduce  the  effective  density  of 

the  hear,,  nevertheless  tests  have  shown  that  an  appreciable  fraction  of 

the  injected  barium  does  progress  up  the  field  line  toward  the  other 
hemisphere. 

The  leading  edge  of  the  barium  jet  travels  at  about  14  km/sec  and 

the  trailing  edge  at  9  km/sec.  Thus,  at  the  geomagnetic  equator  along 

the  field  line  at  L  -  4,  the  equatorial  enhancement  is  approximately 

15.000  km  long  and  about  15  km  in  diameter.  This  corresponds  to  a  volume 

of  approximately  3  x  102'  cm3.  From  the  intensity  of  the  optical  emission 

the  equatorial  density  is  approximately  2-3  barium  lons/cm3.  Approximate!; 

40-50  minutes  is  required  for  the  jet  to  reach  the  geomagnetic  equator 
along  L  =  4. 

On  the  basis  of  the  foregoing  parameters  and  experimental  evidence, 
a  model  of  the  barium  injection  has  been  constructed  to  evaluate  its  effect 
on  the  cyclotron  resonance  interact  on.  Some  15  minutes  after  the  release 
along  L  •  4  the  primary  jet  is  bounded  between  30”  and  38”  from  the 
geomagnetic  equator  and  has  a  density  of  approximately  30  ions/cm3.  At 
30  minutes,  the  jet  is  bounded  by  8”  and  25”  from  the  geomagnetic  equator, 
and  the  density  has  reduced  to  5  lons/cm3.  The  reduction  in  density  is 
due  to  the  increased  diameter  of  the  jet  as  well  as  its  elongation.  In 
order  for  these  enhancements  to  be  significant  relative  to  the  natural 
Plasma  background,  the  injection  must  occur  during  a  mild  geomagnetic 
substora  when  the  plasmapause  is  inside  L  *  4.  The  ambient  natural  density 

in  the  p’asmatrough  is  assumed  to  be  1  electron/cm3  and  vary  radially  as 

as  -'-3.  These  properties  of  the  barium  shaped-charge  injection  model 
are  illustrated  in  Figure  ‘j. 
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BARIUM  SHAPED  - 
CHARGE  INJECTION 
L  =  4.0  TRAJECTORY 


Fig.  5  Rocket  injection  of  a  barium  shaped-charge  along  the  L 
field  line.  The  barium  jet  is  assumed  to  consist  of  1 
fully  ionized  gas.  The  trajectory  iv  assumed  to  be  in 
plasmatrough,  where  NQ  =  1  electron/!  m- 3  and  N  varies 


=  4 
kgm  of 
the 

as  R"3. 


An  important  aspect  of  the  problem  that  is  not  included  explicitly 
in  the  amplification  calculation  is  the  diameter  of  the  barium  ion  jet 
relative  to  the  size  of  the  VLF  wave  front.  As  noted  in  the  theoretical 
section,  it  is  important  that  any  enhancement  be  comparable  to  or  larger 
than  a  wavelength  in  all  of  its  dimensions.  Using  the  experimental 
observation  of  15  km  as  the  diameter  of  the  barium  ion  jet  at  the  equator, 
the  diameter  at  15  minutes  is  approximately  7.5  km  and  at  30  minutes, 

11  km.  For  a  propagation  frequency  of  6  kHz,  which  is  approximately 
half  of  the  quatorial  electron  gyrof requency  at  L  =  4,  the  free  space 
wavelength  is  about  50  km.  However,  the  index  of  refraction  within  the 
barium  cloud  enhancement  is  approximately  10  in  both  cases  so  that  the 
effective  wavelength  is  somewhat  less  than  the  actual  diameter  of  the 
cloud.  Thus,  to  a  first  approximation,  the  cyclotron  resonance  interaction 
should  not  be  limited  by  the  width  of  the  jet. 

The  amplification  exponent  for  the  interaction  of  VLF  waves  within 
the  barium  jet  at  15  and  30  minutes  is  illustrated  in  Figure  6.  Evidently 
the  interaction  is  strongly  enhanced  as  the  jet  approaches  the  geomagnetic 
equator,  which  is  not  unexpected  due  to  the  smaller  Doppler  shift  that  is 
required.  Furthermore,  the  exponent  is  strongly  enhanced  only  in  the  local 
region  of  the  barium  jet.  Note  the  opposite  behavior  of  the  exponent, 
kj,  for  the  propagation  frequencies  0.3  and  0.5  f®Q  at  15  and  30  minutes, 
respectively.  This  is  readily  explained  by  channes  in  the  relative 
magnitudes  of  the  growth  and  decay  terms  in  expression  (2).  The  hot  plasma 
which  provides  the  energy  for  the  amplification  is  arbitrarily  assumed  to 
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Fig.  6  Amplification  exponent  for  VLF  waves  propagating  through  the 
barium  shaped-charge  injection  along  L  =  4.  The  hot  electrons 
dr©  assumed  to  have  the  distribution  sin^a  (see  Fiqures  1 
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Net  path  amplification  for  VLF  waves  propagating  through  the 
barium  shaped-charge  injection  along  L  =  4.  The  hot  electrons 
are  assumed  to  have  the  distribution  E_2  sin^-a  (see  Figures  1 
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have  the  form  E  sin  a  which  is  a  reasonable  approximation  for  the 
electron  energy  distribution  in  the  plasmatrough  during  a  geomagnetic  storm. 

The  net  power  that  can  be  transferred  from  the  energetic  electron 
distribution  to  VLF  waves  by  the  barium  shaped-charge  enhancement  is 
illustrated  in  Fig  re  7.  Evidently  the  amplification  some  15  minutes 
after  the  injection  is  still  negligible  below  ,  but  after  30  minutes 
the  amplification  has  increased  20-30  db  and  the  bandwidth  has  been 
expanded  somewhat.  At  L  =  4,  it  has  been  noted  that  feQ  =  13.6  kHz,  so 
that  the  30  minute  amplification  band  is  between  4  and  8  kHz. 

Geosynchronous  Lithium  Clouds.  One  of  the  more  promising  methods 
for  inducing  strong  amplification  of  waves  in  the  magnetosphere  is  by 
injection  of  a  cloud  of  lithium  at  the  geosynchronous  orbit.  Lithium  has 
a  relatively  slow  rate  of  ionization  in  solar  uu  'violet  so  that  the 
cloud  of  neutral  atoms  can  diffuse  radially  for  approximately  an  hour 
before  most  of  it  is  ionized  and  trapped  by  the  geomagnetic  field.  Although 
this  slow  rate  of  lithium-ion  generation  does  not  permit  its  application  as 
a  communication  method,  the  diameter  of  the  cloud  makes  it  a  strong  contender 
for  a  quantitative  test  of  the  theory.  Furthermore,  the  large  area  that 
would  be  illuminated  with  VLF  noise  at  the  ground  makes  it  attractive  from 
an  observational  standpoint. 

The  distribution  of  the  geosynchronous  lithium  injection  has  been 
modeled  by  the  solution  of  the  diffusion  equation  for  a  point  source  of 
particles.  The  density,  N,  is  given  as  a  function  of  radial  distance,  r, 
(from  the  source  point)  and  time,  t  (measured  from  an  arbitrary  instant 
of  the  release),  by  the  equation, 
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N  =  N1  exp  (-r2/4D2t2)/(4*D2t2)3/2  (20) 

where  is  the  number  of  lithium  ions  injected  and  D  is  the  diffusion 
coefficient  for  the  cloud  of  particles.  In  such  a  solution  no  attempt 
is  made  to  include  the  dynamic  process  for  ionization  of  the  lithium 
atoms,  instead,  it  is  assumed  that  N  describes  the  current  number  of 
ions  in  the  cloud  at  time  t. 

For  this  application  the  value  of  N-|  is  assumed  to  be  8.5  x  102^ 
molecules  corresponding  to  1  kgm  of  fully  ionized  lithium.  The  diffusion 
coefficient,  D,  is  assumed  to  be  1  km2/sec,  which  gives  cloud  diffusion 
rates  in  good  agreement  with  ooservations.  The  induced  amplification 
has  been  evaluated  for  distributions  corresponding  to  10  minutes  after 
release  and  50  minutes  after  release  of  the  lithium  charge.  The  electron 
density  enhancements  corresponding  to  these  release  times  are  illustrated 
in  Fig.  8  as  a  function  of  geomagnetic  latitude.  The  natural  background 
density  is  arbitrarily  assumed  to  be  1  electron/cm3  which  is  a  reasonable 
limit  at  the  geosynchronous  orbit  (L  =  6.6).  Since  it  requires  approxi¬ 
mately  one  hour  for  sunlight  to  fully  ionize  the  lithium  atoms  in  the 
cloud,  it  requires  at  least  five  times  as  much  lithium  to  achieve  one 
kilogram  of  ionization  in  10  minutes  as  it  does  in  50  minutes.  If  the 
lithium  release  from  the  canister  is  only  10  percent  efficient,  as  some 
researchers  have  suggested  (Cornwall  private  communication),  then  the 
amplification  enhancement  at  10  minutes  requires  a  lithium  payload  of 
50  kgm  plus  a  detonating  charge  whereas  the  amplification  at  50  minutes 
requires  only  10  kgm  of  lithium  plus  charge. 
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Fig.  8  Injection  of  lithium  from  a  point  source  (cannister)  in  geo¬ 
synchronous  orbit  (L  =  6.6).  One  kgm  of  total  ionization  is 
assumed  to  be  diffusing  spherically  from  the  source  with  a 
diffusion  coefficient  D  ~1  Km2/s.  The  ambient  background 
plasma  is  assumed  to  be  in  the  plasmatrough  where  NQ  =  1  cm  . 


The  effect  of  this  geosynchronous  lithium  injection  on  VLF  wave 

amplification  is  illustrated  for  two  hot  plasma  distributions  in 

Figures  9-12.  For  one  case,  the  hot  electron  distribution  is  assumed 
-1.5. 

to  have  the  form  E  sin  a  (see  Fig.  2)  which  is  taken  to  represent 
conditions  at  the  geosynchronous  altitude  during  geomagnetical ly 
quiescent  conditions.  The  second  case  is  modelled  after  experimental 
observations  of  a  typical  substorm  (DeForest  and  Mcllwain,  1971)  and 
is  denoted  by  the  functional  form  (dN/dE)„A„  sin  a  (see  Fig.  2).  The 
VLF  amplification  results  for  these  two  distributions  are  discussed 
together  in  order  to  emphasize  the  physical  aspects  of  the  interaction. 

The  first  and  most  obvious  effect  of  the  lithium  injection  is 
the  concentrated  enhancement  of  the  amplification  exponent,  k^  in 
the  vicinity  of  the  geomagnetic  equator  (see  Figures  9  and  11),  which 
is  entirely  attributable  to  the  geometry  of  the  enhancement.  The 
second  point  of  interest  is  the  similar  magnitudes  of  ki  for  the  two  hot 
electron  distributions,  but  note,  however,  that  the  frequencies  are 
significantly  lower  for  the  distribution  E"1,5.  Outside  the  frequency 
bands  shown,  the  enhancement  in  the  amplification  is  relatively 
negligible.  Since  all  other  parameters  are  held  constant,  the  entire 
effect  is  attributable  to  the  shape  of  the  energy  distribution,  and 
points  out  the  need  for  detailed  observations  of  the  phase-space 
distribution  for  an  adequate  quantitative  test  of  the  cyclotron 
resonance  amplification. 

The  net  path  amplification  from  these  two  hot  electron  distributions 

are  illustrated  in  Figures  10  and  12.  The  distribution,  E'1'5,  provides 

strong  amplification  at  frequencies  below  0.35  f  e. 

co 


After  10  minutes 


AMPLIFICATION  EXPONENT, -ckj/w 


38. 


9 


LATITUDE  (deg) 


Fig.  9  Amplification  exponent  for  VLF  waves  propagating  through  the 
geosyncrhonous  lithium  injection.  The  exponent  is  sharply 
peaked  in  the  equatorial  enhancement  which  has  a  mean  diameter 
of  about  6400  km  at  50  minutes  after  injection.  The  hot 
electrons  are  assumed  to  have  the  distribution  sin  a. 
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Fig.  10  Net  path  amplification  for  VLF  waves  propagating  through  a 

geosynchronous  lithium  injection.  The  enhancement  is  entirely 
attributable  to  the  lithium  bubble.  The  hot  electrons  are 
assumed  to  have  the  distribution  E"1-5  sin  a. 
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Fig.  11  Amplification  exponent  for  VLF  waves  propagating  through  a 

geosynchronous  lithium  bubble.  The  hot  electron  distribution 
Is  based  on  Model  'A'  for  a  typical  geomagnetic  substorm 
(see  Fig.  2). 
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Fig.  12  Net  path  amplification  of  VLF  waves  propagating  through  a 

geosynchronous  lithium  bubble.  The  hot  electron  distribution 
is  based  on  Model  'A'  for  a  typical  geomagnetic  substorm. 
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the  injection  does  provide  about  30  db  of  enhanced  amplification  at 
0.1  -  0.2  f£o,  but  after  50  minutes  no  appreciable  enhancement  is 
discernible.  This  is  to  be  contrasted  with  the  amplification  for  the 
substorm  "A"  where  the  ambient  amplification  is  shifted  to  the  frequency 
band  0.4  -  0.55  f®Q.  This  frequency  shift  and  the  reduction  in  ambient 
amplification  by  a  factor  of  20  is  attributable  to  the  depletion  of 
low-energy  particles  in  F|jflT  for  distribution  "A".  For  this  substorm, 
the  amplification  after  10  minutes  is  extremely  high,  on  the  order  of 
200  db.  Such  an  amplification  level  could  only  be  realized  if  the  VLF 
source  is  low-power  incoherent  particle  emissions;  normal  VLF  emissions 
subject  to  this  much  amplification  would  become  nonlinear  and  distort 
the  hot  plasma  phase-space  distribution  suppressing  the  net  amplification. 
From  practical  considerations,  the  most  important  result  is  the 
amplification  obtained  50  minutes  after  injection  during  a  substorm. 

For  this  case,  the  amplification  has  been  enhanced  some  40  db  at  a 
frequency  around  1.5  kHz.  Such  an  effect  should  be  readily  discernible 
from  ground-based  observations  of  VLF  in  the  vicinity  of  the  geosynchronous 
conjugate  points. 
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Energetto  Election  Beams.  The  prospect  of  smpj  ifying~VLF- waves 
by  direct  injection  of  hot  electron  beams  is  attractive  because  the 
net  amplification  is  relatively  insensitive  to  the  state  of  geomagnetic 
activity.  The  state  of  the  natural  background  plasma  is  relatively 
unimportant  since  the  energy  for  wave  growth  is  to  be  derived  entirely 
from  the  injected  beam  of  particles.  Thus,  the  process  is  no  longer 
catalytic  in  the  sense  of  utilizing  the  hot  electrons  that  are  trapped 
in  the  geomagnetic  field,  but  rather  is  a  direct  stimulation  of  the 
cyclotron  resonance  interaction.  The  entire  process  is  quite  analogous  to 
transmitting  waves  from  an  antenna  on  board  a  spacecraft.  The  amount 

of  amplification  is  obviously  proportional  to  the  intensity  of  the 
beam. 

In  the  analysis  presented  here,  the  beam  is  assumed  to  be  launched 
by  a  spacecraft  along  the  geomagnetic  field  line  at  L  =  4  which  is 
assumed  to  be  inside  the  plasmasphere  where  N0  =  232  electrons/cnf*. 

The  natural  hot  plasma  background  is  assumed  to  have  a  distribution  of 
the  form  E  sin  a.  These  quantities  merely  serve  to  define  the  local 

propagation  conditions  and  the  natural  amplification  along  the  propagation 
path. 

The  electron  gun  is  assumed  to  have  a  power  level  of  0.5  kilowatts, 
for  definiteness,  the  gun  is  assumed  to  generate  100  millisecond  pulses 
of  electron  current  with  a  beam  energy  of  10  keV.  At  this  energy, 
the  electrons  have  a  speed  of  6  x  10^  km/sec  and  an  equatorial  gyro 
radius  (at  L  =  4)  of  approximately  1  km,  depending  of  course  on  the 
local  pitch  angle.  Initially,  such  a  beam  would  appear  to  be  a  helix 
that  is  wrapped  around  the  injection  field  line  at  the  appropriate  pitch 
angle.  However,  electrostatic  forces  and  plasma  instabilities  are 
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expected  to  smear  out  this  geometrical  form  within  100  km  or  so  from  the 
injection  point  (Pellat,  private  communication),  and  the  beam  is  more 
appropriately  described  as  a  flux-tube-column  of  electrons.  Assuming  an 
order  of  magnitude  inflation  in  the  cross  section  of  the  beam,  the  volume 
of  the  column  is  2  x  105  km3.  Since  0.8  x  1018  electrons  are  contained  in 
the  pulse,  the  density  of  the  beam  is  0.004  electrons/cm3.  If  the  beam  has 
an  effective  spread  of  2  keV,  this  density  corresponds  to  a  differential 
energy  spectrum  of  0.002  electrons/cm3  keV. 

The  foregoing  considerations  were  employed  to  derive  the  electron 
beam  model  shown  in  Figure  13.  The  pitch  angle  distribution  Hq  is  normalized 
to  unity  and  has  a  pitch  angle  spread  of  10°  corresponding  a  limited  amount 
of  scattering.  The  energy  distribution,  dN/dE,  has  been  normalized 
to  the  value  derived  above  with  an  appropriate  energy  spread  which  is 
attributable  to  scattering.  In  order  to  assess  the  effect  of  injection 
location  and  injection  pitch  angle,  three  values  of  the  equatorial  pitch 
angle  are  arbitrarily  assumed  for  the  beam  model.  For  ao  =  85°  the  beam 
mirrors  at  a  latitude  of  +  3°,  for  55°  the  latitude  is  +  18°,  and  for 
25°  the  latitude  is  +37°.  Of  course  the  beams  must  be  injected  above 
these  corresponding  mirror  point  latitudes  in  order  to  achieve  the  model 
parameters. 

The  size  of  the  beam  cross  section  relative  to  the  wavelengths  of 
interest  determines  the  effectiveness  of  the  interaction.  Near  the 
geomagnetic  equator  at  L  =  4,  the  electron  gyro  frequency  is  13.7  kHz 
and  for  a  plasma  density  NQ  =  232  electrons/cm3.  The  corresponding 
plasma  frequency  is  137  kHz.  For  a  propagation  frequency  at  half  the 
gyro  frequency,  7  kHz,  the  index  of  refraction  is  approximately  20, 
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3  Electron  beam  model  distributions  for  pitch  angle,  H0,  and 
differential  energy,  dN/dE.  The  beam  distributions 
are  sharply  peaked  in  pitch  angle  and  energy  assuming  a 
particle  gun  deploys  the  beam  from  a  satellite  or  very  hiqh 
altitude  rocket. 
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and  the  wavelength  is  2.5  km.  This  is  to  be  compared  with  a  beam  injection 
diameter  of  about  2  km.  If  the  beam  cross  section  expands  by  a  factor 
of  10,  as  assumed,  the  cyclotron-resonance  interaction  between  the  beam 
and  ambient  VLF  signals  could  be  quite  effective.  However,  it  must  be 
remembered  that  the  interaction  flux  tube  maps  onto  a  very  small  region 
in  the  ionosphere  so  that  the  actual  VLF  energy  flux  may  still  be 
extremely  small. 

As  noted  in  the  introduction,  electron  beams  with  the  foregoing 
properties  have  been  generated  experimentally  (Hendrickson  et  al ,  1971; 
Cartwright  and  Kellogg,  1971;  and  Hess  et  al,  1971).  These  beam  experi¬ 
ments  behaved  qualitatively  as  expected,  but  many  questions  remain  about 
the  the  quantitative  aspects  of  their  behavior.  One  important  conclusion 
is  that  the  well-known  two-st.ream  plasma  instability  does  not  effectively 
break  up  the  beam  as  might  be  expected.  Evidently  there  is  enough 
spread  in  the  thermal  energy  of  the  beam  to  quench  the  electrostatic 
waves  that  are  predicted  on  the  basis  of  delta-function  beam  distributions. 
These  experiments  and  their  analyses  are  continuing. 

The  cyclotron-resonance  interaction  between  electron  beams  and 
VLF  waves  has  some  unusual  properties  due  to  the  shape  of  the  phase-space 
distribution.  Adjacent  to  the  central  pitch  angle  and  energy  values  of 
the  beam,  the  distribution  has  exceedingly  large  derivatives  which  play 
a  dominant  role  in  the  integrand  of  the  amplification  exponent,  k^ , 
displayed  in  equation  (2).  Thus,  in  the  vicinity  of  the  resonance,  there 
is  strong  absorption  on  one  side  and  strong  amplification  on  the  other 
due  to  the  change  in  sign  of  the  derivative.  This  behavior  is  clearly 
illustrated  in  Fig.  14  for  the  case  aQ  =  55°  and  EQ  *  10  keV.  For  a 
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Fig.  14  Amplification  exponent  for  VLF  waves  propagating  tnrough  an 
electron  beam  along  the  field  line  at  L  =  4  inside  the 
plasmapause.  The  electron  beam  parameters  are  an  equatorial 
pitch  angle  of  55°  and  a  particle  beam  energy  of  10  keV.  See 
the  text  for  an  explanation  of  the  shape  of  these  curves. 
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given  frequency  the  interaction  is  confined  to  a  narrow  band  of  geo¬ 
magnetic  latitude  where  the  parallel  component  of  the  beam  velocity 
matches  the  resonance  velocity  given  by  equation  (3).  The  fact  that 
the  absorption  occurs  on  the  equatorward  side  and  the  amplification 
on  the  earthward  side  of  the  interaction  region  is  simply  due  to 
the  values  of  the  derivatives  in  ki  at  the  local  resonance  velocity. 

The  location  of  the  resonance  region  is  further  from  the  equator  at 
the  higher  frequencies;  this  is  wholly  attributable  to  a  reduction 
in  the  parallel  velocity  component  of  the  beam  as  it  progresses 
toward  its  mirror  point  and  necessitates  frequencies  closer  to  the 
local  gyrofrequencies  to  achieve  resonance. 

For  the  beam  energy  chosen  in  this  analysis,  the  pitch  angle 
distributions  centered  at  85°  and  25°  do  not  contribute  significant 
effects.  At  85°  the  entire  interaction  is  extremely  close  to  the 
geomagnetic  equator  and  the  parallel  velocity  component  is  relatively 
small  demanding  a  frequency  very  close  to  the  local  cyclotron  frequency 
where  the  amplification  is  also  small.  At  25°  the  beam  has  more 
parallel  than  perpendicular  energy  in  the  region  of  interaction  so  that 
the  wave  energy  is  absorbed  by  the  particles,  but  again  the  interaction 
is  relatively  small . 

The  net  path  amplification  for  the  three  cases  is  illustrated  in 
Fig.  15.  Evidently,  the  interaction  at  25°  and  85°  is  essentially 
neutralized,  that  is,  the  absorption  of  VLF  energy  just  balances  the 
amplification.  However,  at  55°  the  absorption  is  clearly  significant 
above  0.3  fCQ.  For  this  pitch  angle  the  interaction  is  strongly  dominated 
by  the  equatorward  absorption. 


NET  PATH  AMPLIFICATION  (db) 


49. 


Fig.  15  Net  path  amplification  for  VLF  waves  propagating  through  an 
electron  beam  at  L  =  4  inside  the  plasmapause.  Equatorial 
pitch  angles  for  the  beam  are  25°,  55°,  and  85°  and  the  beam 
energy  is  10  keV.  See  the  text  for  an  explanation  of  the 
shape  of  these  curves. 
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These  foregoing  properties  of  the  beam  amplification  provide  an 
unusual  method  for  generation  of  narrow  band  VLF  noise  pulses.  If  the 
beam  is  launched  from  its  mirror  point  toward  the  equator.it  is  possible 
for  VLF  noise  to  build  up  on  the  earthward  side  of  the  interaction 
region  and  subsequently  propagate  toward  the  ionosphere.  Since  the 
equatorward  side  of  the  interaction  region  is  absorbing  energy  more 
strongly,  no  signals  from  the  equatorial  region  can  penetrate  through 
the  beam.  Consequently,  for  any  signals  to  be  available  for  amplification 
the  electron  beam  and  the  natural  background  plasma  must  be  spontaneously 
emitting  copious  amounts  of  incoherent  VLF  noise.  Our  current  under¬ 
standing  of  this  inchoherent  noise  source  (Liemohn,  1965;  and  Trulsen 
and  Fejer,  1970)  suggests  that  it  is  wholly  inadequate.  For  example, 

in  Fig.  14  the  amplification  region  alone  can  provide  about  35  db  but 

-30 

the  incoherent  VLF  signal  is  only  10  watts/Hz  from  each  electron. 
Nevertheless,  it  is  intriguing  to  consider  the  possibility  of  generating 
a  succession  of  descending  tone  pulses  with  an  electron  beam. 


Low-Energy  Ion  Beams.  For  this  case,  the  attractive  feature  of 
injecting  cold  plasma  to  stimulate  amplification  by  catalytic  means 
is  coupled  with  the  advantages  for  direct  control  that  are  offered  by 
the  beam  concept.  In  the  magnetosphere  a  low-energy  (eV)  ion  is  subject 
to  the  forces  of  both  gravity  and  geomagnetic  gradients.  The  individual 
ions  are  reflected  at  low  altitude  by  the  geomagnetic  gradient  and  driven 
toward  the  equator  where  gravitational  forces  take  over  and  pull  them 
earthward  again.  Thus,  these  ions  are  caught  in  intra  hemisphere  orbits 
which  are  nonlinear  but  quite  repetitive.  Such  an  orbit  is  illustrated 
in  Figure  16  for  a  proton  with  energy  0.32  eV,  and  argon  ion  with  energy 
13  eV,  or  a  cesium  ion  with  energy  43  eV. 

Such  a  beam  does  not  interact  with  VLF  waves,  but  it  is  accompanied 
by  a  cloud  of  neutralizing  electrons  that  stimulate 
strong  amplification  by  the  hot  (keV)  electron  plasma.  The  generation  of 
such  an  ion  beam  by  a  particle  gur,  on  board  a  satellite  or  rocket  requires 
a  neutralizing  source  of  negative  charge  such  as  a  hot  cathode  emitting  an 
electron  current.  If  the  ions  follow  a  well-established  trajectory,  the 
electron  gas  will  be  dragged  along  with  the  icns  by  the  ambipolar  electric 
field.  It  is  this  cloud  of  electrons  that  actually  enhances  the  cold 
plasma  density  and  induces  catalytic  amplification  of  VLF  waves  along  the 
field-line  trajectory  of  the  ions.  The  shape  of  the  electron  cloud  is 
assumed  to  emulate  that  of  the  ion  beam  since  the  ambipolar  field  is  localized 
and  the  ions  are  so  much  more  massive. 

Before  undertaking  a  discussion  of  the  beam  interaction  with  VLF 
waves,  the  details  of  the  individual  ion  motion  in  the  gravitational  and 
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Fig.  16  Schematic  diagram  of  magnetic-gravitational  trapping  of  low 
energy  Ions  on  an  intra-hemisphere  trajectory.  Injection  of 
these  low-energy  ions  above  the  Ionosphere  can  produce  a 
sharp  density  enhancement  along  a  portion  of  the  flux  tube. 


magnetic  field  must  be  considered.  The  Important  forces  acting  on  the 
ion  are  the  gravitational  component  parallel  to  the  local  geomagnetic 
field  and  the  field-aligned  geomagnetic  gradient  force, 

fgrav  "  n1  9  s1n  */r2  (21) 

1  e  7»l  B  v.2 

FMAG  a  ~  2  1 - f*  <22) 

me  w 

c 

where  g  is  the  acceleration  of  gravity  at  the  surface  of  the  earth,  *  Is 
the  angle  between  the  radial  direction  and  the  field-line  normal,  and  v(| 
is  the  gradient  operator  along  the  local  geomagnetic  field  direction. 

These  force  components  coupled  with  the  first  adiabatic  Invariant, 

2 

VjYB  =  constant,  and  the  geometry  of  the  geomagnetic  field  are  sufficient 
to  fully  identify  the  Ion  trajectories.  However,  since  the  forces  are 
conservative.  It  Is  entirely  permissible  to  use  the  conservation  of  total 
energy  to  arrive  at  the  end  points  of  the  orbit  without  deriving  the  full 
solution  for  the  trajectory.  The  energy  equation  has  the  simple  form 

V»I  *  vfl  -  Wni  ■  »1h  - 


(23) 
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whoro  subscript  I  refers  to  the  Injection  location  and  M  refers  to  the 
mirror  or  turning  points  of  the  trajectory  along  the  flux  tube.  This 
formulation  of  the  problem  follows  that  used  previously  (Eviatar  et  al , 
1954;  and  LeMarie  and  Scherer,  1970),  but  the  results  are  displayed  in 
an  entirely  different  manner. 

The  solution  for  the  mirror  points  of  typical  mn  trajectories  are 
Hu  trated  n  Fig.  17.  The  Injection  level  Is  specified  at  1600  km 
<R  =  1.25)  along  the  geomagntlc  field  line  L  =  4.  The  mirror  point, 

is  plotted  as  a  function  of  Ion  energy  for  hydrogen,  argon,  and 
cebium  with  the  injection  pitch  angle  a  as  a  parameter.  For  example, 
an  argon  Ion  with  an  energy  of  13  ev  and  an  Injection  pitch  angle  of 
90  Brahes  back  and  forth  between  R  *  1.25  and  R  *  3  as  shown  in 
Fig.  16  as  well.  An  extreme  example  Is  a  4  eV  argon  Ion  Injected 
upward  parallel  to  the  local  field  (a  »  0°)  which  Is  reflected  at 
R  =  1.6  by  gravity  and  falls  directly  back  Into  the  atmosphere  at 
R  1  1.  Similar  cases  have  been  worked  out  for  a  variety  of  injection 
lev°l  ,  and  geomagnetic  field  lines;  all  of  them  have  very  similar  curves 
and  nearly  Identical  energy  regimes.  The  presence  of  a  natural  electric 
field  parallel  to  the  local  geomagnetic  field  can  seriously  distort  the 
shape  of  the  trapping  region,  however,  and  It  may  play  a  significant  role 
during  disturbed  conditions, but  It  must  be  Ignored  here  due  to  lack  of 
evidence. 

In  order  to  Illustrate  the  effectiveness  of  these  low-energy  Ion 
beams  to  stimulate  V! F  wave  amplification,  the  specific  case  of  43  eV 
cesium  ions  with  an  injection  pitch  angle  of  90°  at  the  rocket  altitude  of 
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Fig.  17  Trapping  zones  for  warm  Ion  injection  above  the  ionosphere. 
The  end  points  of  particle  orbits  for  the  various  ions  are 
displayed  as  a  function  of  energy.  The  injection  level  is 
at  an  altitude  of  1600  km  along  the  field  line  at  L  *  4. 


1600  km  will  be  considered  explicitly.  At  the  injection  location  these 
ions  have  a  speed  of  7.5  km/sec  and  a  gyro  radius  of  0.2  km.  Due  to 
scattering  mechanisms,  the  beam  helix  is  expected  to  expand  into 
a  column  about  1  km2  in  cross-sectional  area  which  is  10  times  the  gyro-area. 
If  the  beam  runs  continuously  for  100  seconds,  an  initial  column  length 
of  50-100  km  might  be  anticipated  since  most  of  tne  particle  speed  is 
perpendicular  to  the  field  line.  Thus,  an  initial  plasma  density  enhance- 
ment  of  500  cesium  ions/cm  may  be  achieved  with  a  beam  current  of  only 
0.05  amps  which  corresponds  to  a  gun  power  of  just  2  watts. 

Such  an  enhancement  of  cesium  ions  and  cold  electrons  drawn  up  by 

ambi polar  diffusion  is  not  particularly  effective  inside  the  plasmasphere 

but  is  a  remarkably  strong  enhancement  in  the  plasmatrough.  The  effect 

is  graphically  illustrated  in  Fig.  18  along  the  geomagnetic  field  line 

L  -  4  where  the  ambient  plasma  density  is  normalized  to  unity  at  the 

equator  and  varies  as  R  3.  Such  a  natural  density  corresponds  to  conditions 

during  and  immediately  following  a  geomagnetic  substorm.  Since  the  ion 

beam  is  confined  to  a  geomagnetic  flux  tube,  its  density  falls  off  only 
-3/2 

as  R  .  Consequently,  the  enhancement  is  even  more  effective  at  the 
higher  altitudes. 

Such  a  beam  enhancement  stimulates  amplification  of  VLF  waves  along 
its  field-line  path  at  frequencies  near  the  local  electron  gyrofrequency. 

The  local  amplification  exponent  for  selected  frequencies  is  displayed  in 
Fig.  19  along  the  field  line  L  =  4.  The  hot  electron  plasma  that  supplies 
the  energy  for  wave  amplification  is  assumed  to  have  the  form  E-2  sin2« 

(see  Figs.  1  and  2).  As  noted  in  the  figure,  the  region  of  enhancement 
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Fig.  18  The  net  plasma  density  enhancement  produced  by  a  warm  (43  eV) 
cesium  beam.  The  ambient  background  plasma  is  assumed  to  be 
in  the  plasmatrough  where  NQ  =  1  ion-electron/ cm3. 
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Fig.  19  Amplificatic  i  exponent  for  VLF  waves  propagating  through  the 
warm  (43  eV)  cesium  beam.  The  hot  electrons  are  assumed  to 
have  the  distribution  E”2  sin2",.  Note  these  frequencies  are 
above  the  equatorial  value  fc®. 
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covers  the  geomagnetic  latitude  band  from  25°  to  56°.  The  amplitude 
exponent,  ,  has  been  evaluated  for  frequencies  above  the  equatorial 
electron  gyro  frequency  fc®  but  the  actual  energy  exchange  between  the 
hot  electrons  and  VLF  waves  takes  place  below  the  local  electron  gyro 
frequency.  Since  f*Q  =  13.7  kHz,  the  amplification  actually  occurs 
in  the  LF  band  over  most  of  the  interaction  region. 

The  net  amplification  induced  by  the  cesium-electron  beam  along 
the  field  line  path  L  =  4  is  displayed  in  Fig.  20.  Amplifications  of 
30-40  db  may  be  expected  at  frequencies  between  60  and  150  kHz.  Near 
the  equatorial  zone  where  the  cyclotron  resonance  interaction  is  very 
strong,  amp1 ifications  of  100  db  or  more  apparently  can  be  expected  for 
frequencies  in  the  range  20-30  kHz.  One  should  be  reminded  that  this 
amplification  power  is  derived  catalytical ly  from  the  natural  hot 
electron  plasma  in  the  energy  range  10-100  keV;  it  is  not  derived  from 
the  2  watts  of  beam  power. 

Certain  properties  of  the  cesium  beam  make  it  rather  attractive  as 
a  stimulator  of  VLF  and  LF  rad’o  noise.  The  size  of  the  beam  column  is 
adequate  to  interact  with  an  appreciable  wave  front  of  VLF-LF  energy. 

At  50  kHz  the  free  space  wave  length  is  6  km  and  the  index  of  refraction 
is  10-30  so  that  the  effective  wavelength  is  significantly  less  than 
1  km.  The  gvroradius  for  cesium  ions  at  R  =  2-3  is  0. 6-1.0  km  and 
any  scattering  that  occurs  will  expand  the  beam  column  well  beyond  these 
values.  Since  a  column  of  50-100  kilometers  in  length  can  be  anticipated 
for  100  second  bursts  from  the  ion  gun,  there  is  an  adequate  interaction 
zone  extending  over  many  wavelengths.  A  single  plasma  beam  generates  a 
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Fig.  20  Net  path  amplification  for  VLF  waves  propagating  through  a 
warm  cesium  beam.  The  hot  electrons  are  assumed  to  have 
the  distribution  £-2  sin2a. 
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relatively  narrow  band  VLF-LF  descending  noise  tone  as  it  proceeds  up 
the  field  line  to  its  mirror  point  and  an  ascending  noise  tone  as  it 
returns  toward  the  injection  point.  Similar  tone  patterns  might  be 
generated  on  adjacent  field  lines  by  the  release  of  other  beam  columns. 
The  bandwidth  and  noise  amplitude  obviously  depends  on  the  geomagnetic 
field  line  and  altitude  of  the  injection  as  well  as  the  natural  hot 
electron  plasma  that  is  available  to  interact  with  the  waves. 
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ULF  AMPLIFICATION 

Stimulated  amplification  of  ULF  waves  in  the  magnetosphere  is  more 
difficult  than  its  VLF  counterpart  for  two  reasons.  First,  the  cyclotron 
resonance  interaction  is  inherently  weaker  at  ULF  because  the  index  of 
refraction  for  the  medium  is  relatively  smaller  requiring  a  relatively 
higher  value  of  the  resonance  velocity  (see  Eq.  (3)).  Second,  the 
injection  of  clouds  of  proton  plasma  which  would  provide  the  most 
efficient  means  for  catalytic  stimulation  of  the  interaction  are  simply 
not  feasible;  less  convenient  types  of  ion  injection  must  be  substituted. 
Despite  these  limitations,  significant  amounts  of  amplification  can  be 
stimulated  by  the  techniques  that  will  be  described  in  the  following 
subsections. 

For  ULF  wave  amplification  the  critical  consideration  is  the  pitch 
angle  and  energy  distributions  of  the  hot  protons  (MOO  keV).  Pitch 
angle  distributions  0f  the  form  sinmot  are  quite  applicable  to  the  proton 
plasma  as  well.  Distributions  with  m  =  1  or  2  are  frequently  encountered 
in  the  magnetosphere  and  will  be  used  in  the  analyses  that  follow. 

As  in  the  preceding  electron  cases,  the  energy  distributions  of  the 
protons  may  be  subdivided  into  a  cold  plasma  and  a  hot  plasma.  The  cold 
plasma  contains  most  of  the  charge  density  and  Its  distribution  defines 
the  propagation  characteristics  of  ULF  waves.  The  cold  protons  are 
assumed  to  have  energies  on  the  order  of  0.1  eV.  The  main  source  of  the 
hot  proton  plasma  is  the  ring  current  that  waxes  and  wanes  with  geo¬ 
magnetic  activity.  Some  typical  hot  proton  distributions  as  well  as 
the  cold  plasma  distribution  are  displayed  in  Fig.  21.  The  ring 
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Fig.  21  Typical  proton  energy  distributions  in  the  magnetosphere 
The  calculations  are  based  on  distributions  of  the  form 
E"n  which  are  commonly  encountered  by  spacecraft.  The 
ring  current  distribution  'B'  is  based  on  OGO-5  observa¬ 
tions  (Pizzella  and  Frank,  1971). 
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current  distribution  "B"  is  based  on  data  derived  from  several  satellite 

observations  (Pizzella  and  Frank,  1971)  beyond  the  plasmapause.  The 

distributions  E  n  have  also  been  used  extensively  for  quiet  periods  in 

the  magnetosphere.  The  excessively  high  values  of  these  distributions 

below  1  keV  does  not  play  a  significant  role  in  the  cyclotron  resonance 

interaction  since  the  resonance  energy  at  frequencies  of  interest  is  well 
above  this  level. 

Uthuvn  Jet  in  the  Plasmatrough.  During  strong  geomagnetic  storms 
the  plasmapause  frequently  comes  inside  L  =  3.5  and  the  proton  ring 
current  develops  a  strong  maximum  at  L  =  3.5  -  4.5  which  might  be  tapped 
rur  uLF  amplification.  In  this  case,  the  plasmatrough  has  an  unusually 
low  density  N  proton/cm3)  and  the  hot  proton  ring  current  is  exceptionally 
dense  (>_  0.02  protons/cm3).  Such  conditions  allow  artificial  stimulation 
of  ULF  amplification  by  injection  of  a  cloud  of  lithium  ions  (Cornwall 
a,id  Schulz,  1971).  One  source  for  lithium  clouds  at  L  >  3  -  4  (Hoch, 
private  communication)  is  a  shaped-charge  injection  of  a  neutral  lithium 
jet  from  a  low  altitude  rocket.  By  choosing  a  launch  time  and  location 
that  avoids  sunlight  for  part  of  the  ballistic  trajectory  to  the  equatorial 
region  at  L  =  3.5,  an  appreciable  fraction  of  the  initial  jet  might  be 
deposited  as  an  ion  cloud  in  the  plasmatrough  near  the  ring  current  maximum. 
The  process  is  illustrated  schematically  in  Figure  22.  Owing  to  the  ballistic 
trajectories  there  is  a  moderate  amount  of  focusing  that  can  be  achieved 
by  proper  injection  velocities  and  orientation  of  the  jet. 

A  detailed  model  of  the  injection  characteristics  has  been  developed 
in  order  to  assess  its  amplification  potential.  A  jet  velocity  band  of 
10-15  km/sec  may  be  assumed  at  the  source  (experimentally  observed  for 
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Fig.  22  Shaped-charge  injection  of  lithium  to  the  equatorial  region 
at  L  =  3.5.  The  neutral  lithium  jet  of  approximately  1  kgm 
is  expected  to  follow  a  ballistic  trajectory  across  the 
magnetic  field  until  solar  radiation  ionizes  it  and  deposits 
approximately  10  lithium-ion  electron  pairs/cm-3.  The  ambient 
background  plasma  is  assumed  to  be  in  the  plasmatrough  where 
Nq  =  1  proton-electron  pair/cm-3. 
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barium  injections  by  Peek,  private  communication).  Since  the  jet  is 
slowed  by  the  gravitational  forces  on  the  neutral  atoms,  it  is  estimated 
to  require  30-40  minutes  for  the  lithium  cloud  to  travel  about  2.5  earth 
radii  from  the  injection  point  to  the  equatorial  region  at  L  =  3.5. 

The  velocity  spread  of  5  km/sec  corresponds  to  a  radial  dispersion  of 
about  6,000  km.  If  the  jet  is  sprayed  in  a  narrow  range  of  latitudes 
and  is  carefully  columnated  along  the  geomagnetic  meridian  plane,  the 
lithium  ions  at  L  =  3.5  can  be  confined  to  a  latitude  band  of  10,000  km 
across  the  geomagnetic  equator  and  perhaps  only  2,000  km  spread  in 
longitude.  This  corresponds  to  a  volume  of  1026  cm3. 

Since  the  ionization  of  lithium  atoms  by  solar  ultraviolet  is  a 
continuous  process  with  a  characteristic  ionization  time  of  60  minutes, 
perhaps  only  10%  of  the  initial  jet  is  available  in  the  lithium  ion 
cloud.  Only  10  ions  of  lithium  are  available  in  1  kgm  so  that  100  kgm 
of  lithium  is  required  to  deposit  a  density  of  10  lithium  ions/cm3  in 
the  enhancement  region.  These  numbers  are  obviously  debatable,  but  they 
serve  as  a  basis  for  judging  the  usefulness  of  this  method.  Obviously, 
a  total  payload  weight  of  1,000  kgm  including  the  explosive  charge  is 

not  out  of  the  question  with  certain  rockets  that  might  be  considered  for 
this  experiment. 

According  to  equation  (4),  ULF  propagation  characteristics 
are  sharply  modified  in  a  lithium  ion  cloud.  Below  the  lithium  gyro- 
frequency  the  phase  velocity  is  reduced  appreciably  whereas  above  it, 
the  propagation  is  evanescent.  This  alteration  from  the  normal  proton 
plasma  propagation  characteristics  allows  the  cyclotron  resonance 


interaction  to  proceed  with  lower  energy  particles  which  are  normally 
more  numerous  and  consequently  the  Interaction  Is  appreciably  enhanced. 

In  accordance  with  Fig.  22  the  lithium  Injection  at  L  =  3.5  Is 
assumed  to  have  a  density  of  10  lithium  1ons/tm3  and  the  ambient 
plasmatrough  Is  assumed  to  have  a  density  N  «  Nq  (3.5/R)3  where  N0  » 

1  cm  3.  The  amplification  exponent,  ,  for  ULF  waves  propagating  through 
this  lithium  Ion  cloud  is  Illustrated  In  Fig.  23.  The  hot  proton  plasma 
for  the  ring  current  Is  assumed  to  have  the  distribution  E sln^  a  which 
is  consistent  with  experimental  observations  at  these  relatively  low 
L  shells.  The  curves  for  the  exponent  are  entirely  attributable  to  the 
lithium  cloud;  the  ambient  (cold)  proton  plasma  causes  a  negligible  amount 
of  natural  amplification  at  these  frequencies.  The  sharp  cutoff  In  the 
exponent  at  5°  and  15°  is  due  to  the  propagation  cutoff  above  the  local 
lithium  gyro  frequency.  Evidently,  there  Is  significant  non-equatorlal 
amplification  that  allows  strong  enhancements  at  frequencies  above  the 
equatorial  lithium  gyro  frequency  f£*  =  0.143  f  This  Is  more  vividly 
illustrated  In  Fig.  24  where  the  net  path  amplification  actually  has  Its 
maximum  above  fj^.  Obviously  ULF  waves  above  f£*  cannot  penetrate  the 
evanescent  zone  at  the  geomagnetic  equator,  so  that  the  ULF  band  from 
0.143-0.22  f^  Is  restricted  to  Intra hemisphere  propagation.  Despite  the 
high  lithium  ion  density  that  Is  proposed  for  injection,  the  net  ampliflca- 
t  on  Is  only  a  few  decibels.  Consequently,  this  method  for  stimulation 
of  ULF  amplification  Is  not  particularly  effective. 
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Fig.  23  Amplification  exponent  for  ULF  waves  propagating  through  the 

lithium  shaped-charge  injection  along  L  *  3.5.  The  hot  protons 
ere  assumed  to  have  the  distribution  E-2  sin2a  (see  Figs.  1 
and  21 ). 
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24  Net  path  amplification  for  ULF  waves  propagating  through  the 
lithium  shaped-charge  injection  along  L  =  3.5.  The  hot  protons 
are  assumed  to  have  the  distribution  E_2  sin^a  ('ee  Figs.  1 
and  21). 
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GeoeiinchronouB  Lithium  Cloude.  The  Injection  of  lithium  Ion 
clouds  from  satellites  at  geosynchronous  altitudes  has  already  been 
discussed  by  others  (Cornwall  and  Shultz,  1971;  and  Cornwall,  1972), 
but  it  Is  Included  here  in  order  to  corroborate  their  findings  and 
perhaps  elaborate  upon  them.  The  amplification  analysis  performed 
here  Is  based  on  the  two  cloud  models  depicted  In  Fig.  8.  Details  of 
the  injection  requirements  and  diffusion  properties  are  described  In 
the  section  on  VLF  amplification.  In  one  case  the  1  kgm  cloud  of 
lithium  Ions  Is  confined  to  a  radius  of  0.1  Rg  which  corresponds  to 
10  minutes  of  diffusion  from  a  point  source,  whereas  the  other  case 
has  a  radius  0.5  Rg  corresponding  to  50  minutes  of  diffusion. 

As  before.  It  Is  Important  to  recognize  that  the  lithium 
Injection  alters  the  local  propagation  characteristics.  Above  the 
lithium  gyrofrequcncy,  ULF  waves  are  evanescent  whereas  below  the  phase 
velocity  Is  strongly  reduced.  The  latter  effect  gives  rise  to  the 
enhanced  cyclotron  resonance  Interaction. 

The  ULF  amplification  that  Is  stimulated  by  these  two  Injection 
models  Is  Illustrated  In  Figs.  25-28  for  two  models  of  the  hot  proton 
distribution  function.  The  model  E  sin  a  represents  quiescent 
conditions  In  the  magnetosphere  when  the  proton  ring  current  Is  very 
weak.  The  second  case  is  the  ring  current  model  "B"  which  has  been 
experimentally  observed  during  geomagnetic  storm  conditions;  It  Is 
denoted  by  the  functional  form  (dN/dE ) .,g,,  sin  a  (see  Fig.  21}. 

As  expected,  the  ULF  amplification  enhancement  Is  restricted  to 
the  narrow  equatorial  region  where  the  lithium  Injection  dominates  the 
total  plasma  density.  The  Injection  region  is  restricted  to 
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Fig.  25  Amplification  exponent  for  ULF  waves  propagating  through  the 
geosynchronous  lithium  injection  (see  Fig.  8).  The  exponent 
is  sharply  peaked  In  the  equatorial  enhancement  which  nave 
m^an  diameters  of  1280  and  64  km.  The  hot  protons  are  assumed 
to  tave  *he  distribution  E-l*5  sin  a. 
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Fig.  26  Net  path  amplification  for  ULF  waves  propagating  through  a 
geosynchronous  lithium  injection  (see  Fig.  8).  The  enhance¬ 
ment  is  entirely  attributable  to  the  lithium  bubble.  The  hot 
protons  are  assumed  to  have  the  distribution  E-' sin  a. 
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Fig.  27  Arpl ification  exponent  for  ULF  waves  propagating  through  a  geo¬ 
synchronous  lithium  injection  (see  Fig.  8).  The  exponent  is 
sharply  peaked  in  the  equatorial  enhancement  which  have  mean 
diswlers  of  1280  and  6-00  km.  The  hot  proton  distribution  Is 
based  on  Model  B  for  the  ring  current  of  a  typical  geomagneti 
substorm  (see  Fig.  21). 


NET  PATH  AMPLIFICATION  (db) 


74. 


Fig.  28  Net  path  amplification  for  ULF  waves  propagating  through  a 

geosynchronous  lithium  injection.  The  enhancement  s  entirely 
attributable  to  the  lithium  bubble.  The  hot  proton  distribu¬ 
tion  is  based  on  Model  * B '  for  the  ring  current  of  a  typical 
geomagnetic  substorm  (see  Fig.  21). 


±  3°  geomagnetic  for  the  10  minute  case  and  +  5°  for  the  50  minute  case. 
ci nee  the  interaction  Is  dominated  by  the  lithium  gyrofrequency  cutoff  in 
the  propagation,  ULF  frequencies  above  0.143  fp  simply  do  not  get 
amplified  significantly.  An  important  distinction  betv/een  the  results 
in  Figures  25  and  27  is  the  order  of  magnitude  Increase  in  the  value  of 
k.  during  disturbed  conditions. 

The  dependence  on  the  shape  of  the  energetic  particle  distribution 
nanifests  itself  more  explicitly  In  the  net  path  amplification  Illustrated 
in  Figures  26  and  28.  The  quiescent  conditions  yield  an  amplification  of 
4-6  db  whereas  the  ring  current  conditions  provide  35-55  db.  This  remarkabl 
difference  is  attributable  to  the  higher  value  of  the  ring  current  density 
(Figure  21)  at  those  energies  where  the  resonance  energy  exchange  Is  most 
effective  (2  -  200  keV).  Another  Important  conclusion  is  that  the 
amplification  Is  stronger  at  50  minutes  after  Injection  than  at  10  minutes. 
This  is  attributable  to  the  longer  propagation  path  over  which  amplification 
is  enhanced. 

The  amplification  achieved  by  injection  of  a  lithium  Ion  cloud  at 

geosynchronous  altitude  is  comparable  to  the  earlier  predictions  (Cornwall 

and  Schultz,  1971).  Perhaps  what  has  not  jeen  emphasized  as  strongly  by 

previous  research  is  the  relatively  sharp  amplification  maximum  that  Is 

nearly  centered  on  the  equatorial  lithium  gyrofrequency  fLi  which  has 

co 

the  value  0.0  Hz  at  L  =  6.6.  The  duration  of  the  lithium-ion  cloud 
depends  of  course  on  the  strength  of  the  local  electric  fields  in  the 
magnetosphere,  but  it  Is  expected  to  survive  for  several  hours  as  a 
discernible  entity.  The  ULF  noise  generated  bv  the  cl  cud  will  also  have 
c  distinctive  characteristic.  According  to  the  results  In  Figure  28, 
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the  noise  consists  of  an  ascending  tone  of  Increasing  amplitude 
maximizing  at  f|:J  along  the  Injection  field  line  and  subsequently 
broadening  over  a  narrow  range  of  frequencies  corresponding  to  the 
L  shells  covered  by  the  lithium  cloud.  If  a  convective  electric 
field  is  present,  the  cloud  may  drift  radially  causing  a  change 
in  the  frequency  corresponding  to  the  change  in  the  local  gyrofrequer.cy. 

Energetic  Proton  Berne .  Amplification  of  ULF  waves  by  direct 
injection  of  energetic  proton  beams  Is  obviously  more  desirable  than 
the  preceding  methods  because  It  Is  not  subject  to  the  vagaries  of 
geomagnetic  activity.  The  wave  growth  that  Is  achieved  by  beam 
injection  Is  derived  entirely  from  the  energy  of  the  beam  particles. 

The  process  is  not  catalytic  In  the  sense  that  the  natural  background 
of  hot  protons  In  the  magnetosphere  has  a  negligible  effect.  The 
amount  of  amplification  is  proportional  to  the  Intensity  of  the  beam 
as  expected  and  the  frequency  band  depends  on  the  beam  pitch  angle  at 
the  Injection  location. 

The  beam  is  assumed  to  be  launched  during  geomagnetical ly  quiet 
conditions  along  the  field  line  at  L  =  4  which  is  assumed  to  be  Inside 
the  plasmasphere  where  NQ  =*  232  protons/cm3.  For  completeness,  the 
natural  hot  proton  distribution  is  assumed  to  have  the  form  E"1'5  sin  a, 
although  it  does  not  contribute  appreciable  amplification. 

The  proton  gun  Is  arbitrarily  assumed  to  have  a  power  level  of 
10  kilowatts,  and  the  beam  current  consists  of  100  millisecond  pulses 
of  10  keV  protons.  In  so  far  as  the  author  is  aware,  no  one  has  developed 
a  proton  gun  for  use  on  spacecraft.  The  foregoing  gun  characteristics 


are  probably  very  difficult  to  achieve  due  to  the  high  current  level 
(in  amps).  These  beam  parameters  were  chosen  intentionally  so  that  the 
proton  beam  could  be  modelled  quantitatively  by  the  same  pitch  angle  and 
energy  distributions  that  are  used  for  the  electron  beam  (see  Fig.  13). 
The  protons  have  a  speed  of  1.5  x  103  km/sec  and  an  equatorial  gyro- 
radius  of  approximately  30  km,  depending  as  before  on  the  local  pitch 
angle.  Electrostatic  forces  and  plasma  instabilities  undoubtedly  break 
up  the  initial  helix  trajectory  so  that  the  beam  is  smeared  into  a 
plasma  column.  Assumi,  j  an  order  of  magnitude  inflation  in  the  cross 
section  of  the  beam,  the  column  volume  is  -v  4  x  106  km3.  In  each 
pulse  there  are  1.6  x  1019  protons,  so  that  the  density  of  the  beam  is 
0.004  protons/cm  .  Thus,  for  an  effective  beam  energy  spread  of  2  keV, 

the  differential  energy  spectrum  is  0.002  protons/cm3  keV  (as  shown 
in  Fig.  13). 

The  pitch  angle  distribution  HQ  is  normalized  to  unity  with  a 
spread  of  10“  corresponding  to  a  limited  amount  of  scattering.  As  in 
the  electron  case,  the  pitch  angle  distributions  are  centered  on  85 
55  and  25  corresponding  to  mirror  points  at  +3°,  +18°,  and  +37°, 
respectively. 

The  effectiveness  of  the  beam  for  stimulation  of  the  cyclotron 
resonance  interaction  depends  on  its  diameter  and  length  relative  to 
the  ULF  wavelengths  involved  in  the  interaction.  Near  the  geomagnetic 
equator  at  L  =  4  the  ULF  index  of  refraction  is  500-600  for  propagation 
frequencies  around  3Hz  which  corresponds  to  0.5  fPQ.  The  corresponding 
wavelength  at  3  Hz  is  160-200  km  which  is  comparable  to  the  dimensions 
of  the  plasma  beam  column.  Since  the  beam  should  be  several  wavelengths 
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in  each  dimension,  its  effectiveness  as  an  amplifier  is  very  questionable 
even  at  the  large  beam  power  levels  proposed.  Nevertheless,  It  is  worth¬ 
while  to  investigate  the  amplification  potential  assuming  that  ways  may 
be  found  to  increase  the  size  of  the  beam  while  maintaining  its  internal 
density. 

The  amplification  exponent  for  a  beam  with  aSTREAM  =  55  is  shown 
in  Fig.  29.  The  sharp  peaks  In  k1  are  due  to  the  large  gradients  in  the 
proton  beam  energy  and  pitch  angle  distributions  which  cause  strong 
interactions  in  limited  regions  where  the  beam  energy  happens  to  match 
the  local  cyclotron-resonance  Interaction  energy.  On  the  earthward 
side  of  the  interaction  region,  the  beam  strongly  amplifies  ULF  waves 
whereas  on  the  equatorward  side  It  strongly  absorbs.  Curiously  the 
frequency  at  0.4  is  strongly  amplified  around  the  geomagnetic 
equatorial  region;  evidently,  the  resonance  energy  never  gets  low  enough 
to  permit  significant  absorption. 

This  result  is  more  clearly  demonstrated  In  Fig.  30  where  the  net 

path  amplification  at  0.4  fcP  has  a  sharp  peak  of  50-60  db  for  the  stream 

at  55°.  There  are  other  unexpected  results  as  well.  The  55°  has  another 

maximum  around  0.9  fcP  which  must  be  attributed  to  the  local  interaction 

characteristics  at  a  geomagnetic  latitude  around  20°.  The  beam  with 
STREAM.  Q  0  _ 

°o  ‘  85  has  a  sharP  maximum  at  0.7  f^  which  Is  attributable  to  a 
narrow  equatorial  band  Interaction  where  the  absorption  portion  of  the 
distribution  never  has  an  opportunity  to  enter  the  interaction.  Finally, 
the  beam  at  25°  is  simply  too  weakly  interacting  in  the  equatorial  region 
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Fig.  29  Amplification  exponent  for  ULF  waves  propagating  through  a 
proton  beam  along  the  field  line  at  L  =  4  inside  the  plasma- 
pause.  The  beam  model  has  the  same  distributions  as  for 
electrons  (§ee  Fig.  13).  The  beam  has  an  equatorial  pitch 
angle  of  55  and  a  particle  beam  energy  of  10  keV. 
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Fig.  30  Net  path  amplification  for  ULF  wave:  propagating  through  a 
proton  beam  along  the  field  line  at  L  =  4  inside  the  plasma- 
pause.  Equatorial  pitch  angles  for  the  beam  are  25°,  55°, 
and  85°  and  the  beam  energy  is  10  keV. 


to  provide  much  net  exchange  of  energy.  The  reader  Is  referred  to  the 
section  on  the  electron  beam  anpllflcatlon  of  VLF  waves  for  more  details 
about  the  general  physics  of  the  Interaction. 

If  proton  beams  of  sufficient  energy  can  be  generated  from  rocket 
or  satellite  platforms,  they  may  prove  to  be  useful  generators  of  ULF 
noise.  From  the  cursory  results  displayed  In  Figs.  29  and  30,  it  seems 
clear  that  certain  noise  frequencies  will  be  st'^ongly  amplified  depending 
on  the  source  location  and  Injection  pitch  angle.  Unlike  the  electron 
beam-VLF  noise  Interaction  which  provided  little  net  amplification, 
the  proton  beam  ULF  noise  Interaction  yields  a  net  amplification  that 
may  produce  a  recognizable  signature.  Much  more  research  Is  needed  to 
verify  these  results. 


CONCLUSIONS 


The  foregoing  methods  for  enhancement  of  the  cyclotron-resonance 
interaction  to  amplify  VLF  and  ULF  waves  in  the  magnetosphere  were  chosen 
for  analysis  because  they  showed  potential  for  stimulation  of  amplification. 
Generally  the  models  for  the  interaction  are  based  on  current  state-of-the- 
art  knowledge  about  injection  methods  and  charged  particle  guns.  Conditions 
in  the  magnetosphere  were  arbitrarily  chosen  to  represent  either  a  disturbed 
or  a  quiescent  state  with  an  arbitrary  ambient  plasma  density.  The  parameters 
are  not  necessarily  the  most  likely  conditions  that  will  be  encountered, 
but  they  are  typical  and  provide  a  uniform  basis  for  comparison  of  the 
different  methods  for  stimulation  of  the  cyclotron  resonance  instability. 

A  much  more  thorough  investigation  of  those  methods  that  show  promise  is 
advisable  in  order  to  assess  the  variety  of  results  that  can  be  anticipated 
for  a  wide  range  of  conditions. 

Perhaps  the  most  immediate  question  concerns  the  validity  of 
the  linear  theory  in  these  applications.  It  clearly  is  a  useful  start¬ 
ing  point  for  these  investigations  because  a  wide  variety  of  injection 
methods  may  be  studied  parametrically  for  a  relatively  small  expend¬ 
iture  of  computer  time.  However,  most  of  the  useful  results  do 
penetrate  the  bounds  of  linear  validity  in  some  way.  There  are  at 
least  three  levels  of  nonlinear  sophistication  that  may  be  applied. 

First,  the  Taylor  series  expansion  used  to  derive  equations  (4)  and 
(5)  requires  |  M  «kr  which  is  frequently  violated  in  the  foregoing 
applications.  An  asymptotic  solution  for  |kj  »  kp  is  readily  derived 
and  numerical  solutions  of  (1)  are  possible  when  |k^|  =  kr.  Numerical 
solutions  for  the  acoustic  mode  (Liemohn  and  Scarf,  1964)  suggest  that 


the  qualitative  behavior  for  |k.|  *  kf  is  described  by  the  linear 
expressions  above,  but  the  quantitative  results  are  open  to  question. 

Second,  the  linear  theory  is  based  on  the  assumption  that  the 
phase  space  distribution  FHQT  is  not  perturbed  significantly  by  the 
electromagnetic  fields.  The  extent  of  this  distortion  depends  on  the 
magnitude  of  the  fields  and  the  repetitive  interaction  rate.  These 
effects  are  difficult  to  estimate  quantitatively,  because  the  field 
amplitudes  are  net  well  known  and  the  repetition  depends  on  the  size  of 
the  injection  and  the  local  resonance  energy.  The  quasi-linear  theory 
describes  the  equilibrium  distortion  of  FHQT  rather  well  for  the  natural 
conditions  in  the  magnetosphere,  but  it  is  inappropriate  for  these 
transient  applications.  The  conclusion  is  that  a  full-blown  nonlinear 
solution  is  needed  to  properly  assess  the  linear  results.  Hopefully, 
existing  nonlinear  theory  (Roux  and  Solomon,  1971;  Bud'ko  et  al . ,  1972; 
Istomin  and  Karpman,  1972  a,  b;  and  Brinca,  1972)  can  be  adapted  to 
these  applications. 

Based  on  the  linear  analysis,  a  few  cases  show  considerable 
promise  for  stimulation  of  significant  amplification;  one  or  two  may 
be  appropriate  for  further  study  if  a  breakthrough  in  injection  methods 
can  be  achieved,  and  some  cases  simply  did  not  have  the  desired  proper¬ 
ties.  For  example,  at  VLF  the  barium  shaped-charge  injection  gave 
a  moderate  mount  of  amplification  only  after  30  minutes  had  passed  and 
its  effectiveness  was  further  limited  by  the  experimentally  observed 
striations  that  separated  the  jet  of  particles.  Consequently,  it  appears 
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to  be  a  relatively  weak  candidate  using  the  current  technology  that  was 
modeled  in  the  case  presented  above.  However,  the  advent  of  a  larger 
payload  at  higher  injection  velocities  which  may  produce  a  denser  cloud, 
does  offer  real  promise  of  strong  amplification. 

The  geosynchronous  lithium  injection  is  currently  the  most  promising 
method  for  stimulating  either  VLF  or  ULF  radio  noise.  The  conditions 
modeled  some  tO-minutes  after  injection  are  unrealistic  at  the  present 
time  due  to  the  large  payload  required  to  achieve  the  plasma  density; 
however,  the  conditions  modeled  at  EO  minutes  after  injection  appear 
to  be  readily  achievable  with  current  technology  and  a  modest  payload. 

These  cases  do  require  a  nonlinear  investigation,  however,  because 
ki  =  kr  when  amplification  is  appreciable  (above  20-30  db). 

The  electron  and  proton  beams  have  attractive  characteristics 
for  communication  applications.  But  they  do  not  appear  to  be  feasible 
due  to  beam-current  limitations.  At  VLF  the  electron  beam  has  the 
property  that  it  produces  as  much  absorption  as  amplification  so  net 
wave  growth  is  negligible.  At  ULF  the  net  amplification  appears  to 
be  much  stronger,  but  for  significant  amplification  a  major  breakthrough 
is  required  for  the  proton  gun.  Thus  neither  type  of  beam  looks  promising. 
Furthermore,  no  other  parameter  regimes  look  hopeful. 

The  low-energy  cesium  beam  shows  considerable  potenti'l  both  for 
stimulation  of  amplification  and  as  a  flux  tube  for  earthward  propagation. 
This  somewhat  radical  concept  should  not  be  difficult  to  test  because  the 
technology  is  well  established  in  laboratory  plasma  systems.  The  beam 
of  electrons  that  is  assumed  to  follow  the  cesium  is  of  course  somewhat 


uncertain,  but  ambipolar  diffusion  in  the  plasma  media  is  well  recognized. 

An  interesting  sidelight  is  that  the  amplification  exponent  in  Fig.  19 
satisfies  the  requirements  for  application  of  the  linear  theory.  The 
only  remaining  question  is  whether  the  helical  beam  expands  into  a  larger 
plasma  column  as  assumed  in  the  analysis. 

The  lithium  injection  at  L  =  3.5  that  is  proposed  to  simulate 
ULF  waves  fails  miserably  but  it  has  important  consequences  for  under¬ 
standing  the  natural  interaction  at  the  plasmapause.  Ignoring  the  injection 
for  the  moment,  the  important  quantity  is  the  index  of  refraction  just 
inside  and  outside  the  plasmapause.  Inside  the  index  is  at  least  an 
order  of  magnitude  greater  than  outside  because  the  plasma  density  is 
correspondingly  lower  in  the  plasmatrough.  Thus  the  resonance  speed  f0r 

the  cyclotron  resonance  interaction  (equation  (7)  is  an  order  of  magnitude 
larger  outside  than  inside  the  plasmapause,  and  the  resonance  energy  is 
two  orders  of  magnitude  greater.  Consequently,  there  is  no  appreciable 
amplification  in  the  plasmatrough  immediately  outside  the  plasmapause  and 
the  proton  ring  current  is  relatively  unaffected  until  it  diffuses 
radially  inward  across  this  discontinuity  (Cornwall  et  al.,  971).  Returning 
to  the  lithium  injection,  it  is  now  clear  that  huge  densities  of  lithium  are 
necessary  to  induce  significant  cyclotron  resonance  interactions  below 
the  local  lithium  gyro  frequency  just  outside  the  plasmapause.  The  best 
place  to  inject  lithium  ion  clouds  for  stimulation  of  amplification  is 
at  geosynchronous  altitudes  where  only  modest  amounts  of  lithium  are 
required  and  the  proton  ring  current  is  still  adequate. 

Future  studies  must  investigate  the  nonlinear  aspects  of  the 
cyclotron  resonance  interaction  for  those  cases  that  show  promise  of 
stimulating  significant  amplification.  At  VLF,  improvements  in  the 
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shaped-charge  bariun  jet  must  be  investigated  theoretically  for  its 
nonlinear  limitations,  and  the  low-energy  cesium  beam  results  should  be 
independently  corroborated.  At  ULF  the  only  case  that  appears  realisti¬ 
cally  feasible  is  geosynchronous  Iniection  of  a  lithium  cloud.  While 
these  cases  may  not  fulfill  the  immediate  objective  for  communications, 
they  do  represent  demonstrable  ways  to  test  the  cyclotron-resonance 
interaction  quantitatively.  When  this  has  been  verified  more  exotic 
means  of  transmitting  information  can  be  explored. 
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CYCLOTRON  RESONANCE  AMPLIFICATION  OF 
ULF  AND  VLF  WHISTLER-MODE  WAVES  IN  THE  MAGNETOSPHERE 

Computer  Program  Description 


Programmer:  R.  J.  Kerr 
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Cyclotron  Resonance  Amplification  of  ULF  and  VLF  Whistler-Mode  Waves  in  the 
Magnetosphere. 


Mathematical  Models 

Net  Amplitude  Exponent  (dimensionless) 


XNIXY  -  n. 


imag 


ckj/w 


For  other  models,  see  text. 


cyciocron  Kesonance  ampiincation  ot  uu*  ana  VLh  wnistler-Mode  Waves  in  the 
Magnetosphere. 


INPITT  VALUES 


Card  No. 

Columns 

Format 

Symbol 

Definition 

1 

1-2 

12 

MONTH 

Month,  i.e.,  February  enter  02. 

1 

3-4 

12 

NH\Y 

Day  of  month. 

1 

5-6 

12 

NYR 

Year,  i.e.,  1972  enter  72 

2 

1  80 

80U 

-- 

Identification  title  for  the  data  set 
imwediately  following. 

3 

MO 

F10.1 

TH0 

Surface  latitude  (deg). 

3 

11-20 

F10.0 

XNO 

Equatorial  density  (cm**). 

3 

21-25 

15 

NTYPEQ 

Indicator  for  cold  plasma  density 
model.  See  text  for  list  of 
available  models. 

3 

26-30 

FS.O 

YS 

Lower  band  limit  of  range  of  propagation 
frequency/equatorial  cyclotron  frequency. 

3 

31-35 

FS.O 

YE 

Upper  band  limit  (as  above). 

3 

36-40 

FS.O 

DELY 

Step  size  (as  above). 

3 

41-45 

IS 

NALPHA 

Number  of  integration  steps  over  pitch 
angle  range. 

3 

46-50 

15 

NX 

Number  of  integration  steps  over 
latitude  range. 

3 

51-55 

15 

NBETAP 

Number  of  integration  steps  over  per¬ 
pendicular  velocity  range. 

3 

56-50 

Frc0 

PORE 

Indicator  PORE  ■  1  for  electron  plasma 
particles.  PORE  *  1836  for  ion  plasma 
particles. 

4 

1-5 

15 

NUMUQ 

Number  of  cards  to  be  read  for  the 
quiescent  pitch  angle  and  distribution 
values. 

4 

6-10 

15 

NIMJQ 

Number  of  cards  to  be  read  for  the 
quiescent  normal  energy  and  distribution 
values. 

4 

11-15 

15 

NIMB 

Number  of  cards  to  be  read  for  the 
streaming  normal  energy  and  distribu¬ 
tion  values. 

Magnetosphere. 

Input  Values,  cont. 


Card  No. 

Columns 

Format 

Symbol 

Definition 

1 

10-20 

15 

NIMJS 

Number  of  cards  to  be  read  for  the 
streaming  normal  energy  and  distribution 
values . 

s 

1-10 

F10.3 

ALPHAQ(1) 

First  value  of  quiescent  pitch  angle 
array  (deg) . 

5 

11-25 

F15.3 

“»(i) 

First  value  of  distribution  array 
(relative) . 

1  *  \IMIQ 

1-10 

F10.3 

ALPHAQ 

(NUMHQ) 

NUMHQ  value  of  quiescent  pitch  angle 
array. 

l  *  NUMHQ 

11-25 

F15.3 

^(NUNtlQ) 

NlMiQ  value  of  quiescent  distribution 
array. 

5  ♦  NUMIQ 

1-10 

F10.3 

^(1) 

First  value  of  quiescent  normal  energy 
array  (LOG^  (keV). 

S  ♦  NUMHQ 

• 

• 

11-25 

F15.3 

^(i) 

First  value  of  distribution  array 
(LOG10  (cm'3). 

• 

• 

4  ♦  NUMHQ 

♦  MJMNQ 

1-10 

F10.3 

(NUMNQ) 

• 

NUMNQ  value  of  quiescent  normal  energy 
array. 

4  ♦  NUMHQ 

♦  NUMNQ 

11-25 

F15.3 

^ (NUMNQ) 

NUMNQ  value  of  quiescent  distribution 
array. 

It  values  tor  streaming  arrays  are  to  be  read,  they  are  punched  in  the  same 
manner  as  given  for  the  quiescent  arrays  and  follow  immediately  behind.  Values 
for  pitch  angle  and  normalized  energy  must  be  in  ascending  order. 

lo  run  more  than  one  set  of  data,  repeat  instructions  starting  at  Card  #2. 
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SAMPLE  OUTPUT 


1=1.  J  N  =  1.6  tLtCTKON  TEST 

jU  <  FACl  LATITUOt  (THg)  =  bj.  Ju  PLASMA  PARTICLE  (PORt)  =  ~1 

SuUhIuRaAl  JlNSaTY  (  XNu)  =  232.  gU  ULNS1TY  MODEL  (NTYPLQ)  =  2 

PRJPAGATIDn  FRtQ. /eQUATORIAL  CYCLOTRON  FREQ.  <!  OELY )  =  .100 
UAQ  LI11TS  (YS)  =  ,1JJ  (YL)  =  •  9  0  U 


gUItSCLNT  01 STR1BUTIUNS 


PITCH  ANGlE 

distribution 

( ML  PH  A  Q ) 

(HQ) 

1  0  .  w  u  0 

\  li  D  ‘ 

1 2  .  C  v  0 

•  ii  o  j 

It •  0  J  0 

.  I5u 

lo  lUtii 

.21l 

1 3  •  0  J  0 

.  2oG 

20.  L  j|} 

.31, 

2  2 .  U  0  0 

.  3yC 

Et.OJG 

.390 

c  o .  0  j  0 

«  h  3C 

23.00c 

.h7J 

30 • u  g  l 

.5  00 

32. U  g  0 

•  6  30 

■J  4  .  L  J  0 

.  y  59 

3  O  »  Cob 

•  58o 

3  3  .  u  0  u 

•  b  lb 

iiii  GliO 

.  b*t3 

1 2  «  u  ^  C 

•  bo9 

‘*4.0.0 

•  b  9y 

to. 0  G  U 

.  719 

1 8  •  0  u  u 

•  7  +  3 

y  J  •  0  J  0 

,7  ot 

5  2 . oil 

•  I  03 

yt . L  G  0 

.60  0 

*30 . 0  u  0 

.  32  1 

63  ,  l uG 

,  3 .5 

oO . o  w  o 

.  Sob 

G2.0tu 

.833 

0  4  t  (J  u  U 

.  899 

oOtCwO 

.  9  it 

o  3 .  w  u  0 

.927 

7  6 .  j  u  u 

,94b 

It, C JO 

.  951 

Ik, 0 0  u 

.  961 

7b. OuC 

.  970 

/  6 ,  u  o  u 

.978 

oU.UGO 

.  985 

o  2  •  G  G  0 

.  990 

oh. o  t  u 

.994 

OO  •  boil 

.  998 

^  i  t  u  v  l 

.9  09 

^  0  •  L  J  O 

1.  UOG 

NORM.  ENERGY 

Distribution 

(LUG  (c. /  MCQ)  ) 

( LOG ( DNQ/ ( OE/MCQ)  )  ) 

"5. 0  00 

4.750 

-4.  5  Lb 

4.000 
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ELECTRON  BEAM 
a0  =55°  EQ  =  10  keV 
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PROPAG.  F-etJ./tUTRL.  CYCLOTRON  FRtQ.  (LOCAL )  =  *3000 
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4*  7  3+02+02 

4..1  +  88E-G3 

-3.2265E-Q5 

4jl.7_865E-(L2 

•  0433 

1.3  192c  +  jj 

+  .  6729E+G2 

3.2387c  “0  3 

-3.8732E-05 

3.7767L-02 

•  0  65b 

1.  03432+ JJ 

4. 58822+02 

“1  .8+302-03 

-3.9715E-05 

-2.26872-02 

*  o  866 

1*  0541E  +  0  j 

4.  48262+02 

-3 .3o+b2-G  3 

-4.6085E-05 

-4. 22042-02 

.1083 

1  *L  JZ87.E 

4*  3§8.3_E+42_ 

-6.6772c-04 

-5.0457E-05 

-9.1/66t-03 

_*ma 

1.1085c+Q J 

4.  21782+02 

“1 .58962-05 

-5.43782-05 

-9. 32262-04 

.1516 

1.14392+JJ 

4.  C638E+02 

•2 *635+2-06 

-5. 6.282 1- Q5 

-8*07352-04. 

*  1 732 

1*  18542  +  J  J 

3.  89882+02 

0. 

-6. 2214E-05 

-9.O362E-04 

.1949 

1.23352+JJ 

3.  7255E+U2 

u  • 

-6*  3036E-05 

-9.6519E-04 

.2165 

i. 2  38  9E  +  3  J 

3. 54632+02 

0. 

-6.31+7E-05 

-1*  0241E-0  3 

•  2382 

1*  3.524c +  0J 

3.  3fi35Jt+li2 

JU._  . 

-**1553E-Q5 

-1*  0b22£-03 

.2598 

le  42512+J J 

3.1791E+02 

0. 

-5*6b4bE-05 

-1.04492-03 

.2815 

1*  5081c  *  J  J 

2.  99502+02 

u  • 

4277.2-05 

-1*67482-03 

j_3  031 

l.o029£+ Jj 

2  •  8 127  2+  02 

0  • 

-4.9bb7E-05 

-1.06052-03 

•  3248 

1*  7 1132+ JJ 

2.  63362+02 

0. 

-4. 35212-05 

-l.GGb3t-Q3 

•  3464 

1* 83542+ j  J 

2  •  4587  2+  02 

0  • 

-3*  91962-05 

-9.85742-04 

.3681 

1.977  71  +  J  J 

2*28892+02 

0  • 

-3j>4397E-05 

_ -9*44952-04 

.  jl.3897  _ 

2. 14142  +  OJ 

c.  12492+02 

J . 

-2.9170E-05 

-8. 7927E-04 

.4114 

2.3  3Q5E  +  TJ 

1.  96732  +  02 

0. 

-2.53b0E-05 

-0*425.32-0-4 

*4330 

2. 54982  + J  J 

1.  8 165 E+  02 

J . 

-2.1449E-05 

-7. 89292-04 

.4547 

2. 3  G54t  +  J  J 

1.  67272+02 

G  t 

-1.73bbE-05 

-7.11292-04 

.4763 

3.  1G5  02  +  J  J 

1. 5360  2+02 

0. 

-1.4007E-05 

-6*  4199E-04 

.4980 

3.  +584C  ♦jjl  J 

1  •  4  057 £+ 02 

G  • 

-1.1228E-05 

-5*79132-04 

.5196 

3.87832+ jj 

1 . 2  8+6  2+02 

u  • 

-8 • 41842-06 

-4.91582-04 

.5413 

*♦.  38122  +  1  J 

i.  1698t+w2 

G  f 

-5.9+41E-06 

-3.«  3557.2-0  4 

.  .5219 

4. 96892 +  J  J 

1.  L 6222  +  02 

0  . 

-3.92o8E-0b 

-3.00132-04 

.5846 

5 . 7  33  5c + JJ 

5.6173  2+Gl 

0  . 

-2.3847E-06 

-2*10772-04 

•  6Q62 

o.o  +  59E  +  'Jj 

b  .  6  8 17  2+  u  1 

0  . 

-1. 230uE-0b 

-1.2b892-04 

.6279 

7.7.90JE  +  Jj 

7  .  81  +  5  fc  +  y.l 

J... 

-4..672bE-17 

-5.31552^05  ... 

.6495 

9. 24G52  +  J  J 

7, 0142  c+  01 

u  • 

-1. 1588E-07 

-1. b79b2-05 

.6712 

1. 11092+ Jl 

b .  2  796  c+Ui 

0. 

-1*64092-08 

-2.87ML-06 

.6928 

i • 3563c ♦ J  t 

5  •  6(3  9o  E+Gl 

u . 

-9. 82032-10 

-2.1039E-07 

.  7  1+5 

1. ob53E ♦ J  1 

5 «  0  0  34  c+  0 1 

u . 

-1,46792-11 

-3. 9164E-09 

.7361 

2 . 1 377  c  +  u  1 

4 . +6U9E+  0 1 

3  . 

0. 

0. 

.7578 

2. 778»c+ Jl 

3  •  98£9  t+  0 1 

o  • 

0. 

0. 

.7794 

3.71882+01 

3.5723  2+01 

3  • 

0. 

0. 

.8011 

5.1bLlc  +  Jl 

3.2  354  c+  ol 

u  • 

3* 

C. 

.8227 

7. 4 92  IE ♦ Ji 

2. 98olt+  u  1 

u . 

0. 

0. 

.8444 

1.  15382+  J2 

2. 65592+31 

G  . 

0. 

0. 

•  8660 

p  y  ,  J 


r  i./iuT- 1 


o YC luT-xun  FREQ 


(lOCAL) 


.4000 


OuA.  3  /  x  2  f  R  L  4 

ryF  RAC  T  1  Vc. 

LUCAL 

AMPLITUDE  exponent 

SINE 

i  4  0  L  A  *  *  2 
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•  5413 

* .  3o 12  2  *  J  j 

9«  01732+ui 

0. 

-7  •  3p9*tt-  Go 

-3. b74GE-u4 

.5629 

M.  Jao  9c.4  J  J 

t .  lb  3fc  l+  0 1 

8. 

-b.0b93t-Gb 

-2.8287c-04 

.5846 

? •  7  30  62  +  J  J 

; .  37itt2+gl 

0. 

-3. 2636E-C6 

-2. 12982-0“ 

.6062 

o.o4C9t+o) 

6.63932+01 

J  • 

-1. 8b39t-0b 

-1. 4068 t-04 

.6279 

7.  7  9c  jL*-  j  j 

j  •  96m2  2+  41 

u  • 

-8.7o7lt.-G7 

-7 . 84872-05 

.o495 

it  i  4u  3C.  «■  j  .1 

v  •  3 h  3  c.+  J 1 

u  • 

-  2 • ott52E-0  7 

-2.8733t.-05 

.6712 

1.  llu9E*  'll 

4  .  7  7  3b  2+  iil 

U  • 

-o. 1797E-G8 

-b. 6993c-0  6 

.6928 

l»3p63t+ol 

4 ,  2 b  3b  c.+  j  1 

C*  • 

->t.*tbH9t-G9 

-7 . o634t-G7 

•  7 14b 

1 . 66532  +  J 1 

3.  7  9  99  t+ ox 

4  • 

- 1 , lb/ 2t-i0 

-2. P0b52-G8 

.7361 

2.  1  37  7t  <■  j  i 

3*  3oa2c+ux 

n 

U  • 

J. 

t  • 

.7578 

2.  7 7 ti 3l  f  81 

3.  0  2  2d  2+  u  x 

u  • 

J. 

ll  • 

.7794 

3.  7  1862  <■  J  l 

c«7113u^al 

C  « 

u  • 

l . 

.8011 

3,lCtlc.4Jl 

2.  *6  ol t+  o  1 

6  • 

j  • 

G . 

.8227 

7 . 492xc  j  1 

2.  2o7y  c.+  o  x 

V  • 

j  • 

c  • 

•  84*4 

1.  1638c.  *■  j2 

2  .  16  63  l+  0  x 

J  • 

J . 

4. 

.8660 

100 


P'<OP4*.  F  J./cUTR-.  CYCLU,rRoN  FRlG.  (LOCAL)  =  •  & 0 0 0 


l-U^Al  3/iOfKL  B 

XcFRAC  TiVc 

LOCAL 

amplitude  exponent 

SINE 

IND£X**2 

S  TRcAMi N& 

QUIlSCeNT 

NELT 

(LATITUDE) 

(  B) 

(  X  N  2  ) 

l  X  J  3  ) 

(XJQ) 

(XNIXY) 

(X) 

1 .  9  )  u 

4  .  1  2 12  ,c.+  0  2 

*J  • 

fiL*-Git20E-J&4_ 

1.39.S8fc-M 

_  ji.9  615 

1  •  m  u2 1  c.  *■  1  j 

4.  0wi?3£  +  02 

U  . 

2.77U6E-C4 

1. 3b94t-93 

.w2l7 

1  •  J  0  8  5  t  +  J  J 

? . 9  b  j j  c+  c2 

J  . 

2.&531E-C4 

1.  <L?78t-0  5 

•  0  43.3 

It  J  19  2 1  f  j  J 

3  .  8  c+  <j  2 

U  • 

2. 2O35E-04 

1. 126C L.-0  3 

•  o  650 

1  •  J  343E  +  j  i 

3  .  7  ti  3 1  £+  j  2 

a . 

1.7430E-04 

9. 1658E-04 

•  o  86b 

1.0541E  +  )i 

3  .  bo  22  £+  J  2 

0  • 

1.3359E-04 

7.2839E-04 

.1083 

i.  J787L+ J j 

3. 52  Jgt+w2 

V  . 

9. 3844L-05 

6.3438E-04 

.1299 

1.  il'85t  +  „  J 

3  <  3  b-+3  £+  j  2 

4  »  u  7 i  to  £  - 1 4 

6.7 823E-05 

3.4619L-04 

.1516 

1.1h39£+  J  J 

3. 1977E+02 

3.93i3c-05 

1 . 5623E-0  5 

3A4  92  7k-Q  4 

.1732 

1.1854E+UJ 

3  .  C  242  £+  oZ 

3  •48“fOL."*ii3 

-9.947  0E-06 

2.  3429t-Q2 

.1949 

1.2335E+JJ 

2. 8472£+i;2 

”2 . 1  j  7  d  c  —  u  3 

-2.9715E-05 

-1.641 It-  9  2 

.2165 

1.2389E+ J j 

2 . t  o97 £  +  u  2 

-1 . o 1 99  t-fc  3 

-4. 699GE-Q5 

-1.  2919E-02 

.2382 

1.  3»2hL«-  J  J 

2  .  4  942  c+  u2 

-2.h29h£-£i7 

-4.7 192E-05 

-3.9740L-04 

.2598 

1*4  25  1 L  +  i  J 

2. 322b£+C2 

0  . 

-5. 2or»2t-£l5 

-4.7312E-&4 

.2815 

1.5 )8lt + J  j 

2. 1363  t+  0  2 

0  • 

-  5 • 5U85E-05 

-6. 4544L-04 

.3031 

1.0.29E  + J  j 

1  •  99oo  £+  u  2 

c . 

-6. 4653E-Q5 

-6. 916b£-04 

.3248 

1.7  11 3c.  +  9  J 

1.  S4h0  c+  0 2 

u  • 

-5.UO86E-05 

-b. U498E-04 

.  3  4b4 

1.63546+9 j 

1 . b992  £+42 

0  . 

-4.7  O47E-05 

-b. 24l6t-04 

.3661 

1.9777c+ J J 

1.  6  b  2  h  2  +  02 

G  . 

-4. 284OL-0 5 

-b. 20b7t-&4 

.3897 

2.1414c ► j  j 

1  •  433oL+u2 

u  •> 

-J.6771E-05 

-5. 91h4£-U4 

.4114 

2  •  3  30  5c  +  j  u 

1. 3127L+32 

u . 

-3. 1397E-05 

-6,b281t-(i4 

*4330 

2.5h9oc+Jj 

1 .  1997  c  +  12 

J  • 

-2.7 145E-0  5 

-5. 446UE-04 

.4547 

2.  8u6hL  +  !jj 

1.  u  9+1  c+  L  2 

J  « 

-2.247LE-05 

-6.  ti662c-04 

.4763 

3.  iu5ut+ j  ) 

9 . 9  p  92  £+  0 1 

0  . 

-1. 8127E-05 

-4. ol3U£-U4 

.4960 

3.  456 -+c  +  J  ) 

9  •  o  ho 7  £+01 

u . 

- 1 . 4970E-05 

-4. 3219t-04 

•  5  196 

3.  87tt3t + J  ) 

b«  2 3  t] c £  +•  4 1 

0  . 

-1. 1796E-05 

-3.  8547t-l»4 

.5413 

4.  3612-.  +  J  l 

7  .  h  1 7  7  t+  u  i 

u  * 

-6.7  684E-06 

-3. 317bE-u4 

,5629 

+ .  J88  9c+Jj 

c  .  1 9  +  7  t+  0 1 

L>  • 

“ o . 2o99l-66 

-2.7  u67E-G4 

.5846 

5  .  7  6  L  5  C  t-JJ 

=  .  u  2  37 t+ul 

L»  t 

-4. 1881E-06 

-2.  126&C.-U4 

.6062 

o.  o  4  5  9  E  +  j  j 

3. 41  jCc+ui 

0. 

-2. 4770E-06 

-1.4745L-84 

.6279 

7.7yCuL*')j 

4  .  6s57  c.+  Ul 

U  • 

- 1 • 2b96E-C6 

-8.9462L-05 

.6495 

9. 2hL  5c  +  3  j 

+  •  3  4  3->  £+  y  1 

U  • 

-5 • 1o53£-07 

-4. 3625E-05 

.6712 

1 . 1  It  9£  +  j  1 

o.  o77  7  c+^l 

u . 

-i. 2137E-07 

-1. 2430E-86 

.6928 

1. 35b3t  + J1 

3. 4  5  04 t+ui 

u  . 

-1. 28o6E-08 

-1. 6106E-06 

.7145 

1. oo5  3£  +  J  1 

3,  u  78L  £  +  oi 

Li 

- 4»  7997E-10 

-7 • 4795E-Q8 

.  7  3b  1 

2 . 1  37  7  E  +  Ji 

2  .  7  4  lo  £+  J  l 

U  • 

-2.7  229E-1 2 

-6. 406UE-10 

.7578 

2.  7785C.+  H 

2. 44oBt_+oi 

b  « 

9. 

a. 

-7794 

3.  7 18dc  +  Ji 

2  .  1947  c+  j1 

u  • 

J  • 

0. 

.8011 

i.  lbu  lc  +  J  1 

1. 9  6 6o  L+  1. 1 

i 

^  • 

J. 

0. 

.8227 

7.h921£+j1 

1. 83b3c+01 

0  » 

u  • 

it  • 

•  8444 

1.  l>3o£  <■  )Z 

1 .  7  po4  £  +  0 i 

u  • 

J  • 

l. 

.8660 

1 


101. 


P<0P4b.  F  Uu./tUT-L  .  CYCLuT^Uim  FRtQ.  (LOCAL)  =  .6000 


-Ol-Al  J/.QWl 

J  KLFRAGTli/t 
IUDEX*»2 

LOCAL 
S  TKlAMi NO 

AMPLITUDE.  LXPUNtNT 
QUit.SCE.NT  NlT 

SINE 

(LATITUDE) 

(  B ) 

(XN2) 

(XJS) 

(XJQ) 

(XNIXY) 

(X) 

i .  0  J  L  t  ►  J  J 

4 ,  idOLOd.j 

u  • 

9.6B93E-04 

3.  1697E-0  3 

.0000 

l.uG2it*)J 

t. 1679l+  d2 

u  . 

9.  *+3l7t-04 

3. Ill 3L-0  3 

.0  217 

1*  u  j  6  p  c.  ♦  ^  j 

't  •  1  >J  7>  u  C 

J  • 

8.  OOiiit-04 

2. 95C4i--0  3 

.0433 

1  •  d  t92  L  ♦  J  j 

■*  •  C 1  <j  1 tj  2 

u  t 

7  •  9ut»Bt-04 

2.7 165L-0  3 

•  0o50 

in  J40i.tJJ 

■U  .  0  O  J  2  C.+  L  «. 

j  • 

o.  *+759c-04 

2.  3068E-0  3 

.0  866 

1«oB41l+Jj 

3.  7  2  3d  E*'J2 

d  • 

-/.  9607E-CH 

1. 847  3t-0  3 

.1083 

l.c76/E+jj 

3.  B*+o7  l+  32 

u  • 

3  •  0o7 1E.-Q4 

1.S180L-03 

.1299 

l  •  i  i.  8  "j  l  +  J  .1 

3. 3ddit.  *U2 

d  • 

2 . 0u9lE-04 

l.l7i2t-03 

.  l5lo 

l.l>’»3,dLFj 

0.  1543  c.*  j2 

id  *3 

1.023 7 t- 04 

o. 1539L-04 

.1732 

1.  1  d 5tE  1  j 

2. u+  ut; 

7 .6972i_-( 9 

1. Uo38t-04 

5.  066lt.-04 

.1949 

1.2  >  3 ‘j  L  *■  J  J 

2.7h7o(_+i.2 

u«16^6c.”u3 

4. 91o2c-05 

1. 1600E-G2 

.21b5 

1.2>69c.«-j.> 

;  .  5  t  ’id  c.+  u2 

3. 2323c.  - c  3 

1 .  B4J9t-Ci6 

1. 0247ii-02 

.2382 

1  •  >5  l>  c  4 t,  *  ;  J 

d  •  3  53.1  C.+  u2 

-2  .  d  7r>Bc.  -v  3 

-  2.  ‘♦'+63L-05 

-1, 7tib6t-02 

.2598 

1.  LB  1  l  *■  j  j 

2  •  1  /  j  7  _+  1 2 

-I.30d2c.-u6 

68992.-05 

-3.2739L-04 

.2815 

1  •  d ifclc+Jj 

1 « 99vo  l>  u2 

V  • 

-4.  <.o8lfc-C 5 

-3.2753c.-04 

.3031 

li3.  C  Jt*  J. 

i.  63.J3C.+  L2 

u . 

-j • 0341E-0  5 

“4. 1291L-04 

.3248 

1*7  1 1 3  t.  ♦  j  i 

1  •  o7  59c.+  32 

o  • 

-o. 2885E-C5 

-*».  0u42E-L 4 

•  3464 

1.  )  jIpiL  >  j  i 

1  •  :>  J 1 7  t>  J  2 

U  • 

-»• U464E-05 

-t.  o53liE-04 

•  3681 

1  *  9777  c.  j  j 

3  «  0  9  7  o  ii  +  j  2 

Li  « 

-4.  toaiE-u  5 

-5.  l2b9c.-04 

•  3  8^  7 

2»l*.ltC.^uJ 

1. 2  7  33  t+ jc 

U  • 

“4 • w  756L-05 

-5.1254E-04 

.4114 

2.33i/>i_*,j.) 

1 .  1b83  c.+  c2 

t . 

•  3. 40O4E-O  5 

-4.  ob/bc-04 

•  433Q 

2.j4  9at<-)J 

1.  c822E*o2 

u  • 

-2. 9435E-05 

-4.b747c-04 

•  4547 

Jit  :5>tc  +  )  j 

9.  b4-»2c!+u1 

0  . 

-2.  *f243c-05 

-4.  35u9c-C4 

•  4763 

i.loBot+jJ 

b  •  0  44u  E+  ill 

If  • 

-2 , 0242E-0  5 

-4.  121 5t-04 

.4980 

3  •  4  t>8  4  L  J  J 

7.  81  64  u*  L  1 

i  t 

-1.6530L-05 

-3*  635 9L-0h 

•  5196 

)i  )  /b it  j  . 

7  .  L5o4t,+  dl 

J  • 

-l.iu55E-05 

-3. 4b95E-u4 

.5413 

4.  3ol2t  F  J  ), 

3  .  3  39  3  t+  w 1 

w  • 

-9.9653E-06 

-3.  0499E-34 

•  5629 

4.938  it  Jj 

5 .  7  2  c.+  J 1 

u  • 

-7.15O2E-06 

-2,  537  3E-G4 

•  5  846 

J  «  7  30  tt  +■  j  j 

5. lObBE+m 

c  • 

~4.  83o3E-06 

2  «  0  u  L  7 1  -  0  4 

.6062 

7  %  7 

4 , 6  j  3*3  l.4 u  1 

U  • 

-3. 14o9E-06 

-1.  53ulc.-t4 

•  6279 

4.  1 1  32  L+  d  1 

J  • 

-1. 7028E-06 

-9. 8233E-05 

•  6495 

9.  2  4L.dE*  u  - 

3.  b772c.+  dl 

J  « 

-7 , ?l*7t-07 

-5.1942E-05 

.6712 

1  •  1 i l  it  +  /  l 

3 . 2  7  3  l>  u  1 

t  • 

-2.  12  Jut-07 

-1* 7763E-C5 

•  6928 

1.  idb 3E  *  .1 

2 . 91  do  c.+  oi 

w  • 

-2.  o37b£-08 

-2. 9211E-06 

.7145 

1.o8b3E*j1 

c.d97ot*Li 

0. 

- 1. 4L18E-09 

-I.6U14E-07 

.7361 

2. 1  37/  t.  *  j  1 

2.  dll'ic+ol 

u  • 

-1. BOJlt-ll 

-2.5829E-09 

.7578 

2.778  •;  t  ^  / 1 

2. u  6  29  2+  u 1 

U  • 

0. 

t  • 

•  7794 

3 .  7  i  t>  5  u  *  ■  i 

1  «  O  B  J-r  2*  j  1 

w  • 

J. 

U  • 

.8011 

P  t  4  t'  •/  1  Ci  ^  /  i 

7 i litit*  ji 

ltO/aiuti.. 

w  « 

3. 

u  • 

•  6227 

i«  '■>  4  3  9  t.+  u  ± 

•  • 

J  • 

L  • 

•  8444 

l  •  1  C  3  3  L  ♦  .  ? 

1  •  “»  tr  1  U  W  ♦  w  1 
«  *  .... 

L  t 

‘j  ....  .  • 

l  % 

.  86oU 

102 


P<OPAG.  F  <■  J./cuTKL.  CYCLCHkON  FKLQ.  (LOCAL)  a”  ,7iiOC 


uUlAl  0/;'QrKL 
(B) 

1  •  J  J(j  v  t.F  o  J 
1. Ju21L+ j J 
1.0o85L4')'i 
1.  0192E  +  J  j 
1.  J  34 3t  *■  j  J 
1  •  0941c  *■  J  j 
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1.6c29t+- J  J 
1.7il3L<- J  J 

1.  6  39*L  *■  J  J 
1.9777c  *  J * 

2.  l4ltE*  Jj 
2.  3309c.  *■  J  j 
2. 949ot  *  J  ) 

2.  b  o94t  *■  J  j 
3.1u5dc*Jj 
3  •  4  s  8  4  c  *  0  j 

3.  6?83c  *•  J  J 
4. 3oi2c  «•  J  j 

9669c  +•  j  j 
5.73(j9E4jj 
o  •  b459c  <■  j  w 
7.790wt^Jj 
9.2hC5E*  Jj 
1. llo9t*  Ji 
1  .  3bo3t  *■  J  1 
it  o693c  * Jl 
2.  1377L  f Ji 
2. 7769c  4  J  1 
3. 7186c  ►  Ji 

■3,  iDOlL  *■  Jl 

7. 4921E  *01 
1#  1936c  *•  j2 


AcFRACTI  Ml 
INDEX**2 
( A  N  2  ) 

4 1  7.853l+  02 
4.  7536L+Q2 
4#  Lb J7  c+  02 
*♦#  5l29c+U2 
4.  31 96  E>  j2 
4.  0 9 23 1+  u2 

3, 6419E+02 
3. 5793 t+02 
3#  3127  c+ J2 
3#  C9J1E  + 02 
2 #  7 9o3 E+  02 
2#  9991L+u2 
2. 328^E*u2 
2. 11 79  l+  u2 
1. 9225t* j2 
1#  7433L  +  02 
1. 579it.+  G2 
It  ‘♦291c*- 02 

1 .  292**  c+02 
1. 16/9c4G2 
It  0  9 -*6  c+  j2 
9 . 9 1q3  c  +  u 1 
6. 5799E+G1 
7  t  7279E+01 
bt  9938L*-ui 
3t  29J2E+U1 
9t  61U9L+U1 
9  t  0  3  J2  c+  0 1 
4t  903lc+0l 
4  •  u2b3c+lil 
3  •  9927  t+u 1 
3 1  2  0  2u  L+  0  1 
2. 8902  c+  01 
2t  3347  E+ui 

2.  2p39l+0x 
2. o091c+ul 
It  7867  c+Ci 
It  6  U  44  C>  u  1 
1 . 4d  *3  c+  u 1 
It  3*+3o  C+-  0  1 
i. 2  o  p  2  c  +  0  x 


LOCAL 

amplitude  exponlnt 

SINE 

STREAMING 

QUItSCENT 

NET 

(LATITUDE) 

UJS) 

(XJQ) 

(XN1XY) 

(X) 

u  • 

_3j  228%E-li3_ 

b, 8124E-U  3 

.0000 

U  t 

3.1553E-03 

b«  7 050E-0 3 

.0217 

U  t 

2-t.9192.E-0  3 

6. 333EE-JJ 

tJ)43.3 

u  • 

2.9593E-03 

9t  7447E-0  3 

•  Qb50 

V  • 

2 1 1839E-0  3 

5, 1339E-C  3 

.0  866 

t . 

1. 7386E-G3 

4.3284E-03 

.1083 

o  • 

.  i*32§7E-(J3 

3.5295E-03 

-  .A.239 

0  t 

9.6282E-04 

2.7648E-03 

.  1516 

c  • 

6.  4341E-04 

2.007  3 1- 03 

.1732 

u  • 

4.1G44E-04 

1*  3996E-03 

.1949 

&  t  7 J  J  9  c  —  1 0 

2 1  5294E-04 

9. 4770E-04 

•  2165 

2.7  2jUc-h3 

1# 2432L-04 

1.1760L-02 

.2382 

1 .  4  4 1 5  c.  - L 3 

4t.72t3E“Q5 

_b»  8167E-03 

.2598 

-2  is l79t-G  3 

- 3t bupbt-06 

-1.2626E-02 

.2815 

-It2up3t-u8 

■ 2t  9129E-0 9 

-l.b513L-04 

.3031 

0  . 

-3.6136E-Q5 

-2.4132E-04 

•  3246 

U  • 

-4.6372E-Q5 

-3.2846E-Q4 

.  3  4b4 

C.  9 

-3. -271E-05 

-3.9949E-04 

.  3661 

d  t 

-4. O283E-09 

-4.3146E-04 

.3897 

J  • 

“4#  3930E-05 

-4.3848E-04 

.4114 

j  • 

-3.  6360E-.05 

-4t.139.tt-.U4  . 

_ ait  .310 

c . 

“ 3 1 2U69E-05 

-4.1373E-04 

.4547 

M  • 

-2.7  J45E-05 

-3.9o67E-04 

.4763 

0  . 

-2.2G88E-Q5 

-3.6999E-04 

.4980 

U  t 

~1#  773*tE-05 

-3.3989t-04 

.5196 

u  • 

-1.4996E-05 

-3.2176E-U4 

.5413 

d  9 

-It  130 lt-05 

-2. 8801E-04 

.5629 

t  • 

-6t  2272E-06 

-2.4377E-04 

.5846 

0  • 

-5 • pp17  E-06 

-1.9248C-04 

.6062 

U  8 

-3. 6189E-06 

-1.4789E-0H 

•  b279 

o . 

-2,U»74t-Q6 

-9._9997t-0  5 

.6495 

o . 

-9.  9099E-07 

-9. 7793E-U9 

.6712 

L  • 

-3.  2162E-07 

-2.L770E-05 

.6928 

a  • 

-9t  Z'tJbE-OS 

-4.9b38t-06 

.7145 

O  • 

-3. 2872L-Q9 

-3.5769E-07 

.7  361 

o . 

-9. 1H24E-11 

-7.1208E-09 

.7578 

0  • 

0. 

0. 

.7  794 

0  • 

J. 

0. 

.8011 

u  • 

0 1 

t  • 

.8227 

0  . 

Jt 

t  • 

.  8444 

O  • 

U  • 

bt 

.8660 

rj  iv 


103 


P  <0+  A  b «  (■  L.  oYbLUTnibN  F  Kc  U  •  (LOCAL)  =  .6000 


CK  A  L  J/.JFkL  'J 

R..FR4CT1  Vw 

LOCAL 

AMPlXTODL  tXPONtNT 

bINc 

lr-IDtX**2 

CTkcAMINC 

QUlcSCtNT 

NET 

(LAriTUOE) 

IR) 

(  A  N  2  ) 

(XJb) 

(XJQ) 

( XNIXY) 

(X) 

1.  J  jJ 

g  4  27  7  o  c  + w  2 

0. 

6 • 0  664E  — 0  3 

l.o 9446 -02 

•  U  Go  J 

1  .  J  i,  2  1  L  +  )  1 

o.  2l-»2  t*  w2 

u  • 

6. 73096-03 

1.  L 866^  ~m  2 

.0217 

1.  j  ilnt  *  j  j 

04  L  3 13  l+  L  2 

W  • 

6. 2s55E-G  3 

1.  L597t~t-  2 

.  0433 

i.  L  19  ^  L  *  ,  J 

’j  %  7  *+  79  c.+  b2 

U  . 

7.p279E- 03 

1.  0157t-G  2 

.0650 

1  .  J  j>tOL  +  J  1 

t?4  39uoL+j2 

J  • 

0 *63796 "03 

9.42962-03 

.  0  866 

I.JjhIc  +  jJ 

t. 93322+02 

wf  • 

6 4  42586-03 

6.4847  E-C 3 

»  1 0  8  * 

1 .  j  t  6  i  c  +  J  j 

44  job',  c+  J 2 

U  • 

4*  21306-03 

7.  2259c-G3 

.1299 

ii  Kb  jl*  J) 

44  1 4  5  i  c+  Cl  2 

J  • 

3.u625E-03 

5. 85156-03 

.1516 

1»1h39L+Jj 

3i7itulc4ijt 

-  • 

2.12146*03 

4.48812-03 

.1732 

1  t  1  ^  5  t  k  *■  J  J 

3  4  36  j  3  L+  J 2 

J  • 

1.37476-03 

3.25792-03 

.1949 

l«2i3;jc*-'Jj 

2  •  uliibE+iiL 

J  • 

8. a4l7E-G4 

2.27602-03 

.  2 165 

iiKoJt*  j; 

2  4  d  92-»  t+  j2 

2«2f*ufu*t 

4.^2376-04 

2. 04196-0  3 

.2382 

1.3u24c+j- 

2  •  *t  u  69  c.+  u  c 

4  *H  3 

2. 39766-04 

1.56286-02 

.2598 

x.4c51c+j  i 

2. xn9l l+02 

i .  j4ob£-04 

-1. &3112-0  2 

.2815 

I.;ubltKj 

i. yiyoc+oC 

-1 #3  2  d  f c*i  3 

1. 73466-05 

-5.D9702-U3 

.3031 

JJ 

1. 7  !.j3  c+ut 

L  • 

- 1. 74306-05 

-o. 58122-05 

•  3  248 

1.  7113c.  +  j  J 

1 4  b  3  3b  c  +  l)  2 

U  • 

-3.322bt-05 

-1»  b5?3L.-04 

.3464 

1  •  3  3  5  •+  l  +  j  : 

1  4  3  7  1 7  c+  0  c 

t  . 

-4. 2bJ9t-G5 

-2. 70102-04 

.3631 

1  •  9  t  if  1  L  +  J  J 

1  «  c 2  7-*  L+  U2 

^  • 

-  4  •  7i<J  26-05 

-3.  3933c-ii4 

.3  697 

K  KltL>  j  i 

1 4  293/6+J2 

u  • 

-  4. 4oob£-05 

-3.66176-04 

•  4  114 

3  •  o  30  j  c  +  J  j 

9 4  8379 t+0i 

J  • 

-4. 09756-05 

-3.62856-04 

.4330 

2 .  :j  4  9  3 l  +  jj 

a. bi Jit+ui 

0. 

-3.43566-05 

-3. obb7E-04 

•  4547 

2. 

7  4  6  8  96  l+  u  1 

u » 

-2.69436-05 

-3. 53442-04 

.4763 

3«j.u6jL  +  Jj 

7  4  \jo*n  c  +  u  1 

J  • 

-  2 • 3y37£- 0  5 

-3.  296bc-04 

.4980 

3.  ‘♦3(3  HC.  +  J  J 

6  4  3  2  27  6+  J1 

c  • 

- 1. 9b306- 05 

-3. Io76t-04 

.5196 

3.6763L+jJ 

5  •  o  :>o3  6+  J 1 

J  • 

-1.56236-05 

-2. 91372-04 

.5413 

t  .  3 ol2  c  +  j  J 

K j5o/ c+  *  x 

J  • 

- 1 • 24456-05 

-2.6943L-04 

.5629 

4 .  9  o  6  9  c  +  J  - 

h  4  5 169  c+  u 1 

a  • 

-  6 . 95o4£-0  6 

-2. 2627  2-0  4 

•  564o 

3(7K3LKJ 

•*  4  i.  3  J9  l+  u  1 

v  • 

- o. 36ol£-06 

-1.69382-04 

.6062 

b .  o  h  5  9  l  +  j  . 

2  •  5  9  34  t+  x  1 

j  • 

-4. Io29£-06 

-1. 4592c-L4 

.  b  279 

7  .  7  J  1,  Jtbj  J 

3.  1999L+ux 

J  • 

-  2 . 34416-06 

-1. u 62 5 2-04 

.  b  4 1  -j 

i  .  C  n  1'5  L  ’>  J  J 

2  4  64646+  J  1 

li  • 

-1.27776-06 

-b. 417  46-05 

.6712 

i.lKjtt.l 

64  6  2  93  c+  g  x 

r 

w  • 

-4. 4490 6- C7 

-2.  717  O2-L  p 

.6923 

1. 5  bb  3  c  +  J 1 

2  4  24  7  y  6+  u  i 

U  • 

-i.02o 26-07 

-7. 71822-0  b 

•  7  145 

liOO^OL  ^  j  k 

i.99o.,cKl 

r. 

w  • 

-0. 59346-09 

-b. 2J082-07 

.7361 

2.l3/7u+jk 

1 # 7 1 i 1l+ux 

U  • 

-1. 34J16-10 

-1.60496-08 

.7578 

2.7  /OJtKl 

X  4  ’J  O  3  i  L  J  J  k 

6 

)  • 

J  t 

0. 

.  7794 

3 1  7 lo ol  +  J 1 

1  .  42  Jkl  t+  Jl 

it  • 

Vi  • 

U  • 

.8011 

5. logic  ►  Ji 

1  .  2  6  7  3  c+  u  1 

G  . 

0. 

l  • 

.  6227 

7.  4  92it  Jl 

1  4  1  O'Jb  t+  ■.  1 

L  • 

u . 

it  • 

•  8  444 

li  l53at.  KJ 

i  .  1  3  3  .  t+  u  1 

*4  « 

J  • 

O* 

•  8660 

104 


h=1.j  N  =  l»S’  ELECTRON  TtST 


0  P  A  0 «  F<aj./tQTkL  CYClGTKJN  FREQ 

■  1  j  [t 

•  2  o  U 

•  3  W  U 

•  >♦  U  l) 

.  5  L  0 

iuui] 

•  /  'J  J 

•  o  U  U 

•  9  u  o 


NtT  AMPLIFICATION  (DECIBELS) 
i. 736663961420t+ul 
1.850691 97 obbyt +01 
-4.+423b8976719t+0l 
- 1 .  2 64 7  77 7 3 1  <*2  9 E  ♦  Ud 
-3.60  67 9124972b t+0  0 
“5«  26  98  5121 5.5  6  9t +  01 
-7.692933324534E+C1 
“1.34 542245 7571L+Q2 
-2.4387448b8t>75E  +  02 


PLmSMA  D  i  N j IT  Y 
UUI£SCt.Nr/OOLD 

3.6bijjw“C3 


PITCH  ANGLE 
QUIESClNT  norm 


PLASMA  DENSITY 
STREAMING /COLO 

1.3612E-15 


PITCH  ANGLE 
SIR t AWING  NORM. 

5.2957t-01 


% 


<  -r  *1 


I 


1%. 


PROGRAM  CYCLES  (INPUT, 
UIHf  NSION  A(100) 

DIMENSION  HQ { 5  0) 

iALpHftS(50) ,EHCQ(50) ,  EMC 
2SPL  CNQ ( 4 ,5C ) 

COMMON  HQ,HS , XNQ.XUS ,Al 

1  SPLCNS.RAO.PI ,TH0,XN0 

2  ^.XNZ.B.XN.OFG.OETASl 

3  NTYPt  G,NTYPE,(10ETAPtn 


PROGRAM  FOR  HAROLD  LIEHOHN 
PROGRAM  MODIFIEO  2/28/72 


OJT?UT , TAPEI.INPUT ,TAP£3*0UTPUT ) 

.R<  100)  ,XSAV£(100) 

’??i?C',«XN‘J(50)  »XKS<50)  »ALPHAQ<50)  t 
i  (  50)  ,S  *LCNS(50)  ,SPLCHG(4,50)  .SPLCHSC4 ,50)  , 


SPLCNQ, 

MCO  2  r 


ACT  1445  TP41233 

ROBERTA  KERR 


WHISTLER  mode  amplificati on 


npa  c r  =  1 

RAD=  .  fl 74532925 
PI  =  3. 141592654 
0EGs57, 2957795 


>2000) MONTH, NUAY ,NYR 


OF  ENERGY  OISTRIB( XNQ,  XNS) 


REAP  <1 
TI TLE  CARO 

5  READ  (1,1030) 

IF(FPF,1)  5000,6 

6  CALL  SPRINT 
WRITF  ( 3,1030 ) 

INITIAL  VALUES  ANO  CONTROL  PARAMETERS 

READ  (1,10  40)  THO.XNO, NT YPEQ,YS,Y£tDELV , M ALPHA  .NX  ndctio  Dnoc 
WRITf  (3,1050  )  THO  *  PORE,  NA  L  PHA ,  XNO, NT VPEQ, NX, DEL Y.NBET AP * YS »VE 

READ  “ITCH  ANGLE  OISTRI B(HQ, HS)  ANO  LOGIO 
READ  (1 ,1000) NUMHQ,NUKNQ,NUHHS, NUNNS 
IF(NUMHQ.LE.  0)  GO  TO  20 

r^r  co’!S10)  (ALPHAQ<I>  HQU),I*l,NUHHO) 

CALL  SPLICN ( ALPHA  Q , HQ ,NUMHQ , SPLCHQ) 

I F  (  NUMNQ,  LE  •  rj )  GO  TO  21 

?f.\r  tl  ’J?10)  f-MC;j(I)  ,XNQ(1)  ,I  =  1,NUNNQ) 

CALL  SPLI^N  (tMCQ,  XNO  , NJNNQ ,  SPL CNO) 

IF ( NUMHS.LE.  9 )  GO  10  22 

PAM  ‘i:*1;10’  »HS  ( I)  ,I  =  1,NUNHS) 

CALL  SPLICN  (  ALPHA  S,HS,NUMHS  ,  SPLCvHS ) 

IF( MUMNS.LE. 0 )  GO  TO  23 

PA|A|°  ^di’J?10>  <u',CS<I}  tXNS(I)  ,1*1  .NUNNS) 

HprTF{3a020?MCS’XNS,NUMNS’SPLCNS> 

MMA  X  =  MAX0 (NUMHQ , NUMHS , NUMNQ, NUMNS) 

00  ?5  1*1, MMA X 

WRITE (3, 1001 ) 

WRITE (3,  10  02 1 
WRITE (3, 1003) 

WRITE ( 3, 1004) 


20 


21 


22 


2  3 


2  5 


IE(J  LE.NUMMQ) 
IF ( I .LE. NUMNQ) 
IFC  T  .IE.NUMHS) 
IFII.LL. NUMNS) 
CONTT^Ut 


ALPHAQ(I) , MG  < I) 
EMCQ(I)  *  XNG Hi 
ALPHAS!!)  ,HS(I) 
EMCS(I) ,XNS(I) 


F  I  £  L  n  L  1  f‘f  PARAMETERS 


TH i P  =  Th3*  RAO 
COSTn  =  COS(TMOR) 

SINTH  =  SIN ( THOR ) 

HP1  =  S.65C*J4*  SQKT(XNQ) 
wC',  =  r.SE  +  Q6#oQSTH**6 
wc 

1  P ( NUMNS .LE • 3  »  00  TO  26 

xt  l  =  S'JRT  ( i .  - 1 . /( ( io.  *•  £NCS<  1) /PORt*l. )  •♦2)) 

•r  r  f  *-■?=  jrJRI  (l.-i./(  (lO.**,EHOS(  NUHNS)  /P0RE*1. )  **2) ) 
NURHA'_!.>  PITCH  ANi.Lfc  UiSTRi  3  JT  1  ONS 
2  f  Hi  »  f  =1 

CHI  SJ  MP<  0.  ,  Pi/2.  .NALPHA  ,C0) 
co=ro*A  .*PI 
CQ=  1  ./oCi 

I'(  M'HmS.LE. 0)  GO  TO  160 

NTvfr  =  2 

CALI  SI  MPT  ALPHAS ( 1) *«AO,  ALPHA S( NUHHS) • RA C ♦ NAL PH A ♦ CS ) 
rs=r< *l.‘Pi 

f  S  -  *  /  c  s 

SLT  P’ffAbATION  FW£CUENCY 

16  0  NUMVtr 

Y  r  Y  <* 

computt  c v ci ut run  resonance  amplification  at  y 

17  0  Y2=  v  *  y 

CALL  SPRINT 
NRTTf  CJ.UIO)  V 
WRTTF  (  3 , 1 1  0 v> ) 

♦10  CALL  S1MPXT0. , SI  NTH , NX,  FI ) 

NUM  Y  =  NUMY  ♦  t 

ATNLfMY)  =-2.  0  *0 6* Y*  GOS  TH* • A* F  I 

XSAVF (NUMY) H 

Y  =  y  +  Ot  L  Y 

IF  » Y-YE) 170 , 170, loO 

EVALUATE  MINIMUM  RESONANCE  VE.OCITY  COMPONENT 
IttO  YArY** 

IF (PPFE.GT. 13. )  GO  TO  362 
Z  A  =  Y  A 
GO  TO  363 
362  ZA=1 

J63  XN?=1.+WP0**2/(HC02*2A*(1.-YA) ) 

XN=  SORT ( XN2 ) 

YA2=Y A*YA 

Cl=  1  /(XN2*YA2) 

GET  A  R=  1 .  -  SQRT  (i.-(l,  *Ql)  M  1.  -1./YA2)  I 
GFT  A  c  s-3ETAR/(XN*U.*Cl>  ) 

NORMALIZE  HOT  PlASMA  DISTRIBUTIONS 

NTYF  -  =  5 

IF  <nL TA%-. 001) 500,502,501 


hk;. 


iJO  LALl  SIMP(3ETAR,. JOi,NBETAP,£PSLQl) 

>02  STRTs.001 
GO  TO  503 
501  STRT  =  BE  TAR 
EPSLOl  =  0.  0 

50  3  CALI  SIHP(STRT,  .  999,  NBE  T  A  P,  EPSL  Q2) 

£PSLNQs£PSLQl ♦EPSLQ2 
IFCNUMHS.GT.  3 )  GO  TO  50  4 
EPSINS=0. 

CS  =  C. 

GO  TO  505 
5Q4  NTYP<-  =  6 

CALL  SIMP(BETAS1,9ETAS2 .MBETAP .EPSLNS) 

»0  5  CALL  SPRINT 
NRI  Tf  (  3 ,10  3C  > 

MRITr (3,1060) 

WRIT* (3,1070)  <XSAVE(1),A(I)  ,I*1,NUMY) 

WRITE (3,1121)  EPSL NO, CQ, EPSLNS, CS 
GO  TO  5 
5100  STOP 
130  0  F OPNAT  (415) 

1)01  FORHflT(F10.3,Fl5.3) 

103  2  FORMAT  (lHf  25X  2F15.3) 

100  3  FORMAT  (1H*  6  5X  2F15.3) 

1004  FORMAT  (lHf  95X  2F15.3) 

1010  FORMAT  (2F10  •  0  ) 

1320  FORM AT (  1H0  ldX  ’QUIESCENT  DISTRIBUTIONS*  42X  ’STREAMING  OISTRIBUTI 
-IONS*//*  PITCH  ANGLE  DISTRIBUTION  NORN.  ENERGY  OISTRIBUTIO 

-N*  11X  *PITCH  ANGLE  DISTRIBUTION  NORN.  ENERGY  DISTRIBUTION* 

-  /  ?X  *  ( AlPHAQ )  •  10  X  '  <  MQ)  •  T\  *(  LOG  (E/MCJ)  5  *  2X  •  UOG(DNQ/(  OE/NC 

-Q)))*  5X  •  ( AL  PHAS )  •  10X  MHO)*  7X  • ( LOG ( E /WCS ) ) •  2X  *  (LOG  (ONS/  (OE/ 
-NCS ) > ) •  /  ) 

1)30  FORMAT  (<)0H 

1  ) 

104  0  FO°Mf T  (2F 10.0,15, 3F5,Q, 315, F5.0) 

1050  FORMAT (1M  2X  *  SURFACE  LATITUDE  (THO)  **  Fs.Zv4X  *PLASHA  PARTICLE  ( 
-POOF)  F5.0,  1  OX  *INU3R4TION  STEfS*  5X  *PITCH  ANGLE  (NALPHA),  =• 
-15/  *  EQUATORIAL  DENSITY  (XN0)  =  ’  F6.24X  ’DENSITY  MODEL  CNTYPEQ) 
-s*  15,39*  ’LATITUDE  (NX)  **  15  /  1QX  ♦  RuPAGATlON  FREQ, /EQUATORIAL 

-  CYCLOTRON  FREQ,  (OtLY)  =•  F5. 3,21X  • PERFENDICULAR  VELOCITY  ( NBETA 
-P)  *•  15  /  10X  * BANO  UNITS*  2X  *(YS)  *’  F6.3,4X  *(VE»  «•  F5.3) 

1160  FORMAT  ( 1H0  ’PRuPAG.  FREG./EQTRL  CYCLOTRON  FREC.*  5X  *NET  ANPLIMCA 
-TION  (OECIBELS)*  » 

1170  FORMAT (  15X  F 8.3,  25X  E14.6  ) 

1100  FO»MAT(lH0  *  LOCAL  B/EQTRL  B*  2X  ’REFRACTIVE*  6X  ’LOCAL  AMPLITUDE 
-EXPONENT*  ilX  ’SINE*/  1 9*  6HINOEX**2  5X  ’STREAMING*  4X  *QUIESCENT* 

-  8X  ’NFT*  5X  ’(LATITUDE)*/  8X  *(B)*  SX  *(XN2)*11X  *(XJS>’  8X  *(XJQ 
-)*  7X  *  (XNIXY)  •  6X  *  ( X)  *  ) 

1111  FORMAT  (1H0  ’PROPAG.  FREQ./EQTRL.  CYCLOTRON  FREQ.  (LOCAL)  »•  F7.4i 
1 12 1  FORMAT (1M0  ’  PLASMA  DENSITY*  7  X  ’PITCH  ANGLE’  7X  ’PLASMA  OEMS  XT Y’ 

-  YX  ’PITCH  ANGLE*  /  *  QUIESCENT/COLD’  5X  ’GLIESCLNT  NORM.’  5X  *ST 
-RFAmjng/colO’  5*  ’STREAMING  NORM.’  //  E15.4,3E20.4) 

2000  FORMAT  (312) 
i'N  3 


109. 


alJ'JPouTINE  SI  MPX  (  A  »8  »N«  FI ) 

J  iMuG°MI0N  OF  AMPLITUDt  tXPONENT  ALONG  FIELD  LINE 

■% 

IF( A . G_. l.E-5)  GO  TO  3 
AU^F  a  l.E-5 
GO  TO  4 

3  AU”F  =  A 

4  IF  (  M )  it  ,10,20 
13  FI=f 

RFTl'PN 
2  3  FN=N 

OX  r  (9-AUS£»/FN 
10*  =  :?. *0* 
ht  r  xnIxy  < AUSE) 

A I  =  C  . 
oi=r . 

*  a  AUS-r-DX 
a  4= A  USE 
nn=  r/ 2 

■jn  ?:  JsltNM 
x=<4inx 
Xfl=Xf  »TQ> 

A I  =  A  !  ♦  X  Nl  X  Y  (  X  ) 

IF ( J  t-G.MN)  oO  TU  30 
31  =  P  T ♦XNIXY <  <A) 

JO  CO’JT T HUE 

FI  =  XNlXY(B)  ♦  HI 

Fir  OX*  (FI  A*  *2.  *911/ 3. 

Rr  Tl/FN 
ENn 


no. 


FJN^riuN  aNXaYCX) 

KlAl  LAMbUA 

UlMtNiiON  HJt -jd)  ,HS  (50)  ,  XNUC5UI,  XNSC5U)  ,ALPHAQOO)  , 

5  0)  ,tMCQ(50)  ,  cNCJ(50l  ,  SPLlNSC  501  ,  SPLCHQ  C  4 , 50  *  ,SPlCHi ( 4 , 5 0*  , 
2SPlCNU(  4  ,501 

COMMON  HO,HS,XNu,XNS  ,ALPHAU, ALPHAS, fcMCU* fcHCS, SPLCHQ, SPtCHS ,SPLCNu, 

1  »PLCNa,MAO,Pi, THO , XNO , LG , CS , hPU , THOR, COSTH ,SINTH, HU), Y , Y2 , WCO 2 , 

2  B2  ,  XN2  ,u,XN,JEL,  BE TAol, BETAS  2,  POKE  ,  NUMHU, NUHHS , NUMNQ , NUMNS , 

8  NT  YPt.U»NrYPc,Nt>tTAP,NPAGt  , MONTH, NDAY,NYR 


AMPLlfUuc  cXPONtNT 

cl/ALUATt  LOCAu  FitLO  Li  Nc  PARAMETERS 


X2=X*X 

d*  sort  ci.*8.*x*>/m.-x2)**j) 

B2=3*8 

iFCNTYPtu.Nt.b)  GO  TO  oOU 
LAMBDA* ASIN(X) 
lOSl4*C  OS  ( LAMBDA  )**4 
80  OtG  -  *5286  RAUIANS 
XA»SINCLMMaOA-.528t»)  **2 
30  U  HP2*WPU*HPU*UX2  , XA, COSTH, COS L4*NTYPEUI 
IF  (PORc-lO.)l,2,2 

1  Z=Y 

nO  TO  J 

2  1-6 

8  XN2*l.t  HP2/CHC02*Z*(d-m  *PORt _ 


LlTHiUM  SHAPtU  LHARG t  INJECTION  '  " 

^ASt  5  For  lithium  Injection  only 

IFCX2.GI  .U.0Q8)  GO  TO  5  .  rrr  n  c  n  cn  .  .  , 

iFCY.LT. 6/7. I  GU  TO  4  atVDt  =  °'5  Re  (t  =  50  mlnutes) 

xnixy=o.  Delete  for  other  models. 

RETURN 

4  xN2  =  XN2  +  ((1.06MlXP(-(X2/<.1o*C0STH**4*  »  )  >  >  •  (RP0/HCQ  )  **Zi  *  1 7  .♦PORE 

_  -  >/C7.*Y*(8-7.*Y) ) 

5  XN=  SORT  (XN2 ) 

INTcGRATE  BtT A  PEfiPcNJi CUe AR  OVER  DISTRIBUTIONS 
ocT APM=  SORT  U.- V2/82  +  Y2/ (82* C i. +62/ CXN2*Y2) ) I) 

NTYPE=8 

IF ( POKl.oT . 10. i  GO  TO  8 
CAwL  SIMFC0* ,dtTAPM, NdcTAP,XJQ  > 

GO  TO  9 

8  CAlL  SiMPCU.  ,.0l,NdtTAP,XJUl) 

CAll  SiNPC.Ul  ,Bl  IAPM,NdtTAP,XJQ2i 
XJU=XJwi+XJu2 

H  IFC  NUNNS. oT.  0)  GO  TO  20 
10  X J8  -  0. 


111. 


gu  r  o  o  u 

iu  Ucf  AP1=  di_T  AN  (dfc.  T  ASl  ,d,¥  ,XN2I 
dcf  AP2  =  olI  AN  (del  AS2,d,  i ,  XN2)‘ 

1F( dtrAPl.tM.-l. .OP. dcTAP2.tQ .-1 »J  GOTO  10 
NT\rPt=4 

oAll  ii mp (dt tapi ,0cTap2,nbcTap,xjs) 
iw  L.ilK  Y=*»  .  <m*vtPt/  (XN2*T2*rtCQ2)  *  (Xji^XdU) 

.vni  a  ys2  Nlxr 

r  » lujtii  d,  an2  ,  xji  #xju»^nixv  ,x 

louu  FuP.ImT  (  (.Id.  A  ,clfa.A,  3cl3*‘*,F10 .4) 

r.L.  I  J  M 
CNJ 


1 1 2 


S U '3  F OUT  I N Z  S  .1  H P  *  A  ,  o  ,  N  ,  F  i  ) 

I.* 

G  INTEGRATION  OVER  9ET  A  PERPENDICULAR  FOR  AMPLITUDE  EXPONENT 

G 

IF(A.GE.i.E-7)  GO  TO  3 
AUSE  =  l.E-5 
GO  TO  4 

3  AUSE  =  A 

4  IF  (N)  10,10,20 
10  FI  =  C 

RETURN 

20  IF  (ABS(A)  +  ABS  (  0)  )  30,10,30 
3  0  FN=N 

DX  =  ( B-AUSE) /FN 
TDX=  2  #  *  OX 

FI  =  GRANF(AUSE) +GRANFO) 

AIS  P . 

QI  =  0- 

X  =  AUSE-OX 
NN=N/2 

DO  40  J=1 , NN 
X=X+TOX 

40  AI  =  A  T  +GRANF  (  X  ) 

X  =  AUSE 
NM=NN- 1 
DO  5"  J - 1 , N M 
X=  X  +  T  D  X 

50  BI=PI +GRANF ( X ) 

FI=PX*< FI +4,  *AJ>2,*Bl)/3. 

RETURN 

END 


w  u  c>  (.<».<  O  () 


113. 


FUNCTION  GRANF(X) 

OIMFNSION  HQ (50) ,HS(50) »XNQ(50) , XNS « 50 ) , ALPHAQ ( 50 ) ♦ 

1ALPHAS(5C>  ,EHCQ<5  0> ,wMCS(50)  ,SPLCNS(5U) , SPLCHQ ( 4 ,50 ) ,SPLCHS(4 ,50) , 
2SPLCNQ ( 4,50) 

COMMON  HQ.HS',  XNQ ,  XNS ,  AL  PH  AQ  ,  AL  PH  AS  ,EMCQ,  E  MCS ,  SPLCHQ,SPLCHS  ,SPLCNQ , 

1  SPLCNS,RAD,PI,TH0, XN0 , CQ, CS , WP0, TH0 R , CD STH ,  SINTH , WCC  »  Y , Y2, WC  02* 

2  9?,XN2,B,XN,OEG,  BETAS1 , BETAS2 , PORE ,NUMHQ  ,  NUUHS ,NUMNQ , NUMNS, 

3  NTYPEQ,NTYPE,NB£TAP,NPAG£,HONTH,NOAY,NYR 


EVALUATION  OF  VARIOUS  INTEGRALS 

GO  TO  (10, 23  , ‘,0,40,180,190)  ,NTYPE 

N TYPE  1  OR  2  PITCH  ANGLE  OISTRIB  INTEGRAXD 
10  CALL  SPLINE (ALPHAQ, HQ, NUMHQ,SPLCHQ,X* DEG, YY, OX) 
GO  TO  30 

20  CALL  SPLINE( ALPHAS, HS ,NUMMS , SPL CHS , X* DEG , YY , OX) 
30  GRANF=YY*  SIN(X) 

RETUDN 

N  TYPE  ?  OR  h  AMPLITUDE  EXPONENT'  INTEGRAND 
43  n=Nt YPE "2 
X2  «  X*X 

Cl=l  -<1.+B2/(XN2*Y2))M1,-(1.-X2)/(Y2/E2)  ) 
IF(C1,LT.0.)  C1=0. 

QET/IP*(1,-  SQRT  (Cl)  )/(XNMl«+B2/  (XN2*  Y2) )  ) 

BET  A  =  SQRT  (  X2  +  BET  AR**2) 

Cl=l .-BETA**2 
IF  (Cl) 90,103 ,100 

90  WRTTF  (3,1010) X, Cl, BtTA.BETAR, 32, XN2 
<•  C1  =  -C1 

100  GAMMA= 1 , /  SQRT ( Cl ) 

ALPHAS  ATAN(-XZBETAR) 

ALPHA  s  ABS(ALPHA) 

XXs (GAMMA-1.  )  *PORE 
IF,(X«GT.0.)  GO  TO  120 
WRITF  (  3,1020  ) X, Cl, GAMMA, XX 

GRANFsQ.O 

RETL'PN 

120  XX  =  ALOGIO(XX) 

BSQP= 1 • /  SQRT ( 9) 

U=BS  OR*  SIN( ALPHA  ) 

C  Cl=  ASIN(U) 

GO  TO  (  130, 140)  , N 

130  CALL  SPLINE  (  EMCQ , XNQ,NUMNQ, SPL CNQ ,XX , Y Y  ,  ZZ) 
IF(YY.GT.-10. )  GO  TO  133 
GRANFs  q.o 
RETURN 

1  133  CALL  SPLlNE(ALPHAQ,HQ,NUMMQ,SPLCHQ,Cl*OEG,H,OH) 

134  C=CO 

GO  TO  150 

140  CALL  SPLINE  ( EMCS  ,XNS ,NUMNS , SPLCNS , XX  ,  Y Y ,  ZZ) 

CALL  SPLINE ( ALPHA S ,HS  ,NUMMS , SPL CHS, Cl* DEG, H,CH) 

c=cs 


C  J  o 


114. 


150  G=(ir’.**YY)/<XN0*GAMMA**2*BETA>*  PORE 
IF  <C1*0EG-1J .) 160,170, 170 
160  H=0. 

DH=T 

170  F=G*C*H 

DF=G*C*0H*DEG 

GAMMA4=GAMMA**4 

GRANF  =  GAMMA4* (2. *X-X2*X) *F-GAMMA  4* ( i.-B/ (GAMMA* Y) ) *X2/BETAR*DF 
RETURN 

N TYPE  5  OR  6  ENERGY  0ISTRI8  INTEGRAND 
180  GAMMC=1./  SQRT ( 1. -X* X) 

XX  =  ALOGIO ( (GAMMA-1 .)*PORE) 

CALI  SPLINE (EMCQ,XNQ,NUMNQ,SPLCNQ, XX, YY, OH) 

GO  TO  200 

190  GAMMA  =  1./  SQRT  ( l.-X*X) 

XX  s  ALOG10 ( (GAMMA-1 . )*PORE> 

CALI  SPLINE (EMCS,XNS,NUMNS,SPLCNS, XX, YY, OH) 

20  0  GRANFs  10.**YY  *  GAMMA**3  *  X  /  XNO*  PORE 
RETURN 

1  30  0  FORMAT  (2H015E15.6) 
niO  FORMAT  (2H026E15.6) 

132  0  FORMAT  (2HQ34E15.6) 

ENO 


Simrui  INE  SPLICN(X,  Y  ,M,C) 

DIMr  MS I  ON  X(50),Y(5C) ,0(50 ,P(50) ,E(50) ,9(50) ,Z(50) , 

-  CCt  ,f-0),A(50,3) 

CURVE  FTT  COEFFICIENT  ROUTINE 


00  1^  K=1,MM 
0(K) =X (K+l)-X (K) 

P(K) =0(K)/6. 

IP  E  (K )  =  ( Y  (Kfl)-Y(K)  )/0(K) 

00  ?r  K=2,MM 

20  9<K)  =t(K)-t(K-l) 

A(l,2)=-1.-0(1)/D(2) 

A (1 «  ? ) s0 (1) /O ( 2 ) 

A(?,7)=P(2)-P(l)*A(l,3) 

A(?,?)=2.*(P(i)+P(2))-P(l)*A<l,2) 

A(?,3)=A(2,3)/A<2,2) 

B(2)=B(2)/A<2,2) 

00  3 r  K--3,MM 

A(K.?)*2.MP(K-1)  +P  (  K ) )  -P  (  K**l)  * A (K-l »  3 ) 
9  (K) -B(K)-P(K“1) *  9  ( K-l) 

A (K »  3) =P(K) /A ( K , 2 ) 

Q=0(M-2)/0(M-l) 

30  B(<) -B(K)/A(K,2) 

A(Mtl)«l.*Q+A(M-2,3) 

A(M,?)=-Q-A(M,l)*A(M-lt3) 

B('i)=B(M-2)-A(M»l)*B(M-l) 

Z(M) =B(M)/A(H,2) 

MN=  N- 2 

DO  U'  I =1 1 MN 
K=^-T 

40  Z(<)sB(X)-A(Kf3)*Z(K+l) 

Z(1)=-A(1,2)*Z(2J-A(1,3)»Z(3I 
00  Fr  K=1,MM 
0=1 ./(6.*0(K) ) 

C(1,K)=Z(K)  *Q 

C(2«K)>Z(K*1)  *Q 

C(3 f K)>Y(K)/0(K) -Z(K)*PCK) 

50  C(4 ,K) =Y(K*1) /0<K)-Z(K*1) *P(K) 

RETURN 

END 
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SUBROUTINE  SPLINE ( X , Y ,M , 3 , X INT  ,  YINT,  OT  ) 

DIMENSION  X ( 5  0 ) *  Y  (50)  *0(4*50) 

0 

C  CUR VF  FIT  ROUTINE 

C  THIS  SUBROUTINE  USES  SPLINE  COEFFICIENTS  FROM  SPLICON 

■> 

:  FIND  WHAT  TWO  VALUES  OF  XINT  IS  BETWEEN 

*% 

'U 

IF  ( X  INT  «GE»  X  ( 1)  )  GO  TO  50 
TINT = Y ( 1 ) 

OT  =  C  « 

RETURN 

50  IF  (  XINT.LE.  X  (  M)  )  GO  TO  60 
YINT=Y (M) 

OT  =  0 . 

RETURN 

6  C  OO  1?C  J=2,M 

IF(XINT#LE»X(J-1) .OR.XlNT.GT.X(Ji)  GO  TO  ICO 
MU=J 
ML=J-1 
GO  TO  110 
10  0  CONTINUE 

WRITE ( 3  »10  00) 

STOP 

*■> 

C  SOLVE  FUNC1ION  FOR  OEPENOENT  VARIABLE  ( YINT) 

3  OT  IS  THE  DERIVATIVE  OF  Y  =.F(XINT) 

C 

110  Xl=  (X  (MU)-XINT)  **2 
X2=  (XINT-X(ML) ) **  2 

YINT  =  (X(MU)-XINT)  ♦  <C  ( 1,  ML )  * X 1 +  C  (  3 ,  ML)  ) 
YINT=YINT+(XINT-X (ML) )» (C(2,ML)+X2*C(4«ML> ) 

12  0  DT  =  -3.*C(l,ML)*XH-3.*C(2,ML)*X2-C<3tML>  +C  (4»ML> 

RETURN 

130  0  FORMAT  (13H10UT  OF  RANGEE15.6) 

ENO 


I 
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FUNCTION  CI(a2,XA,C0STH,0JSL4,NTYPEQ> 

C  CL  0  PLASMA  DENSITY  MOOELS 
REAL  LAMBDA 

GO  TO  (10, 20, 30, 40*50,  fail,  70,  dO>  ,NTYPEQ 
8-FIELT  MOOEl 
10  Q=l./((1.-X2)**3> 

RETURN 

H  V  ORQST  AT  I C  EQUILIB  MOOEL 

20  Q=<  EXP(X2/(l.-X2)))**(3.*COSTH**2) 

RETURN 

ECTRL  INJECTION  MOOEL 

3  0  Q=1 ./ < (1.-X2) **3> ♦68.0* (EXP  <-<X2/< .f l*COS Th**4> ) ) ) 

RETURN 

4  0  Q*l./  ( (1.-X2) **3) ♦  &. 55* (EXP(-(X2/(.04*COSTH**4) ) ) ) 

RETURN 

50  Q=l./((1,-X2)**3) *1. G6* (EXP<-(X2/< .16*COSTH**4) ) ) ) 

RETURN 

NONE  OUA  TORI AL  BURST 

6  3  U*l./  ((l--X2)**3)+8«55* ( * XP ( -C0SL4*XA / ( , 0 4*CCSTH**4> ) ) 
RET  UDN 

HARIUM  SHAPED  CHARGE  MOOEL  -  FIELO  ALIGNED  STEP  FUNCTION 
3n  OrG  -  x  -  36  OEG  SET  XN0*1. 

70  Q=1  ,  /  (  ( 1  -X  2 )  **  3 ) 

IF<  X2.GE..25.ANO.  X2  .  LE.  ,  3  7  9  0*,)  Q=Q+30. 

RETURN 

3ARIUM  SHAPED  CHARGE  MODEL  -  FIELO  ALIGNED  STEP  FUNCTION 
6  OEG  -  X  -  25  OEG  SET  XN0=1. 

50  Q=l./(<1.-X2>**3> 

IF(X?.GE..O1037.ANQ.X2.Lt.. 17361 )  Q*Q+5. 

return 

.  END 


CiOtJOU 


11 


<>. 


FUNCTION  BETAN(BETA,B,Y,XN2) 

'V 

'  INTEGRATION  LIMITS  FOR  BETA  PERPENDICULAR 

C 

BSM2=8ETA**2 
Ci=  SORT (1.-BSM2) 

FACT=  9SM2-(l.-Cl*B/Y)**2/XN2 

IF( FACT. LT.O.  .OR. FACT. GE. 1. )  GO  TO  5 

BET  A  NsSQRT (FACT) 

RETURN 
5  BETANs-1. 

RETURN 

END 


SUBROUTINE  SPRINT 

OIMFNSION  HQ (50) ,HS(5Q) ,XNQ(50) .XNS (50) *ALPHAQ(50) , 

1 ALPHAS (50 ,EMCQ(50> »EMCS(50) ,SPLCNS(50) ,SPLCHG(4,50) ,SPLCHS(4 ,50)  , 
2SPLCNQ  (4,50)  *  * 

COMMON  HQ.HS tXNQ, XNS .ALPHA Q, ALPHAS .1 MCQ » E MCS. SPLCHQ .SPLCHS .SPLCNQ » 

1  SPLCNStRADtPItTHOfXNOfCQ.CSfWPO.THOR.CCSTH.SINTH.NCO.YtY2.WC02. 

2  92.XN2.B.XN. DEG. BETAS1.BE TAS2. PORE, NUMHQ , NUMHS »NUMNQ, NUMNS. 

3  NTYPEQ.)  fyPE.NBETAP.NPAGE , MONTH, NC AY ,NYR 


OATE, MONTH, YEAR  ANO  PAGE  NO 

WRITr  (3,1000) MONTH, NDAY.NYR.NPAGE 
NPAGE=NPAGE*1 
RETURN 

130  0  FORMAT  ( 1H11 2 , 1H/ 12 , 1H/ 12 , 4 OX 8HPAGE  NO. 14/) 

END 


